\a 

o- 


i  M  ! 

J.  M I  ?  P 


WL  OF  SPACE  STRUCTURES I 


The  Charles  Stark  Draper  Laboratory,  Inc. 


Sponsored  fey 

DeHFonse  Advanced  Research  Projects  Agency  (DOD) 
ARPA  Order  No.  365S 


APPROVE  FOR  PUBLIC  RELEASE;  DISTRIBUTIOH  UNLIMITED 


The  views  and  conclusions  contained  In  this  document  are  those 
or  the  authors  and  should  not  be  Interpreted  as  necessarily 
representing  the  official  policies,  either  expressed  or  Implied, 
of  the  Defense  Advanced  Research  Projects  Agency  or  the  U.S. 
Government. 


o~ 

CD 

CD 

UD 
_ J 

Lu 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Griff iss  Air  Force  Base,  NY  13441 


hr  J’eeni  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and  is 
v  if  1°  16  ational  Technical  Information  Service  (NT1S) .  At  NTIS  it  will 

le  "0  the  general  public,  including  foreign  nations. 


•-TR-82-.1  j.'L 


V°1  II  (cf  two)  has  been  reviewed  and  is  approved  for  publicatior 


APi’RC.H 


RICHARD  CARMAN 
Project  Engineer 


'  P  ■  ■  / 

APPROVED :  C.  _  4^;,  OU/r 

LOUIS  E.  WHITE,  Lt  Colonel,  USAF 
Acting  Chief,  Surveillance  Division 


FOR  THE  COMMANDER: 


/? 


JOHN  P.  HUSS 


Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  Is  no  longer  employed  by  your  organization, 
please  notify  RADC  (0C3E  )  Griffiss  APB  NT  13441.  This  will  assist  us  in  M 
maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  rep  't  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  be  returned. 


ACOSS  ELEVEN  (ACTIVE  CONTROL  OF  SPACE 
STRUCTURES),  Volume  II 


Thomas  H,  Brooks 
Virenda  N.  Mahajan 
Jiguan  G.  Lin 
Glen  J.  Kissell 
Daniel  R.  Hegg 
Ronald  K,  Pearson 
James  D.  Turner 
Hon  M.  Chun 
Eliezer  Foger 


Contractor:  The  Charles  Stark  Draper  Laboratory,  Inc. 

Contract  Number:  F30602-81-C-0180 

Effective  date  of  Contract:  27  April  1981 

Contract  Expiration  Date:  27  April  1984 

Short  Title  of  Work:  ACOSS  Eleven 

Program  Code  Lumber:  1E20 

Period  of  Work  Covered:  May  1981  -  November  1981 

Principal  Investigator:  Dr.  Keto  Soosaar 
Phone:  617-258-2575 


Project  Engineer: 

Phone: 


Richard  Carman 
315-330-3148 


Approved  for  public  release;  distribution  unlimited 


This  research  was  supported  by  the  Defense  Advanced  Research 
Projects  Agency  of  the  Department  of  Defense  and  was  monitored 
by  Richard  W.  Carman  (OCSE) ,  Griffiss  AFB  NY  13441  under 
Contract  F30602-81-C-0180. 


UNCLASSIFIED 


WOlg ,TY  CLASSIFICATION  OF  THIS  PACK  fWl—  Ow.gnltftO, 


REPORT  DOCUMENTATION  PAGE 


V.  REPORT  NUmIZX 

RADC-TR-82-131.  Vol  II  (of  two) 


2.  GOVT  ACCESSION  HO, 


«.  TITLE  f«*W  SuMllliJ 

ACOSS  ELEVEN  (ACTIVE  CONTROL  OF  SPACE 
STRUCTURES) 


S.  TYPE  OF  REPORT  4  PERIOD  COVI  EO 

(Interim  Report 
May  1983  -  November  1981) 


?.  AUTHORfaJ 

Thomas  H.  Brooks 
Virenda  N.  Mahajan 
Jiguan  G.  Lin 


(.  PERFORMING  OTO.  REPORT  NUMBER 

|cSDL-R-1536 


Glenn  J.  Kissell 
Daniel  R.  Hegg 
Ronald  K.  Pearson  (Cont’dj) 


».  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

The  Charles  Stark  Draper  Laboratory,  Inc. 
555  Technology  Square 
Cambridge  MA  02139 


It.  CONTROLLING  OFFICE  NAME  ANO  ADORES* 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Blvd 

VA  22209 


.Arlington 

I*.  MONITORINC 


MONITORING  AGENCY  NAME  4  AOORESyif  dlUtnnl  tnm  Controlling  O/HcaJ 

Rome  Air  Development  Center  (OCSE) 

Griffiss  AFB  NY  13441 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT**  CATALOG  NUMBER 


4.  CONTRACT  OR  GRANT  NUMBER^; 

F3G602-81-C-0180 


1 


10.  PROGRAM  ELEMENT.  PROJECT,  1  »«. 
AREA  4  WORK  UNIT  NUMBERS 


62301E 

C6550104 


II.  REPORT  DATE 


May  1982 


II.  NUMBER  OF  PAGES 

IL44 


IS.  SECURITY  CLASS,  (ol  thl a  rt, 

Unclassified 


isa.  DECL  ASSI fTcAt!on7 Ooii  '.  IsOING 
kj/A 


schedule 


l«.  OISTHIBUTION  STATEMENT  (o!  IMa  Rapo/I) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (at  lha  aAalracf  aniararf  In  Block  30,  II  d  .sianr  from  Ropotl) 


Same 


14.  SUPPLEMENTARY  NOTES 

RADC  Project  Engineer:  Richard  Carina,  (OCSE) 


1*.  KEY  WOAOS  (Conttou#  on  tarot  a  a  tida  it  nacaaaary  and  IdantHy  by  bloc*  nwnfctr) 


Simulation 

Draper  Integrated  Simulations 
Surveillance  Systems 
HALO  Optics 


x  HAL' 


Deformable  Mirror 
Order  Reduction 
Design  Freedom 
Optimal  Control 


Large  Angle  Slew 
Maneuver 


ABSTRACT  fConffnuo  on  ravataa  *ldo  tt  nacaaaary  and  IdantHy  by  blash  m*nb# t) 

The  Active  Control  of  Space  Structures  (ACOSS)  Eleven  contract  includes 
three  distinct  areas  of  endeavor:  Simulations  Extension,  HALO  Optics  and 
ACOSS.  This  report  covers  the  work  performed  in  each  of  these  areas 
between  May  and  November  1981.  Simulations  Extension  and  HALO  Optics  are 
presented  in  Volume  I  and  ACOSS  is  presented  in  Volume  I] 


DD  , 


FORM 
JAN  73 


1473 


EDITION  OF  I  NOV  «8  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (tflian  Data  Entorod) 


UNCLASSIFIED 


security  classification  or  this  wkihw  cm  tntmo 


Item- 7  -  Authors  (Cont’d) 

James  D.  Turner 
Hon  M.  Chun 
Eliezer  Foger 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OR  Tuir  RAOEOWiwi  D*«  SnttnMQ 


ACKNOWLEDGMENT 


This  report  was  prepared  by  The  Charles  Stark  Draper  Laboratory,  lac., 
under  Contract  F30602-81-C-018Q.  This  research  was  b  ipported  by  ihe  Advanced 
Research  Projects  Agency  of  Che  Department  of  Defense  and  monitored  by  the 
Rome  Air  Development  Center. 


The  program  manager  is  Keto  Soouaar  and  the  project  leaders  are  Mr. 
Thomas  H.  Brooks  for  Simula*"^;5*  Extension,  Dr.  Virendra  N.  Mahajan  for  HALO 
Optics,  and  Mr.  Robert  R.  Str>  t.ce  for  Active  Control  of  Space  Structures 
(ACOSS).  The  technical  director  for  ACOSS  is  Dr.  Jiguan  G.  TAn.  The  authors 
of  this  report  are:  Mr.  Thomas  H.  Brooks  (Section  2),  Dr.  Virendra  N.  Mahajan 
(Section  3),  Dr.  Jiguan  C.  Lin  and  Mr,  Glen  J.  Kissel  (Section  4),  Dr.  Jiguan 
G.  Lin  (Section  5),  Dr.  Daniel  R.  Hegg  (Section- 6) ,  Dr.  Ronald  K.  Pearson 
(Section  7),  Dr.  James  D.  Turner  and  Mr,  Hon  M.  Chun  (Section  8),  and  Dr. 
Eliezer  Fogel  (Section  9).  Assistance  from  Dr.  Yu-Hwan  Lin  (Section  4)  and 
Mr.  Warren  J.  Jasper  (Section  5)  is  acknowledged  gratefully. 

Publication  of  this  report  does  not  constitute  approval  by  the  Defense 
Advanced  Research  Projects  Agency  or  the  U.S.  Government  of  the  findings  or 
conclusions  contained  herein.  It  is  published  for  the  exchange  and  stimula¬ 
tion  of  ideas. 


Accession  Vor_ 
jty~v 

r  •;  TAB 

^  .-yvv'VUnGOd 

j  ideation-- 


D\<:irlbuVi«o/ _ _ 

J  'Avf.il  *j*i/** 
Inlrt  \  Special 


TABLE  OF  CONTENTS 


Section  page 

4  COMPENSATED  TRUNCATION  OF  MODAL  MODELS  FOR  DESIGN 

OF  CONTROL  SYSTEMS . . . . .  1 

4.1  Introduction... . 1 

4.2  Reduced-Order  Modeling  for  Control-System-  - 

Designs  by  Modal  Truncation. • . . . . .  1 

4.3  Compensation  for  Truncation .  2 

4.4  A  Completion  of  the  Synthesis .  3 

4.5  Filters  for  Spillover  Reduction . 6 

4.6  Conclusion... . 7 

LIST  OF  REFERENCES . 8 

5  ENSURING  FULL-ORDER  CLOSED-LOOP  STABILITY  IN  THE 

REDUCED-ORDER- DESIGN  OF  OUTPUT  FEEDBACK  CONTROLLERS . . .  9 

5.1  Introduction .  9 

5 .2  Problem  Formulation . . . . . . . . .  3  9 

5.3  Full-Order  Closed-Loop  Asymptotic  Stability  Conditions . 12 

5.4  Principle  of  the  Reduced-Oroer  Design  Technique .  12 

5.5  Discussions . 13 

LIST  OF  REFERENCES .  14 

6  DESIGN  FREEDOM  AND  THE  IMPLEMENTATION  OF 

SUBOPTIMAL  OUTPUT  FEEDBACK  CONTROL .  15 

6.1  Introduction . 15 

6.2  Design  Freedom. . 15 

6.3  Case  Study:  Numerical  Implementation  of 

Suboptimal  Output  Feedback  Control . . . . .  18 

6.4  Some  Facts  from  Linear  Algebra. .  20 

6.5  Conclusion . . . , . . . .  27 

6.6  Addenda... . . . 28 

LIST  OF  REFERENCES . 29 

7  STOCHASTIC  OUTPUT  FEEDBACK  COMPENSATORS  FOR 

DISTRIBUTED  PARAMETER  STRUCTURAL  MODELS .  34 

7.1  Introduction, . 34 

7.2  Notion  of  Weak  Convergence  and  Weak  Compactness . 35 

7.3  General  Existence  and  Convergence  Results . 37 

7 .4  Application  to  Fixed-Form  Compensator  Design  Problem .  48 

7.5  Summary  of  Convergence  Conditions. .  60 

7.6  Conclusions . . . .  . .  64 

LIST  OF  REFERENCES . .  ...  65 

8  LARGE- ANGLE  SPACECRAFT  SLEWING  MANEUVERS . 66 

8.1  Introduction. . . . 66 

8.2  Free  Final  Angle  Optimal  Slewing  Maneuvers... .  66 


iii 


TABLE  OF  CONTENTS  (cont.) 


Section  Page 

8.3  Control-Rate  Penalty  Technique  for  Optimal 

Slewing  Maneuvers. . . . 76 

8.4  Single  Axis  Retargeting  Maneuvers  for 

a  Rigid  Spacecraft......... . 89 

8.5  Preliminary  Study  of  Slewing  Maneuvers  for 

the  ACOSS  Model  2 . . .  99 

8.6  Example  Maneuvers . . . 99 

LIST  OF  REFERENCES . i08 

9  ORDER  REDUCTION  BY  IDENTIFICATION— SOME  ANALYTICAL  RESULTS . 114 

9 . 1  Introduction,, . . 114 

9.2  Least-Squares  Identification . ..115 

9.3  Second-Order  System......... . 121 

9.4  nth  Order  Systems . ....124 

9 . 5  Conclusions . . 127 

LIST  OF  REFERENCES . .,..128 


iv 


LIST  OF  ACRONYMS 


AVD  Air-Vehicle  Detection 

CSDL  The  Charles  Stark  Draper  Laboratory ,  Inc. 

DARPA  Defense  Advanced  Research  Projects  Agency 

DIS  Draper  Integrated  Simulations 

DISCOS  Dynamic  Interaction  Simulation  of  Structures 
DMA  Defense  Mapping  Agency 

FFT  Fast  Fourier  Transform 

GSS  Generic  Scene  Simulation 

HAC  Hughes  Aircraft  Corporation 

IFOV  Instantaneous  Field  of  View 

IMSL  International  Mathematical  and  Statistical  Libraries 

IR  Infrared 

LOS  Line  of  Sight 

LS  Least  Squares 

lsc  Lower  Semicontinuous  Function 

LSS  Large  Space  Structures 

OTF  Optical  Transfer  Functions 

PRA  Photon  Research  Associates 

PSF  Point-Spread  Function 

TPBVP  Two-Point  Boundary  Value  Problem 

wise  Weakly  Lower  Semicontinuous  Function 

wsisc  Weakly  Sequentially  Lower  Semicontinuous  Function 


v 


•MBWBJBSJWo K 


1 


PREFACE 


Volute  2  of  this  report  describes  the  work  CSDL  has  done  to  investigate 
spacecraft  control  theory.  Each  of  the  six  sections  devoted  to  ACOSS  reports 
on  a  different  aspect  of  that  work. 

Section  4,  "Compensated  Truncation  of  Modal  Models  for  Design  of  Control 
Systems i"  describes  the  selection  process  necessary  in  large  space-  structure 
(LSS)  control-system  design  using  a  truncated  finite-element  model.  The  trun¬ 
cated  model  must  be  selected  properly  and  compensated  explicitly  for  control 
and  observation  spillover,  so  the  control  system  designed  through  this  method 
can  perform  satisfactorily  when  implemented  on  the  structure.  Proper  selec¬ 
tion  requires  correct  classification  of  structural  modes  in*o  "nrimary"  and 
"secondary"  modes.  Explicit  compensation  for  truncation  includes:  placement 
of  actuators  and  sensors,  synthesis  of  the  actuator  and  sensor  influences  once 
they  are  placed  on  the  structure,  and  filtering  of  the  actuator  inputs  and 
sensor  output*. 

Section  5,  "Ensuring  Full-Order  Cloeed-Loop  Stability  in  the  Reduced- 
Order  Design  of  Output  Feedback  Controllers,"  builds  on  the  studies  performed 
during  ACOSS  6  that  established  various  conditions  necessary  to  ensure  full- 
order  closed-loop  asymptotic  stability  and  robustness  with  reduced^-order  de¬ 
sign  of  velocity  and  displacement  output  feedback  controllers.  Currently,  the 
work  in  this  area  concentrates  on  how  to  apply  such  results,  to  large  flexible 
space  structures  and  how  to  develop  a  reduced-order  design  technique  that  will 
ensure  full-order  closed-loop  asymptotic  stability. 

The  study  includes  preliminary  development  of  computer-aided  design 
software  and  acceleration  output  feedback  control. 

Section  6,  "Design  Freedom  and  the  Implementation  of  Suboptimal  Output 
Feedback  Control,"  discusses  the  freedom  inherent  ii.i  design.  The  section 
states  that  often  this  freedom  is  sacrificed  purposely  when  simplifying  as¬ 
sumptions  are  made  to  avoid  theoretical  or  computational  difficulties.  Since 
it  is  difficult  to  consider  this  topic  without  referring  to  specific  applica¬ 
tions,  the  section  uses  controller  design  as  an  example  where  work  is  being 
done  to  discover  and  exploit  the  freedom  of  choice  in  design.  Then,  the  sec¬ 
tion  uses  suboptimal  output  feedback  control  as  a  case  study  which  is  relevant 
to  ACOSS  development. 

Section  7,  "Stochastic  Output  Feedback  Compensators  for  Distributed 
Parameter  Structural  Models,"  presents  recent  progress  on  the  stochastic 
output  feedback  design  problem  for  distributed  parameter  plants.  The  results 
presented  are  an  extension  of  work  done  under  the  previous  contract. 

The  concepts  developed  are  general  enough  to  apply  to  a  wide  variety  of 
fixed-form  coapensator  design  problems,  and  current  studies  are  aimed  at 
specializing  the  results  to  the.  optimal  output  feedback  compensator  design 
problem.  The  procedure  developed  will  be  applied  to  the  design  of  velocity 
feedback  controllers  for  a  vibrating  string.  The  results  of  this  simple  test 
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should  provide  insight  into  the  impact  of  various  modeling  assumptions  oh  the 
convergence  of  the  design  procedure  described. 

Section  8,  "Large-Angle  Spacecraft  Slewing  Maneuvers, ”  further  develops 
work  that  was  reported  in  the  previous  ACOSS  contract.  Specifically,  the  sec¬ 
tion  presents  techniques  for  improving  the  optimal  torque  profiles  by  allowing 
the  solution  process  to  determine  the  optimal  terminal  boundary  conditions  and 
by  developing  a  control-rate  penalty  technique  for  producing  smooth  control 
profiles.  Several  example  maneuvers  are  provided  to  demonstrate  the  practical 
application  and  utility,  of  the  techniques  presented. 

Section  9,  "Order  Reduction  by  Identification — Some  Analytical  Re¬ 
sults,"  attempts  to  characterize  control  designs  that  will  guarantee  stability 
using  a  reduced-order  model.  This  lind  of  design  compromise  is  practiced 
regularly,  but  no  one  has  verified  the  validity  of  such  an  approach. 

The  least  squares  (LS)  method  is  used  in  this  analysis  because  it  is  a 
relatively  robust  identification  scheme  and  analytical  expressions  for  order 
reduction  already  exist  for  it.  The  results  of  the  analysis  show  that  a  re¬ 
duced  order  controller  can  be  built  using  the  LS  method  of  identification.  It 
is  planned  to  demonstrate  the  practicality  of  this  approach  on  Draper  Model  t‘2 
in  the  near  future. 


SECTION  4 

COMPENSATED  TRUNCATION  OF  MODAL  MODELS 
FOR  DESIGN  OF  CONTROL  SYSTEMS 


4. 1  Introduction 


Finite-element  models  of  a  large  space  structure  are  too  large  for  the 
design  of  its  control  pystems,  let  alone  implementation  of  "modem"  control 
systems  in  space.  The  model  must  be  truncated  substantially  to  reduce  it  to 
a  reasonably  low  order.  A  proper  truncation  requires  proper  structural-mode 
classifications  according  to  their  influences  on  required  performance  and  sen¬ 
sor  outputs  and  their  responses  to  probable  disturbances  and  actuator  inputs. 

The  truncated  modal  models  must  explicitly  be  compensated  for  their  truncation 
to  prevent  or  reduce  control  spillover  and  observation  spillover. 

4 . 2  Reduced-Order  Modeling  for  Control-System  Designs  by  Modal  Truncation 

In  practice,  not  all  structural  modes  that  are  "molded"  by  a  large 
finite-element  model  are  to  be  calculated  because  of  the  increasing  computa¬ 
tional  expense  and  inaccuracy.  Yet,  the  number  of  modes  commonly  calculated 
(using  NASTRAN,  for  example)  are  very  large  for  a  typical  large  flexible  space 
structure.  Only  a  very  limited  number  of  these  calculated  modes  can  be  used 
to  form  a  reduced-order  model  required  for  design  of  the  structure's  control 
systems,  however.  Some  ideas  for  appropriate  reduced-order  modeling  are  being 
formulated  here  at  CSDL.  The  following  (here  and  Section  4.3)  is  a  preliminary 
sketch. 

Among  the  calculated  modes,  some  are  to  be  classified  as  "primary  modes" 
and  some  as  "secondary  modes".  Those  modes  that  will  influence  the  specified 
performance  of  the  structure  significantly  (e.g. ,  the  line-of-sight  error  and 
defocus)  and/or  will  be  influenced  significantly  by  probable  disturbances  on 
the  structure  (e.g.,  initial  disturbances  due  to  maneuvers,  sinusoidal  or  ran¬ 
dom  disturbances  caused  by  on-board  equipments)  are  classified  as  primary  modes. 
Some  primary  modes  are  regarded  as  critical  if  a  certain  critical  level  is 
exceeded.  With  respect  to  a  configuration  of  actuators  and  sensors  placed  on 
the  structure,  secondary  modes  are  those  nonprimary  modes  which  either  can  be 
influenced  significantly  by  the  actuation  £r  can  influence  the  sensing  signifi¬ 
cantly.  Useful  ranking  techniques  are  being  formulated  and  studied. 

Naturally,  those  calculated  modes  which  will  not  only  be  influenced  sig¬ 
nificantly  by  the  probable  disturbances  but  also  influence  the  specified  per¬ 
formance  significantly  should  be  "modeled"  (i.e.,  retained  in  the  reduced-order 
model)  for  design  of  control  systems.  If  possible,  all  other  primary  modes 
should  be  modeled  as  well.  On  the  other  hand,  it  is  obvious  intuitively  that 
any  calculated  mode,  be  it  primary  or  secondary,  also  should  be  modeled  for  de¬ 
sign  of  control  systems  if  it  can  be  influenced  strongly  by  the  actuation  or  if 
it  can  influence  the  sensing  strongly.  A  recent  research  into  the  linear- 
quadratic  regulator  design  technique  for  application  to  large  space  structures 
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has  indicated  that  if  strongly  influential  or  strong] y  influenced  secondary 
modes  are  not  modeled,  spillover  can  become  a  very  serious  problem.  To  sum 
up,  some  or  all  primary  modes  and  certain  secondary  modes  are  to  be  lrodeled 
for  design  of  control  systems;  all  other  modes  are  to  be-  neglected. 

4.3  Compensation  for  Truncation 


Any  reduced-order  model  carefully  derived  a3  such  is  still  a  coarse 
approximation  of  its  original  full-order  model.  The  truncation  must  be  ad¬ 
justed  and  compensated  first  so  that  the  .closed  loop  performance  of  the  result¬ 
ing  reduced-order  model  (with  the  control  sys t  is  thus  designed)  can  closely 
represent  the  closed-loop  performance  of  the  Lall-order  model  (with  the  same 
control  systems).  In  principle,  actuators  on  the  structure  should  be  configured 
so  that  influences  on  modeled  primary  modes  are  much  stronger  than  on  any  other 
calculated  modes  and  that  the  number  of  secondary  modes  is  reduced  to  a  minimum. 
Similarly,  sensors  also  should  be  configured  so  that  influences  by  modeled  pri¬ 
mary  modes  are  much  stronger  than  by  any  other  calculated  modes  and  that  the 
number  of  secondary  modes  is  reduced  further.  The  first  step  in  compensating 
for  truncation  is  to  adjust  the  configuration  of  the  actuators  and  sensors 
properly,  since  an  initial  intuitive  configuration  generally  is  not  proper  with 
respect  to  such  specific  modeling  requirements.  It  is  ideal  if  control  spill¬ 
over  and  observation  spillover  can  be  prevented  by  proper  placement  of  a  proper 
number  of  actuators  and  sensors.  The  "placement  step"  in  the  three-step  spill¬ 
over  reduction  technique,  which  resulted  from  the  ACOSS-4  study  [4-1],  can  be 
used  to  generate  insights  of  ideal  locations  and  directions  for  actuators  and 
sensors. 

Because  of  practical  constraints  on  the  number,  type,  location,  and 
direction  of  the  actuators  and  sensors  allowed  on  the  structure,  spillover  may 
not  be  prevented  completely.  Assume  that  actuators  and  sensors  have  been 
placed  on  the  structure  and  that  the  modeled  modes  for  control  design  have  been 
determined.  Then,  the  second  step  in  the  compensation  is  synthesizing  the  in¬ 
fluences  of  those  actuators  and  sensors  properly  so  that  spillover  concerning 
a  judicious  selection  of  secondary  modes  is  prevented.  The  "synthesis  step"  of 
the  aforementioned  three-step  technique  is  applicable.  While  ideas  and  methods 
for  selecting  such  secondary  modes  are  to  be  formulated,  a  preliminary  develop¬ 
ment  of  the  synthesis  techniques  has  begun  [4-2]  and  is  being  completed  (see 
Section  4.4). 

The  third  step  in  the  compensation  is  to  filter  the  actuator  inputs  and 
the  sensor  outputs  so  that  control  spillover  and  observation  spillover  of  other 
unmodeled  modes  are  reduced  appropriately.  Low-pass  or  band-stop  filters  are 
to  be  provided  for  proper  attentuation  of  unwanted,  spill-causing  frequency 
components  in  the  inputs  and  outputs,  as  was  considered  as  the  "filtering  step" 
of  the  three-step  spillover  reduction  technique.  It  is  logical  that  control 
systems  should  be  designed  with  these  filters  considered  as  integral  parts  of 
the  reduced-order  model,  as  was  suggested  earlier  [4-1].  On  the  other  hand,  it 
also  sounds  reasonable  intuitively  that  posterior  attachment  of  low-pass  filters 
to  a  control  system  which  is  designed  using  only  the  unfiltered  reduced-order 
model  could  improve  the  full-order  closed-loop  performance  of  the  control  system. 
This  common  intuition  is  being  examined  (see  Section  4.5). 
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where  the  synthesizer  T  is  m  *  m* ,  Q2  is  m  x  (m  -  p)  and  T  is  (m  -  p)  *  m'. 

Here  m  is  the  number  of  (physical)  actuators,  m*  is  the  number  of  columns 
chosen  for  T  and  represents  the  number  of  synthetic  actuators,  and  p  is  the 

T 

rank  of  the  control  influence  matrix  4>0B_  of  such  secondary  modes. 

The  elements  of  matrix  T  are  free  parameters.  The  number  m'  of  synthe¬ 
tic  actuators  can  be  arbitrary,  but,  without  being  redundant,  m  -  p  can  be  con¬ 
sidered  an  upper  limit.  To  exercise  the  complete  synthesis  process,  m’  is 

A  A 

chosen  to  be  unity.  So  T  is  a  vector  of  free  parameters.  Consider  choosing  T 
so  that  the  synthesized  primary  control  influences  take  on  some  desirable  values 
(see  Reference  4-2,  Eq.  (5-16)).  Of  course.  Reference  4-2,  Eq.  (5-8)  will  al- 

A  A 

ways  hold  no  matter  what  the  choice  for  T  is.  Here  T  will  be  cho3en  to  make  the 
elements  of  the  synthesized  primary  control  influence  vector  all  l*s,  if  possi¬ 
ble.  The  relevant  equation  (Reference  4-2,  Eq.  (5-16))  is 

$pBpr  -  v2?  (4-2) 

T 

where  the  primary  control  influence  matrix  B_  is  p  *  m,  p  being  the  number  of 

r  r 

primary  modes,  and  V.  is  p  x  (m  -  p).  Note,  however,  that  p  is  usually  very  much 

A 

greater  than  m  -  p,  so  in  general  one  cannot  solve  Eq.  (4-2)  for  T  such  that 
Vp  U..-.S  on  the  desired  values  exactly. 

A 

Consider  using  a  linear  least  squares  technique  for  finding  F  so  that  tl\e 
error  I  U<!>p'Bpr)cjegire(j  ~  V2r||  An  International  Mathematical  and 

Statistical  Libraries  (1MSL)  linear  least-squares  subroutine  is  used  to  determine 

A 

a  minimizing  T.  A  synthesizer  T  is  then  obtained  using  Eq.  (4-1). 

The  complete  synthesis  approach  to  spillover  prevention  has  been  applied 
to  Models  1  and  2.  The  results  of  post-multiplying  the  entire  control  influ¬ 
ence  matrix  (secondary,  primary,  as  well  as  tertiary  parts)  of  both  models  are 
presented  in  the  following. 
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Figure  4-1  illustrates  the  results  for  the  example  considered  in  Sec¬ 
tion  5.5.1  of  Reference  4-2.  Recall  that  the  desired  values  for  the  elements 
of  the  synthesized  primary  control  influence  vector  are  all  l's.  Here  Modes  1 
and  2  yielded  poor  synthesized  control  influence,  while  Modes  4  and  5  synthe¬ 
sized  quite  well.  Interestingly,  Modes  1  and  2  also  liave  been  shown  to  have  a 
poor  "modal  degree  of  controllability".  The  synthesized  tertiary  influences 
were  clearly  not  of  the  order  of  the  secondary  control  influences  (as  was  de¬ 
sired),  but  ali  were  less  than  1.4.  These  are  indications  of  either  poor 
placement  of  actuators  or  poor  selection  of  secondary  modes  for  spillover  pre¬ 
vention. 
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Figure  4-1.  Synthesized  control  influence 
matrix  for  Model  1. 


This  same  synthesizer,  T,  was  used  to  post  multiply  the  control  influ¬ 
ence  matrix  of  the  perturbed  Model  i.  The  result  i3  shown  by  Figure  4-2,  where 

$  denotes  the  mode-shape  matrix  of  the  perturbed  model.  Some  reduction  was 
achieved  in  the  secondary  influences  though  the  reduction  was  not  uniform  or 
adequate.  The  problem  with  poor  placement  or  poor  selection  was  compounded  by 
parameter  variation. 
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Figure  4-2.  Synthesized  cortrol  influence 
matri:  for  perturbed  Mode*  1. 


This  same  synthesis  process  was  also  applied  to  Model  2  where  the  78  x  19 
elastic  control  influence  matrix  was  to  be  synthesized.  The  synthesizer  was 
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computed  from  the  results  of  the  second . synthesizer  example  in  Reference  4-2, 
Section  5.5.2.  The  synthesized  control  influence  matrix  for  Kodel  2  is  shown 
in  Figure  4-3.  Notice  that  again  some  of  the  primary  modes  yielded  synthesized 
values  close  to  1  while  others  did  not.  A  number  of  the  tertiary  modes  synthe¬ 
sized  to  small  values  as  was  desired.  Note  that  the  secondary  modes  chosen  for 
spillover  prevention  were  originally  selected  based  only  on  knowledge  of  the 
44  lowest  frequency  elastic  modes  and  that  the  selection  was  based  on  a  high 
"model  degree  of  controllability".  In  the  synthesized  secondary  partition,  the 
third  mode  (mode  nine)  was  several  orders  of  magnitude  above  the  other  secondary 
modes,  this  is  because  mode  nine  corresponds  to  the  row  where  the  potential 
pivots  were  too  small  to  eliminate  in  the  secondary  partition  of  the  control  in¬ 
fluence  matrix  (see  Reference  4-2,  Figure  5-3). 
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Figure  4-3.  Synthesized  control  influence  matrix  for 
Model  2. 
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4.5  Filters . for  Spillover  Reduction 

The  reduction  of  spillover  by  filtering  the  actuator  inputs  and  sensor  ■ 
outputs  appropriately  vris  originally  discussed  in  [4-1] .  The  general  effects 
of  the  filters  on  the  transient  and  steady-state  responses  have  been  investi¬ 
gated  and  are  reported  briefly  here.  Some  preliminary  results  on  posterior 
addition  of  filters  to  Model  1  are  discussed  also. 

One  may  consider  using  a  filter  to  reduce  control  spillover.  Then  one 
is  concerned  with  the  effect  of  the  filter  on  steady-state  response  and  transient 
response  with  a  general  sinusoidal  input  to  a  very  lightly  damped  system.  Two 
simple  filters  are  being  examined:  first-  and  second-order  low-pass  filters  of 
the  following  kind 


G  (s) 
c 


1 

TS  +  1 


(4-3) 


G  (s) 
c 


1 

2  2 
s  +  ds  +  c 


(4-4) 


The  time  response  (magnitude  and  phase  angle)  of  a  mode  tc  an  input 
where  filtering  is  included  has  been  determined.  The  symbol  manipulation 
language,  MACSYMA,  was  useful  for  some  of  this  work.  The  input-frequency- 
dependent  effect  of  the  filters  on  the  steady-state  responses  is  as  well  known. 

On  the  other  hand,  the  filters  reduce  the  magnitude  and  phase  angle  of  the 
transient  response  of  each  mode  uniformly  for  all  input  frequencies.  The  re¬ 
duction,  of  course,  varies  with  the  natural  frequency  of  the  mode. 

The  common  intuition  of  posterior  attachment  of  low-pass  filters  to 
vibration  controllers  was  studied  as  a  means  of  reducing  control  spillover  to 
high-frequency  modes  using  Model  1.  A  modal  dashpot  design  of  velocity  output 
feedback  controllers  (based  on  critical  Modes  1,  2,  4,  and  5)  was  considered  be¬ 
cause  of  its  commonly  recognized  robustness  in  closed-loop  stability  [4-2], 

[4-3].  CSDL  found  that  adding  a  second-order  low-pass  filter,  with  its  natural 
frequency  lying  between  Modes  5  and  6,  would  destabilize  the  12-mode  closed- 
loop  system;  intuitively  such  a  filter  was  reasonable  for  reducing  control  spill¬ 
over  to  Modes  6  end  beyond.  It  was  found  also  that  adding  a  filter  with  a  much 
higher  natural  frequency  (e.g.,  between  Modes  8  and  9)  would  not  cause  closed- 
loop  stability. 

Specifically,  the  damping  ratio  to  be  achieved  for  the  primary  modes 
with  this  output  feedback  controller  is  0.1.  Two  cases  are  considered.  In 
the  first  case,  the  feedback  gain  matrix  is 

Gv  -  -(^BF)//(2Zp0p)(Cvtp)//  (4-5) 
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and  the  control  u  is 


u  *  Gvy  (4-6) 


The  six  actuators  and  the  six  sensors  are  colocated,  respectively.  In  the  second 
case,  the  six  (physical)  actuators  are  combined  into  four  synthetic  actuators  so 
that  all  spillover  is  prevented  to  modes  three  and  six.  The  synthesizer  Ti  "hich 
is  (6  x  4) ,  is  chosen  so  that  the  synthesized  primary-mode  control  influence 


is  the  identity  matrix.  The  six  (physical)  sensors  are  also  synthesized  in  the 
same  way  but  "dualized".  So,  for  the  four  pairs  of  synthetic  actuators  and 
synthetic  sensors,  the  feedback  gain  simply  becomes 

Gy  “  -2Zp«p  (4-7) 


The  intuitively  positive  effects  of  filtering  both  actuator  inputs  and  sensor 
outputs  (toward  reducing  closed- loop  spillover  and  thus  allowing  the  damping 
ratio  to  be  achieved  by  the  primary  modes)  will  be  demonstrated  in  a  later 
report. 

For  now,  the  potentially  negative  effects  of  filter  insertion  on  system 
stability  is  investigated.  (Low-pass  filters  also  could  represent  the  dynamics 
of  the  actuators.)  The  first-order  filter  above  was  considered  for  these  values 
of  t;  0.05,  0.222,  0.25,  0.2,  1.0,  and  10.  The  system  remained  stable  for  this 
range  of  t  in  both  cases,  and  the  closed-loop  poles  of  the  12-mode  model  (plus 
filters)  moved  toward  the  ju  axis  as  r  increased.  The  stability  of  the  system 
was  not  preserved  when  the.  second-order  filter  was  inserted.  Two  values  of  c 
were  chosen  (c  =  4  and  c  *  6),  and  d  was  chosen  to  make  the  filter  critically 
damped.  When  c  was  4,  the  closed- loop  system  became  unstable  in  both  cases. 

With  c  «*  6,  however,  the  system  remained  stable  in  both  cases. 

4.6  Conclusion 

Truncated  modal  models  of  a  large  space  structure  need  to  be  selected 
properly  and  compensated  explicitly  for  control  spillover  and  observation 
spillover,  so  that  the  control  systems  thereby  designed  can  perform  satis¬ 
factorily  when  implemented  on  the  structure.  Proper  selection  requires  proper 
classification  of  structural  modes  into  "primary"  and  "secondary"  modes.  Ex¬ 
plicit  compensation  for  truncation  includes 
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(1)  Proper  placement  of  actuators  and  sensors. 

(2)  Proper  synthesis  of  the  influences  of  the  actuators  and  the 
sensors  once  they  are  placed  on  the  structure. 

(3)  Proper  filtering  of  the  actuator  inputs  and  sensor  outputs. 

Direct  applications  of  state-of-the-art  design  techniques  will  then  become 
possible,  and  the  resulting  designs  more  effective  in  closed-loop  performance. 
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SECTION  5 

ENSURING  FULL-ORDER  CLOSED-LOOP  STABILITY  IN  THE 
REDUCED-ORDER  DESIGN  OF  OUTPUT  FEEDBACK  CONTROLLERS 


5.1  Introduction 


As  a  result  of  ACOSS-6  studies,  CSDL  has  established  various  conditions 
for  ensuring  full-order  closed-loop  asymptotic  stability  and  robustness  with 
reduced-order  design  of  velocity  and  displacement  output  feedback  controllers; 
for  details,  see  [5-1]  Section  3.  The  emphasis  of  the  current  studies  is  to 
apply  such  results  to  large  flexible  space  structures  and  develop  a  reduced- 
order  design  technique  that  will  ensure  full-order  closed-loop  asymptotic 
stability. 

The  study  includes  a  preliminary  development  of  computer-aided  design 
software.  Acceleration  output  feedback  control  is  included  in  this  develop¬ 
ment  also.  The  full-or<Kr  closed-loop  asymptotic  stability  conditions  are 
therefore  extended  to  cover  the  reduced-order  design  of  acceleration  output 
feedback  controllers. 

Some  important  theoretical  or  technical  issues  that  need  to  be  ad¬ 
dressed  in  depth  in  later  work  are  discussed  briefly  at  the  end  of  this 
section. 

5.2  Problem  Formulation 

Consider  the  following  finite-element  representation  of  large  space 
structures 

Mq  +  Dq  +  Kq  =  f  (5-1) 

where  vector  q  denotes  the  L  generalized  coordinates,  vector  f  the  L  genera¬ 
lized  forces;  matrices  M,  D,  and  K  denote  the  mass  (or  inertia),  the  damping, 
and  the  stiffness,  respectively.  As  usual,  M  is  real,  symmetric,  and  positive 
definite,  while  both  D  and  K  are  real,  symmetric,  and  nonnegative  definite. 

Let  there  be  m  force  (,or  torque)  actuators  for  control  of  structural  motions 

f  *  B  u  (5-2) 

F 

where  vector  u  denotes  the  m  actuator  inputs,  one  for  each  actuator,  and  L  x  m 
matrix  BF  denotes  the  actuator  influence  coefficients.  Also  let  there  be 

acceleration  sensors,  4^  velocity  sensors,  and  displacement  sensors  for 

measurement  of  structural  motions 
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where  vector  y^  denotes  the  acceleration-sensor  outputs*  vector  y^  the 
velocity-sensor  outputs,  and  y^  the  displacement-sensor  outputs;  the  * 

L  matrix  C^,  the  Jly  x  L  matrix  Cy,  and  the  JtD  *  L  matrix  CD  denote  the  influence 
coefficients  of  the  acceleration,  velocity,  and  displacement  sensors,  respec¬ 
tively. 

Now,  consider  the  following  form  of  output  feedback  control 

("a  '  -V* 

j  “v  ■  -Vv 

(  "d  *  -Vd  (5'4) 

where  the  m  *  matrix  G^,  the  m  x  fcy  matrix  Gy,  and  the  m  x  jid  matrix  G^ 
denote  the  feedback  gains. 

Partitioning  the  actuator-influence  matrix  appropriately  yields  the 
following  expression  of  the  control  force  applied  to  the  structure 

UA 

V  ‘  lBA  i  Bv  !  V  ■  Va  +  Vv  +  Vd  (5'5) 

'v 

Substituting  Eq.  (5-2)  through  (5-4)  in  Eq.  (5-1)  yields  the  following  expression 
for  the  closed-loop  system 

(M  +  BAGACA)q  +  <D  +  VVVq  +  (K  +  BDGDCD)q  =  °  (5"6) 

It  follows  that  this  output-feedback  control  results  in  an  alteration  of  the 
mass,  damping,  and  stiffness  properties  of  the  structure.  The  feedback  gains 
Ga,  Gy,  and  Gj.  can  be  designed  using  any  technique  for  satisfying  any  speci¬ 
fied  performance  requirements,  but  the  asymptotic  stability  of  the  resulting 
closed- loop  system  must  be  ensured,  particularly  for  large  precision  space 
structures, 

A  common  practical  condition  for  asymptotic  stability  is  that  the  aug¬ 
mented  mass  matrix  (M  +  B^G^C^) ,  the  augmented  damping  matrix  (D  +  ByGyCy) , 

and  the  augmented  stiffness  matrix  (K  +  era  all  positive  definite. 

Symmetry  of  both  mass  matrix  (M  +  B^G^C^)  and  the  stiffness  matrix  (K  + 
BDGDGD^  are  always  assume<l  in  the  literature.  Such  positive-definite  condi¬ 
tions  are  not  satisfied  easily  for  large  space  structures  since  matrix  products 
BaGaCa,  ByGyCy,  and  B^G^C^  practically  can  never  be  made  positive  definite. 

Moreover,  the  design  of  control  systems  for  large  space  structures  commonly  is 
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based  on,  not  the  whole  model,  Eq.  (5-1),  but  a  truncated  version.  Thus, 
full-order  closed-loop  asymptotic  stability  is  even  more  difficult  to  ensure 
without  any  special  efforts  in  the  reduced-order  design. 

Expressed  in  modal  coordinates,  the  full-order  model,  Eq.  (5-1)  and 
(5-3),  become 


n  +  An  +  £n  =  4  f 


(5-7) 


yA  -  CkH 

yv  ■  Cv«; 

yD  -  V" 


(5-8) 


where 


(5-9) 

(5-10) 

(5-11) 

(5-12) 


To  enable  the  design  of  a  control  system,  this  modal  model  is  commonly  truncated 
to  reduce  its  order.  Let  the  undamped  natural  modes  retained  be  denoted  by 
subscipt  M  Then,  the  design  of  feedback  gains  G^,  Gy,  and  is  based  only 

on  the  following  reduced-order  model 


"k  +  Vm  +  Vk  *  V 


(5-13) 


CA$MnM 

cvVm 

cdVm 


(5-14) 


The  problem  ic  to  ensure  that  the  resulting  full-order  closed-loop  sys¬ 
tem,  given  by  Eq.  (5-6)  or  equivalently  by  Eq.  (5-2),  (5-4),  (5-5),  and  (5-8) 
combined  with  Eq.  (5-7),  is  asymptotically  stable  while  the  feedback  gains  GA> 

Gy,  and  GD  are  being  designed  using  a  reduced-order  model  such  as  given  by 

Eq.  (5-13)  and  (5-14).  For  comments  on  the  underlying  difficulties  in  ensuring 
full-order  closed-loop  asymptotic  stability  with  reduced-order  design  of  velo¬ 
city  and  displacement  output  feedback  control,  see  Reference  5-1,  Section  3. 
Those  comments  also  apply  to  reduced-order  design  of  acceleration  output  feed¬ 
back  control. 
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Full-Order  Closed-Loop  Asymptotic  Stability. Conditions  ' 

The  stability  conditions  previously  established  for  reduced-order  de¬ 
sign  of  velocity  and  displacement  output  feedback  can  easily  be  extended  to 
cover  the  design  of  acceleration  output  feedback.  Let  the  modeled  modes 
contain  those  that  can  influence  the  specified  performance  (say, . line-of-sight 
accuracy)  of  the  structure  and/or  the  sensor  outputs  significantly  and  those 
that  can  be  influenced  significantly  by  disturbance  on  the  structure  and/or 
the  actuator  inputs.  Assume  without  loss  of  generality  that  all  rigid  modes 
are  modeled  modes  and  that  all  unmodeled  modes  are  elastic  modes,  each  of 
which  has  some  positive  amount  of  inherent  damping.  Following  the  same 
development  as  (5-1]  Section  3.3,  the  following  conditions  can  be  stated. 

The  full-order  closed-loop  system  given  by  Eq.  (5-6)  is  asymptotically 
stable  if  the  following  three  conditions  are  all  satisfied. 

(1)  The  acceleration  output  feedback  gain  GA  ensures  that  the  product 
BaGaCa  is  both  symmetric  and  nonnegative  definite. 

T 

(2)  The  velocity  output  feedback  gain  ensures  that  +  ^B^G^Cy*^ 
is  positive  definite  and  ByGyCy  is  nonnegative  definite. 

(3)  The  displacement  output  feedback  gain  G  ensures  that 

T  U 

+  *mBDGDCD*M  is  Positive  definite  and  B^G^C^  is  both  sym¬ 
metric  and  nonnegative  definite. 

A  few  remarks  are  necessary.  Condition  (1)  ensures  both  the  symmetry 
and  the  positive  definiteness  of  the  augmented  mass  matrix  (M  +  BAGACA)» 

Condition  (2)  ensures  the  positive  definiteness  of  the  augmented  damping  matrix 
(D  +  ByGyCy) .  Condition  (3)  ensures  both  the  symmetry  and  the  positive  de¬ 
finiteness  of  the  augmented  stiffness  matrix  (K  +  BpGpCjj). 

5.4  Principle  of  the  Reduced-Order  Design  Technique 

Condition  (3)  takes  on  the  most  general  form  and  is  the  strongest  among 
the  three  conditions.  The  discussions  on  the  design  technique  will  focus  on 
ensuring  Condition  (3)  in  the  reduced-order  design  of  displacement  output  feed¬ 
back  control.  Thus,  the  following  reduced-order  design  technique  also  can  be 
used  or  adapted  for  reduced-order  design  of  velocity  and  acceleration  output 
feedback  control. 

Step  1?  Let  E^  be?,  a  nonnegative  definite  matrix  such  that 

SM  &  SM  +  E*  (5-15) 

is  a  positive  definite  matrix  representing  desirable  closed-loop  stiff¬ 
ness  of  the  modeled  modes,  and  set 

VdWm  “  ZM  (5-16) 
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This  ensures  that  +  ^BDGDCD$M  is  positive  definite  and  BDGDGD 
is  nonnegative  definite. 

Step  2:  Let  the  feedback  gain  0D  be  determined  as  a  product  of  .three 
matrices: 


GD  “  LUSDRD 


(5-17) 


To  ensure  the  symmetry  of  bdgqgd»  set 


Vr>  -  <W‘ 


(5-18) 


and  restrict  Ew  to  be  a  symmetric  matrix, 
M 


Step  3:  Solve  for  L^,  SD>  and  Rjj  from  the  following  matrix  equations 


k  i  7 


(5-19) 


te(«  -  £m 


(5-20) 


5.5  Discussions 


Matrix  Eq.  (5-19)  can  be  rewritten  3S 

AX  =  0 


(5-21) 


with  A  representing  the  known  part 


K !  c3 


and  X  the  unknown  part 


It  is  now  of  the  same  form  as  the  one  in  the  synthesis  of  actuator  influences 
for  prevention  of  control  spillover.  The  analytical  expression  of  the  gen-  • 
eral  closed-form  solution  already  has  been  given  in  Eq.  (2-30)  and  (2-31)  of 
Reference  5-2,  Section  2,  and  that  expression  can  be  used  directly.  A  pre¬ 
liminary  computer  program,  similar  to  that  used  in  the  synthesis  of  actuator 
influences,  has  been  coded  in  PL/I;  a  Gaussian  elimination  method  is  used. 
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Similarly,  Eq.  (5-20)  can  be  rewritten  as. 

AXB  =  C  (5-22) 

x  * 

with  A,  B,  and  C  representing  the  known  parts  (^B^Lp) ,  (R^C^^) ,  and 

respectively  while  X  represents  the  unknown  part  Sp.  Note  that  LD  and 

become  known  if  Eq.  (5-19)  is  solved  first.  The  analytical  expression  of 
the  general  closed-form  solution  recently  reported  in  Reference  5-1,  Sec¬ 
tion  3.7.5  can  also  be  used  directly.  A  preliminary  computer  program  using 
a  Gaussian  elimination  method  has  also  been  coded. 


A  preliminary  combination  of  these  two  programs  to  execute  Step  3  of 
the  design  technique  has  been  completed  and  debugged  with  Mcdc.l  No.  1.  To 
further  develop  this  design  technique,  several  theoretical  or  technical  issues 
need  to  be  addressed  in  depth.  The  general  solution  for  Eq.  (5-19)  contains 
a  matrix  of  free  parameters.  A  first  issue  to  be  addressed  is  how  to  handle 
such  free  parameters:  Is  it  useful  at  all  to  retain  the  free  parameters  in 
matrices  Lp  and  R_?  What  are  the  tradeoffs  if  the  free  parameters  are  passed 

over  to  matrix  S^? 

Both  Eq.  (5-19)  and  (5-20)  require  performing  column  and  row  operations 
on  the  known  matrices.  A  second  issue  is  the  impact  of  different  operation 
selections  on  the  feedback  gains. 

At  the  end  of  Step  3,  Condition (3)  is  ensured,  but  the  solution  Sp, 

and  hence  the  feedback  gain  G^,  is  in  general  a  function  of  free  parameters  in 

the  general  solution  of  Eq.  (5-20).  The  use  of  a  general  solution  (instead  of 
a  pseudoinverse  solution)  is  advantageous  because  it  preserves  the  available 
free  parameters.  A  third  issue  is  how  to  use  these  free  parameters  to  satisfy 
specified  performance  after  the  appropriate  stability  conditions  have  been 
ensured . 

The  design  technique  presented  in  this  section  does  not  require  collo¬ 
cation  of  actuators  with  sensors  or  the  symmetry  of  feedback  gain  matrix  G^, 

unlike  common  practice.  Positive  definiteness  of  gain  G  is  not  required 
either. 

When  this  reduced-order  design  technique  and  the  appropriate  computer- 
aided  design  software  are  developed  satisfactorily  with  respect  to  Model  No.  1, 
CSDL  will  demonstrate  their  application  to  the  much  more  complex  Model  No.  2. 
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SECTION  6 

DESIGN  FREEDOM  AND  THE  IMPLEMENTATION 
OF  SUBOPTIMAL  OUTPUT  FEEDBACK  CONTROL 


6 . 1  Introduction 

The  principal  focus  in  this  section  is  on  the  notion  of  recognizing, 
characterizing,  and  exploiting  to  advantage  the  freedom  of  choice  inherently 
available  in  design.  Although  the  specifics  discussed  relate  to  the  design  of 
feedback  controllers  for  linear  multivariable  systems,  such  a  notion  clearly 
has  much  wider  applicabiity.  It  is  unfortunate  that  potential  freedom  in  the 
design  process  sometimes  remains  unrecognized,  much  like  treasure  hidden  in  a 
field  waiting  to  be  discovered.  More  often,  perhaps,  such  freedom  is  pur¬ 
posely  sacrificed  by  making  simplifying  assumptions  to  avoid  aeriiug  with 
certain  theoretical  or  computational  difficulties  that  appear  unpleasant  and 
seemingly  insurmountable.  However,  when  sufficient  motivation  exists,  the 
effort  expended  in  facing  the  difficulties  squarely  often  provides  insight 
that  leads  to  a  much  deeper  understanding  of  the  design  problem  and  which 
reveals  a  wider  class  of  possible  solutions  [6-1],  Such  motivation  certain ’v 
exists  in  designing  active  feedback  controllers  for  vibration  suppression  ui 
large  flexible  space  structures  when  the  precision  required  in  pointinf  and 
tracking  with  the  structure  exceeds  the  present  capability.  Performane 
improvements  associated  with  tapping  previously  unexploited  design  freedom 
could  be  significant. 

The  potential  for  design  freedom  often  arises  in  connection  with  singu¬ 
larities  that  are  allowed  to  remain  in  the  mathematical  model  of  the  system 
being  studied.  Hence,  it  is  difficult  to  make  ideas  on  this  matter  precise 
apart  from  specific  applications.  In  Section  6.2,  a  number  of  specific  recent 
attempts  to  discover  and  exploit  the  freedom  of  choice  in  controller  design 
are  discussed  to  bring  out  some  important  general  ideas.  Special  attention  is 
then  given,  in  Section  6.3,  to  a  particular  design  approach  (suboptimal  output 
feedback  control  [6-2])  that  has  relevance  in  large  flexible  structure  appli¬ 
cations  which  has  been  enhanced  significantly  by  exposing  previously  hidden 
potential  for  design  freedom.  The  current  status  of  efforts  to  exploit  this 
freedom  to  advantage  is  described. 

6 . 2  Design  Freedom 

For  clarity  of  discussion,  some  notation  is  developed  briefly.  Consider 
the  standard  modal  representation  for  a  flexible  structure 


n  +  2Znn  +  fi2.;  *  (4*TB_)u 

r 

y  =»  (Cp$)n  +  (Cy$)n 
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where 


n  5  (n. . n  )  is  the  vector  of  modal  coordinates 

1  n 

8;  n  x  n  is  the  diagonal  matrix  containing  the  modal 

frequencies  to,  , . . .  ,co 
I  n 

Z:  n  x  n  is  the  diagonal  matrix  containing  the  modal 

damping  ratios  associated  with 

1  n 

inherent  material  damping 

$ :  n  x  n  is  the  matrix  with  columns  that  are  the  struc¬ 
tural  mode  shapes  , . . . , <J>n 

B  :  n  x  m  is  the  static  gain  matrix  of  the  actuators 

*  .  T 

driven  by  the  control  u  =  (uiv...,u  ) 

1  m 

Cp:  l  x  n  and  C^:  l  x  n  are  the  static  gain  matrices  of  the  displace¬ 
ment  and  rate  sensors,  respectively 


y  =  (y1*...,yi)  is  the  system  output 

Superscript  "T"  denotes  matrix  transpose 


Several  recent  investigations  into  controller  design  for  large  flexible 
space  structures  have  uncovered  substantial  potential  for  improved  system  per¬ 
formance  through  design  freedom  which  was  made  available  by  relaxing  restric- 

tions  on  the  rank  of  certain  parameter  matrices  formed  from  $  B  and  (C  ,C  )♦ 

F  P  V 

through  mode  selection.  These  studies  are  described  as  follows: 

(1)  Study  of  a  suboptimal  output  feedback  approach  [6-2]  revealed 
that  the  equation  for  the  feedback  gain  matrix  has  an  infinite 
number  of  exact  solutions  when  a  reduced  order  observation  matrix 
associated  with  (Gp,C^)$  has  less  than  full  rank  [6-3],  Specific 

indications  of  potential  benefits  to  be  gained  by  exploiting  this 
nonuniqueness  were  given. 


A  study  of  techniques  for  alleviating  control  and  observation 
spillover  associated  with  reduced-order  controllers  [6-4]  through 
the  placement  of  actuators  or  the  synthesis  of  their  influences 
revealed  that  spillover  to  certain  classes  of  modes  can  be  pre- 

T 

vented  if  a  reduced-order  matrix  associated  with  4>  B_  is  allowed 

F 

to  have  less  than  full  column  rank. 
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(3)  In  a  study  of  damping  and  stiffness  augmentation  through  output 
feedback  [6-5],  it  was  observed  that  the  equation  for  the  feed¬ 
back  gain  matrix  that  realizes  prespecified  closed-loop  damping 
(A)  and  stiffness  (K)  matrices  has  an  infinite  number  of  solu- 


(Cp,Cy)$  are  allowed  to  have  less  than  full  rank.  Specific  indi¬ 
cations  of  potential  benefits  from  such  freedom  of  choice  are 
given.  The  relaxation  of  the  usual  assumptions  on  the  closed- 
loop  matrices  A  and  K  also  contributes  to  the  freedom  of  choice. 

Similar  investigations  of  freedom  in  estimator  design  have  been  made  [6-6]. 

Freedom  of  choice  in  design  can  be  of  the  most  value  when  particular  as¬ 
pects  of  it  can  be-  characterized  and  systematic  procedures  for  effectively 
exploiting  it  can  be  developed.  This  point  is  worth  clarifying  by  describing 
a  specific  and  significant  result  of  this  nature  in  some  detail.  Consider  the 
capabilities  of  full-state  feedback  in  a  linear  time-invariant  multivariable 
system 


x  =  Ax  *  Bu  ;  u  *  Kx  +  Ln  (6-1) 

y  *  Cx  (6-2) 

where 

x  =  (x  , . .  . ,x  )T 

1  n 

T 

u  =  <V 

y  H  (yi . VT 

(n  represents  an  external  input).  The  property  of  controllability  [6-7]  has 
been  characterized  in  the  frequency  domain  as  the  ability  to  find  a  matrix  K 
such  that  the  poles  of  the  closed-loop  system  matrix  A  +  BK  formed  from 
Eq.  (6-1)  coincide  with  an  arbitrary  preassigned  symmetric  set  of  complex 
numbers  [6-8].  It  is  easy  to  show,  using  the  companion  canonical  form  [6-9] 
for  controllable  systems,  that  if  there  is  only  one  input  (m  *  1),  specifica¬ 
tion  of  all  the  (n)  closed-loop  eigenvalues  determines  the  (n)  elements  of  an 
associated  matrix  K  uniquely.  This  uniqueness  disappears  when  more  than  one 
input  is  present,  which  leaves  open  the  possibility  of  choosing  the  (nm) 
elements  of  the  feedback  gain  matrix  to  achieve  desirable  properties  of  the 
closed-loop  system  beyond  the  assignment  of  eigenvalues.  The  problem  of 
characterizing  this  freedom  of  choice  was  recognized  as  a  significant  open 
question  quite  recently  [6-10]. 


17 


A  very  elegant  characterization  of  this  freedom  was  provided  subsequently 
by  Moore  [6-11]  and  Klein  16-12]  (an  amplified  version  of  [6-11]  can  be  found 
in  16—13] ) .  In  essence,  this  characterization  states  that  the  freedom  beyond 
the  ability  to  assign  eigenvalues  is  the  ability  to  assign  associated  eigen¬ 
vectors.  Although  such  eigenvectors  cannot  be  selected  arbitrarily,  the  re¬ 
strictions  upon  the  class  from  which  the  eigenvectors  must  be  taken  are  quite 
mild.  A  precise  statement  of  this  result  (Theorem  6-2)  is  given  in  Section 


6.6.1. 


Once  such  a  characterization  is  established,  attention  is  focused  on  dis¬ 
covering  the  most  effective  means  for  exploiting  the  indicated  freedom.  Whereas 
the  assignment  of  eigenvalues  determines  the  speed  (i.e.,  natural  frequencies, 
damping  ratios)  of  the  closed-loop  dynamics,  the  assignment  of  eigenvectors  de¬ 
termines  the  shape  of  those  dynamics.  Expressing  the  general  solution  for  the 
closed-loop  dynamics  in  spectral  form  which  explicitly  displays  the  eigenvalues 
and  eigenvectors  provides  some  insight  into  possible  approaches  for  selecting 
candidate  eigenvectors  [6-13],  However,  considerable  ingenuity  is  needed  to 
find  the  most  advantageous  selections;  a  number  of  different  approaches  have 
been  proposed  [6-11,  6-12,  6-13-,  6-14,  6-15,  6-16,  6-17],  and  several  of  them 
proceed  so  as  to  reduce  the  sensitivity  of  the  eigenvalues  and  eigenvectors  to 
perturbations  in  the  closed-loop  system  matrix.  Once  a  set  of  eigenvectors 
satisfying  the  requirements  of  Theorem  6-2  has  been  selected,  the  corresponding 
gain  matrix  can  be  determined  readily  by  following  the  constructive  procedure 
given  in  the  proof  of  the  theorem  [6-13].  In  Section  6.6.2,  a  few  detailed  re¬ 
marks  regarding  the  art  of  selecting  appropriate  eigenvectors  are  given. 


6.3  Case  Study:  Numerical  Implementation  of  Suboptimal  Output  Feedback  Con¬ 
trol 


6.3.1  Background 


When  the  Kosut  approach  to  suboptimal  output  feedback  control  [6-2]  was 
reexamined  recently  to  investigate  its  applicability  to  reduced-order  control¬ 
ler  design  for  flexible  spacecraft  [6-18],  the  potential  for  substantial  free¬ 
dom  of  choice  in  gain  matrix  selection  came  to  light.  Theoretical  development 
characterizing  this  freedom  and  examples  illustrating  the  potential  for  ex¬ 
ploiting  this  freedom  are  given  in  [6—3] .  The  essential  observation  that 
generates  the  added  design  freedom  is  the  algebraic  consistency  of  a  linear 
equation  for  determining  the  output  feedback  gain  matrix,  even  when  the  coeffi¬ 
cient  matrix  is  rank-deficient.  In  such  cases  (i.e.,  rank  deficiency),  a  whole 
family  of  exact  solutions  for  the  feedback  gain  matrix  exists;  parameterization 
of  this  family  may  often  be  selected  to  improve  the  performance  of  the  full- 
order  system  which  is  driven  by  a  reduced-order  controller. 


However,  for  this  design  approach  to  be  feasible  with  large-order  sys¬ 
tems,  a  reliable  computational  procedure  for  dealing  with  rank-deficient  linear 
algebraic  equations  is  needed.  Solution  of  such  equations  is  a  very  delicate 
numerical  problem;  useless  results  are  virtually  certain  if  traditional  solu¬ 
tion  procedures  (e.g.,  Gaussian  elimination)  are  used  [6-19;  pp.  200-218],  In 
this  section,  a  mature  version  of  a  computational  procedure  initially  outlined 
in  [6-20]  for  obtaining  the  general  solution  of  rank-deficient  linear  algebraic 
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equations  is  presented.  This  procedure  takes  advantage  of  the  special  struc¬ 
ture  associated  with  the  underlying  controller  design  problem  and  expresses 
the  solution  in  a  form  particularly  appropriate  for  using  the  design  freedom 
inherent  in  the  extended  Kosut  approach  [6-21]. 

It  should  be  observed  that  the  difficulty  in  solving  rank-deficient 
linear  equations  is  not  algebraic;  the  theory  is  classical  [6-22;  pp.  73-109]. 
Rather,  the  difficulty  is  due  to  inevitable  rounding  errors  that  occur  in 
actual  computations  which  make  the  numerical  determination  of  matrix  rank 
extremely  difficult  [6-23].  Considerable  work  has  been  done  on  this  problem, 
but  very  little  of  it  has  been  reported  in  standard  numerical  analysis  text¬ 
books  [6-24,  6-25].  Attempts  to  extend  Gaussian  elimination  methods  to  treat 
rank-deficient  equations  have  not  proved  very  successful,  especially  when  an 
accurate  determination  of  rank  is  required  [6-26,  6-27].  Considerably  more 
success  has  accompanied  development  of  iterative  solution  techniques,  most  of 
which  are  designed  specifically  for  use  in  problems  with  sparse  coefficient 
matrices  [6-28,  6-29,  6-30,  6-31,  6-32].  The  roost  successful  development  to 
date  has  been  a  direct  (i.e.,  noniterative)  technique  based  on  an  algorithm  of 
Golub  and  Kahan  [6-33]  that  employs  a  decomposition  of  the  coefficient  matrix 
to  display  its  singular  values  [6-34] .  This  technique  has  been  published  as 
an  ALGOL  procedure  [6-35]  and  incorporated  into  widely  available  tobust  mathe¬ 
matical  software  [6-36].  It  is  widely  recognized  as  the  only  reliable  ap¬ 
proach  to  problems  in  which  an  accurate  determination  of  numerical  rank  is 
essential  [6-23,  6-34].  A  precise  development  of  what  "numerical  rank"  really 
means  has  appeared  recently  [6-37]. 

For  the  present  application,  precise  determination  of  the  rank  of  the 
coefficient  matrix  in  the  gain  equation  is  very  important  since  the  number  of 
free  parameters  in  the  general  solution  for  the  gain  matrix  is  proportional  to 
the  rank  deficiency  of  the  coefficient  matrix  [6-3].  Moreover,  the  structure 
of  the  underlying  control  design  problem  suggests  that  the  coefficient  matrix 
generally  will  be  dense  rather  than  sparse.  Therefore,  the  procedure  for 
numerical  solution  of  the  Kosut  gain  equation  is  built  around  the  singular 
value  decomposition. 

6.3.2  Problem  Statement 

The  algebraic  equation  to  be  solved  numerically  has  the  form 

AX  *  B  (6-3) 


With  reference  to  the  underlying  problem  of  reduced-order  controller  design 
[6-3],  matrices  A:  £  x  l  and  B:  i  x  m  have  the  structure 


AT  =  CcPcJ  ,  BT  =  F* pcj  (6-4) 
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where 


Cj,:  £  x  2N  is  the  state-to-output  observation  matrix  in  a  reduced-order 
plant  having  2N  states,  m  inputs,  and  £  outputs 

P:  2N  x  2N  is  the  symmetric  positive-definite  solution  of  a  Liapunov 
matrix  equation  associated  with  the  stable  closed-loop 
system  matrix  of  a  reference  plant  of  order  2N 

F  :  m  x  jj,  is  the  full-state  feedback  gain  matrix  in  that  reference 
plant 

T 

The  gam  matrix  to  be  determined  is  G  *  X  :  m  x  £. 

It  follows  directly  from  the  structure  of  A  and  B  (Eq.  (6-3)  and  (6-4)) 
that  the  gain  equation,  Eq.  (6-3),  is  always  algebraically  consistent,  regard¬ 
less  of  the  rank  of  A  (=  rank  (C^))  [6-3] . In  other  words,  the  best  least- 

squares  approximation  to  a  solution  of  Jq.  (6-3)  is,  in  fact,  an  exact  solu¬ 
tion.  When  the  coefficient  matrix  A  is  rank  deficient,  say  rank  (A)  a  r  <  £, 
the  general  solution  of  Eq.  (6-3)  contains  m  •  (£  -  r)  arbitrary  parameters 
(6-3].  The  problem  treated  is  to  develop  a  procedure  for  computing  the  gen¬ 
eral  solution  to  Eq.  (6-3).  The  procedure  is  designed  so  that  the  free  param¬ 
eters  are  displayed  explicitly  in  a  manner  that  facilitates  their  selection 
for  improving  the  performance  of  the  full-order  (finite-dimensional)  plant 
driven  by  the  reduced-order  controller  in  the  underlying  control  design  prob¬ 
lem.  Before  this  is  described,  a  few  essential  facts  about  generaliz'd  in¬ 
verses  and  the  singular  value  decomposition  are  reviewed. 

6 . 4  Some  Facts  from  inear  Algebra 

6.4.1  Generalized  lnvi  rses 

A  matrix  is  invertible  if  it  is  square  and  nonsingul^r .  Otherwise,  in 
general,  some  restriction  of  its  domain  (viewing  the  matrix  as  a  mapping)  is 
required  before  any  property  of  invertibility  can  be  achieved.  The  notion  of 
a  generalized  inverse  makes  the  appropriate  restrictions  and  the  corresponding 
invertibility  properties  precise. 

The  Moore-Penrose  inverse  of  a  real  matrix  A:  p  x  v  is  the  unique  solu¬ 
tion  of  the  equations  16-38,  6-39] 


AXA  " 

A 

(6-5) 

XAX  » 

X 

(6-6) 

(ax)t  “ 

AX 

(6-7) 
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(6-8) 


and  is  denoted  by  A  :  v  *  u,  It  is  identical  with  the  usual  inverse  of  A  when 
A  is  square  and  nonsingular .  A  key  to  understanding  this  notion  is  the  geo¬ 
metric  interpretation  of  the  composite  maps  A  A  and  AA  as  projections.  In 
particular  [6-38] 

A^A  is  a  projection  onto  R(A^)  along  N(A),  (6-9) 

AA^  is  a  projection  onto  R(A)  along  N(A^) 


where  R(»)  and  N(»)  denote  the  range  and  null  spaces,  respectively,  of  their 
arguments.  That  those  projections  are  orthogonal  follows  from  the  relations 
[6-40] 


R(Af)  =  N(A)^,  N(A!)  *  R(A)^ 


which  can  be  deduced  directly  from  Eq.  (6-5)  through  (6-8);  the  superscript!, 
denotes  the  subspace  orthogonal-complement  relative  to  the  usual  inner  product 
in  finite-dimensional  spaces.  Complementary  (orthogonal)  projections  are 

Iy  -  A  A  and  I  -  Aa\  respectively,  where  I  denotes  the  v  x  v  identity  ma¬ 
trix.  The  appropriate  restriction  of  the  domain  of  A  is  now  evident;  restrict¬ 
ing  A  to  N(A)-**  converts  A^A  to  the  identity  map  on  N(A)"^. 

The  Moore-Penrose  inverse  can  always  be  represented  explicitly.  The  de¬ 
composition  of  a  general  matrix  in  the  form 


where  B:  y  x  p ,  C:  p  x  v,  and  p  =  rank  (A),  is  always  possible  [6-41],  Repre¬ 
sentations  for  and  are  [6-42] 


t  ,  T.-l  T 
B  =  (B  B)  B  , 


cWr1 
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and  finally 


However,  it  should  be  noted  that  the  relation 


(BC)*  -  cV 


is  not  always  true  for  arbitrary  matrices  B  and  C  [6-43] . 

One  of  the  chief  applications  of  the  Moore-Penrose  inverse  is  in  repre¬ 
senting  the  solutions  to  linear  algebraic  equations  of  the  (vector)  form 


Ax  ■  b 


(6-10) 


If  Eq.  (6-10)  is  algebraically  consistent,  then  its  general  solution  can  be 
written  as  [6-38] 


A*b  +  (I£  -  AfA)w 


(6-11) 


where  oi  is  an  arbitrary  vector.  Note  that,  since  b  is  in  R(A),  Eq.  (6-11) 
is  an  orthogonal  decomposition  (cf.  Eq.  (6-9)).  If  Eq.  (6-10)  is  not 

algebraically  consistent,  then  the  product  A^b  is  the  least-squares  approxi¬ 
mation  of  minimum  norm.  Further  details  on  the  Moore-Penrose  and  other  gener¬ 
alized  inverses  may  be  found  in  [6-44]. 

6.4.2  Singular  Value  Decompositiou 

.  T  T 

The  symmetric  products  A  A  and  AA  of  a  rectangular  matrix  A  have  the 
same  nonzero  eigenvalues,  and  each  has  the  same  rank  as  that  of  A.'  The  non- 

T 

negative  square  roots  of  the  (real)  eigenvalues  of  A  A, are  called  singular 
values  of  A.  They  have  importance  for  numerical  linear  algebra  because  they 
are  much  less  sensitive  to  perturbation  in  the  elements  of  A  than  are  the 
eigenvalues  of  A  [6-23] . 

Any  rectangular  matrix  A:  u  x  v  with  rank  p  can  be  decomposed  to  display 

its  singular  values  a  >  ...  >  a  >  0  »  o  ,  ■  ...  ■  o  as  follows  [6-34] 

l  —  —  p  P+1  n 


A  -  UEV 


(6-12) 


where  U:  y  x  y  and  V:  v  x  v  are  orthogonal  (unitary  if  A  is  complex)  matrices, 
and  I:  y  x  v  is  a  partitioned  matrix  with  D  ■  diag  (o  ,...,c  )  in  the  upper 

1  P 

left  p  x  p  block;  the  remaining  blocks  are  zero.  The  simplicity  of  existence 
proofs  for  this  decomposition  [6-34,  6-45]  stands  in  uharp  contrast  to  the 
involved  details  of  reliable  algorithms  for  coiaput ing  it  [6-33],  This  is  a 
reflection  of  the  comparative  difficulty  between  analytical  and  computational 
determination  of  matrix  rank  since  the  number  of  nonzero  singular  values  is 
equal  to  the  rank  of  the  matrix. 

The  columns  of  U  and  V  in  the  singular  value  decomposition  Eq.  (6-12) 

have  an  important  geometrical  interpretation.  If  U  =  [Upl^]  and  V  =  [Vj'.Vj] 

are  partitioned  to  be  compatible  with  E,  it  follows  from  writing  Eq.  (6-12)  in 
the  form  AV  *  UE  that  [6-33] 


Columns  of  U^U^)  form  a  basis  for  R(A)  (R(A)'^") 


Columns  of  Vj(V^) 


form  a  basis  for  N(A)”^  (N(A)) 


(6-13) 


These  interpretations  are  used  to  develop  an  expression  for  the  general  solu¬ 
tion  of  the  Kosut  gain  equation,  Eq.  (6-3),  appropriate  for  the  intended  ap¬ 
plication. 


Finally,  it  should  be  observed  that  the  Moore-Penrose  inverse  is  easily 
expressed  in  terms  of  this  decomposition 


A 


t 


t  T 
VE  V 


(6-14) 


where  the  only  nonzero  block  of  E^  is  the  upper  left  p  x  p  block  D  ^  *  diag 

,  -1  -1, 

,  •  •  •  >Op  / . 

6.4.3  An  Expression  for  the  General  Solution 

In  examining  solutions  of  the  matrix  Eq.  (6-3),  it  is  sufficient  to 
consider  vector  equations  of  the  form  of  Eq.  (6-10),  It  is  assumed  that  Eq. 
(6-10)  is  algebraically  consistent.  It  was  shown  previously  [6-3]  that  the 
general  solution  to  Eq.  (6*10)  can  be  expressed  in  the  form 


x 


Afb  + 


SY 


(6-15) 
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where  the  columns  of  S:  1  x  (1  -  r)  form  a  basis  -tor  HOa)  and-Ys  (*•  “  r)  x  1 
is  arbitrary.  The  classical  expression  Eq.  (6-11)  for  the  general  solution  of 
Eq.  (6-10)  differs  from  Eq.  (6-15)  in  a  qualitative  manner  that  is  important 
for  the  intended  application;  namely,  it  appears  that  there  are  more  arbitrary 
parameters  available.  This  is  misleading. 

Theorem  6-1 

Equations  (6-11)  and  (6-15)  are  equivalent  solutions  of  Bq.  iu-lO). 


Proof 

If  A  has  maximum  rank,  there  is  nothing  to  prove,  since  then  A^A  m  1^ 
and  S  is  the  empty  matrix.  Assume  rank  (A)  *  r  <  Jt, 

(+)  Let  id;  l  x  1  be  arbitrary,  and  investigate  the  solvability  of  the 
linear  system 


SY  -  (1^  -  A^A)w 


for  y:  (#-  ~  r)  *  1.  It  follows  from  Eq.  (6-9)  that  1^  -  A*A  is  a  projection 

1  l-r 

onto  N(A).  Since  the  columns  s  ,...,s  of  S  form  a  basis  for  N(A),  there 

exist  numbers  Y  =  (Y, , . . .  ,Y.  )  such  that 

1  x-r 


£-r  .  . 

ST  5  I  Y-sX  -  (I  -  ATA)ui 
i-1  1  * 


(♦)  Let  y:  (1  -  r)  x  1  be  arbitrary,  and  investigate,  the  solvability  of 
the  linear  system 

(1^  -  A^A)w  »  SY 


for  u.  By  the  definition  of  S:  Sy  e  K(A).  Since  1^  -  aV  is  a  projection 
with  range  equal  to  N(A),  there  exists  such  anw.B 

It  is  clear  from  this  result  that  there  are  at  most  i  -  r  "effectively 
independent"  choices  of  free  parameters  in  the  general  solution  of  Eq.  (6-10). 
The  representation  in  Eq,  (6-15)  is  therefore  sore  appropriate  for  the  appli¬ 
cation  intended,  and  it  is  the  foundation  for  the  solution  procedure  to  be 
developed.  More  precisely,  the  matrix  form  of  Eq.  (6-15)  corresponding  to 
Eq.  (6-3)  i3  used 
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X 


A*B  +  sr 


where  T:  (t  -  r)  *  in  is  arbitrary.  Note  that  the  matrix  of  Eq.  (6-13) 

forms  a  natural  choice  for  a  basis  of  N(A)  and  is  available  directly  from  the 
decomposition,  Eq.  (6-12). 

6.4.4  A  Numerical  Procedure 

In  the  numerical  procedure  to  be  described,  the  focus  is  on  the  case 
that  the  coefficient  matrix  A  in  Eq.  (6-3)  is  rank  deficient.  Otherwise,  the 
solution  is  unique,  there  are  no  free  parameters,  and  standard  linear- 
equation-solving  techniques  can  be  used  unless  there  is  some  reason  to  suspect 
A  to  be  ill-conditioned  with  respect  to  matrix  inversion.  The  rank-deficient 
case  is  the  most  interesting  because  of  its  beneficial  potential  for  the 
underlying  controller  design  problem. 

The  procedure  consists  of  three  basic  parts:  determining  the  numerical 
rank  of  A;  computing  the  particular  solution  A  B;  and  computing  the  general 
solution  A^B  +  sr. 

6. 4. 4.1  The  Procedure 

Step  1.  Determine  the  numerical  rank  of  the  coefficient  matrix 

Tbenote  by  rk  (•)  the  numerical  rank  of  the  indicated  matrix. ) 

}  n 

A  T 

(a)  Form  the  product  A  3  C  PC„  (cf.  Eq.  (6-4)). 

L  Ci 

(b)  Compute  the  singular  value  decomposition  of  A,  obtaining  matrices 
U  and  V,  and  singular  values  Oj,...,o^  (cf.  Eq.  (6-12)). 

(c)  Make  an  intelligent  judgment  as  to  which  singular  values  °r+l *  *  * ’ 
o ^  should  be  considered  zero  (i.e.,  those  considered  to  consist 
simply  of  an  accumulation  of  roundoff  errors)  (cf.  Note  1). 

(d)  Tentatively  assign  rkn(A)  M  r. 

(e)  Estimate  the  numerical  rank  of  the  sensor  matrix  Cc  by  repeating 
steps  corresponding  to  Items  (b)  and  (c).  If  rkn(Cc)  *  rkn(A), 

8 top  execution  and  investigate  (cf.  Note  2). 


Form  the  product  B  ■  C  PF  (cf.  Eq.  (6-4)). 

Ci 
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(g)  Estimate  the  numerical  rank  of  the  augmented  matrix  [AlBj  by  re¬ 
peating  steps  corresponding  to  Items  (b)  and  (c).  If  rkn([A!B]) 

A  rk^A),  stop  execution  and  investigate  (cf.  Note  3). 

(h)  When  the  estimates  for  numerical  rank  in  Steps  (e),  (f),  and  (g) 
agree,  assign  to  rV  (4)  their  common  value. 

Step  2.  Determine  the  particular  solution  A*B. 

(a)  Compute  A^  using  the  matrices  U  and  V  from  Step  1(a)  and  the  diag¬ 
onal  matrix  ®  di*g  (o.  , ...,o  *),  where  rk  (A)  (cf.  Eq. 

U  l  rQ  u  n 

(6-14)). 

(b)  Compute  A^B  (cf.  Note  4). 

Step  3.  Determine  the  general  solution  A^B  +  Si 

(a)  Select  a  matrix  T:  (i.  -  r^)  x  m  of  parameter  values  (arbitrary). 

(b)  Choose  a  matrix  S:  t  x  U  -  rQ)  with  columns  that  form  a  basis 
for  N(A)  (cf.  Note  5). 

(c)  Compute  the  general  solution  X  =  A^B  +  Sr  (cf.  Note  6). 

6. 4. 4. 2  Notes  Referenced  in  the  Procedure 

(1)  This  task  is  nontrivial,  and  an  understanding  of  the  underlying 
physical  problem  is  usually  essential  [6-46] .  The  purpose  of 
Steps  1(e)  and  1(g)  is  to  contribute  some  of  this  insight  to  the 
final  assignment  of  a  value  to  rk^(A). 

(2)  Theoretically,  rank  (C^,)  3  rank  (A)  [6-3],  In  most  nses,  com¬ 

putation  of  the  numerical  rank  of  C_  will  be  subject  to  less 

error  than  computation  of  the  numerical  rank  of  A-  In  some 

cases,  it  may  be  possible  to  determine  rank  (C_)  by  inspection, 

c 

which  then  determines  the  numerical  rank  of  A. 

(3)  Theoretically,  rank  ([AtB])  =  rank  (A)  [6-3]. 

(4)  In  actual  computation,  A^  need  not  be  calculated  explicitly;  the 

solution  A  B  is  more  efficiently  computed  as  the  product  (cf.  Eq. 
(6-14)) 
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(5)  The  basis  chosen  need  not  be  a  mutually  orthogonal  one.  However, 
an  orthogonal  basis  is  readily  available  as  the  set  of  those 
columns  of  V  (cf.  Eq.  (6-12))  which  correspond  to  the  zero  singu¬ 
lar  values  of  A  (cf.  Eq.  (6-13)). 

(6)  For  a  given  control  design  problem,  Steps  1  and  2  of  the  proce¬ 
dure  need  only  be  done  once.  Step  3  may  be  computed  repeatedly 
with  different  choices  for  T  and  S  until  satisfactory  control 
system  performance  improvement  is  achieved. 

6.4.5  Summary 

In  order  to  be  able  to  take  advantage  of  the  design  freedom  inherent  in 
the  extensions  [6-3]  to  the  Kosut  method  of  suboptimal  output  feedback,  a 
numerical  procedure  has  been  developed  to  compute  the  general  solution  to  the 
gain  equation  in  the  rank-deficient  case.  The  procedure  is  based  upon  a  com¬ 
putation  of  the  singular  value  decomposition  for  the  coefficient  matrix  and 
expresses  the  solution  in  a  form  that  facilitates  exploiting  the  available  de¬ 
sign  freedom.  'The  design  freedom  consists  of: 

(1)  A  choice  of  the  basis  for  the  null  space  of  the  coefficient 
matrix. 

(2)  A  choice  of  free  parameters  proportional  in  number  to  the  rank 
deficiency  of  the  coefficient  matrix. 

One  candidate  for  a  choice  of  the  basis  vectors  is  provided  explicitly  by  the 
procedure.  Choice  of  the  free  parameters  may  be  motivated  by  consideration  of 
the  full-order  closed-loop  system  equations  that  contain  the  reduced-order 
controller  under  study. 

6.5  Conclusion 


Considerable  effort  has  been  devoted  to  characterizing  the  design  free¬ 
dom  with  suboptimal  output  feedback  that  was  recently  uncovered  and  to  devel¬ 
oping  methods  for  exploiting  that  freedom.  In  the  present  section,  this 
effort  has  been  placed  in  a  much  broader  context  by  showing  the  close  rela¬ 
tionship  with  other  significant  work  relating  to  design  freedom  and  indicating 
the  rapid  development  of  interest  in  such  matters.  We  believe  that,  in  re¬ 
sponse  to  increasingly  stringent  spacecraft  performance  specifications,  an  in¬ 
creasing  amount  of  attention  will  be  given  to  discovering,  characterizing,  and 
developing  systematic  approaches  to  exploiting  available  freedom  of  choice  in 
design. 
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6.6  Addenda 


6.6.1  Characterization  of  Design  Freedom  with  State  Feedback  (m  >  1) 

The  system  of  Eq.  (6-1)  is  being  considered.  In  order  to  facilitate  a 
concise  statement  of  the  main  result,  several  preliminary  ideas  are  neces¬ 
sary.  First,  recall  an  alternate  characterization  of  controllability  in  the 
frequency  domain: 

Fact  6-1  [6-47;  pp.  71,  161].  The  pair  (A,B)  is  controllable  if  and  only  if 
the  polynomial  matrix  X  S(X)  *  [XI^  “  AlB]  ^as  linearly  independent  rows  for 
each  complex  number  X.fl 


Second,  note  a  particular  property  of  such  polynomial  matrices: 

Fact  6-2  [6-47;  p.  194].  If  (A,B)  is  controllable,  there  exists  a  matrix 
function  X  which,  for  each  complex  number  X,  forms  a  basis  for  the 

null  space  of  S(X).I 

The  main  result  is  stated  in  terms  of  the  n  x  m  matrix  function  N(»): 
Theorem  6-2  [6-13] 

Assume  that  the  pair  (A,B)  is  controllable  and  that  a  set  of  distinct 
complex  numbers  X^,  i  *  l,...,n  (containing  the  conjugates  of  each  of  its 

elements)  is  given.  Then  there  exists  a  real  full-state  feedback  matrix  K 
(cf.  Eq.  (6-1))  such  that 


(A  +  BK)vl 


X.v  ,  i  *  1 , . . .  ,n 


if  and  only  if  the  (eigenvectors)  v1  satisfy  the  following: 

(1)  The  collection  {v1}/^  is  linearly  independent  over  the  complex 
field. 


(2) 


The  collection  of  related  pairs  { (X^,v1)}  are  such  that 
v^  **  vl  if  Xj  *  X^  (where  " — "  denotes  complex  conjugation) 


(3) 


For  each  i:  v  c  span  [N(X^)]. 


Mcveoever,  any  \  K  is  unique. 
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This  result  can  be  modifiH  appropriately  to  eliminate  the  need  for 
assuming  controllability  of  the  pair  (A,B)  [6-11]  or  for  assuming  that  the 
assigned  closed-loop  eigenvalues  are  distinct  [6-12], 

6.6.2  Exploitation  of  Design  Freedom  with  State  Feedback  (m  >  1 ) 

A  spectral  characterization  of  the  closed-loop  output  dynamics  of 
Eq.  (6-1)  and  (6-2)  is  given  readily.  For  simplicity,  assume  distinct  and 
stable  closed-loop  eigenvalues  and  a  step  input. 

Fact  6-3  (6-13] .  Assume  the  eigenvalues  of  A  +  BK  are  distinct  with  negative 
real  parts.  Let  n(t)  =  ng  =  constant.  Then 

y(t )  -  y  =*  CVeAtv"1x(0)  +  CVeAtA'1v"1BLn  ,  t  >  0 

“  s  ~ 

and 

y  =  -  CVA”1V~IBLn 

00  s 

where  A  =  diag  (A J  displays  the  eigenvalues,  and  V  =  [v*I,...,I  vN]  displays 
the  eigenvectors,  respectively,  of  A  +  BK.flj 

The  reader  is  referred  to  (6-13]  for  specific  suggestions  on  how  to  use 
this  representation  effectively  as  a  guide  in  selecting  the  eigenvectors  to 
shape  the  dynamics  of  the  response  y(t)  -  y  .  In  particular,  the  potential 

for  eliminating  effects  of  a  given  mode  upon  as  many  as  m  -  1  of  the  output 
components  is  of  interest.  It  should  be  noted  that  much  more,  freedom  of 
choice  really  exists  than  is  actually  exploited  in  [6-13],  since  all  but  one 
of  the  eigenvectors  are  assumed  fixed, 
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SECTION  7 


STOCHASTIC  OUTPUT  FEEDBACK  COMPENSATORS 
FOR  DISTRIBUTED  PARAMETER  STRUCTURAL  MODELS 


7.1  Introduction 

This  section  presents  recent  progress  on  the  stochastic  output  feedback 
design  problem  for  distributed  parameter  plants.  The  results  presented  are  an 
extension  of  those  described  in  Section  2  of  Reference  7-1.  The  desire  to  de¬ 
sign  optimal  output  feedback  compensators  for  distributed  parameter  flexible 
models  motivated  consideration  of  the  following  problem. 

Given  a  cost  functional  J(*)  defined  over  some  set  T  of  admissible  com¬ 
pensator  designs  for  the  unknown  distributed  parameter  plant  P,  determine  the 

'ft 

optimal  design  K  that  minimizes  J(»)  over  T. 

To  accomplish  this,  define  a  sequence  of  finite-diraensionai  approximat¬ 
ing  plant  models  {P^}  along  with  an  associated  sequence  of  cost  functionals 

{ Jn( - )}  and  determine  the  optimal  compensator  designs  that  minimize  J^i*) 
over  F.  For  "reasonable"  choices  of  P,  F,  J(*),  and  the  approximating  scheme 
defining  Pn  and  J^C • ) ,  one  would  expect  that  the  optimal  compensators 

* 

should  converge  to  K  as  n  +  «.  However,  it  is  not  obvious  exactly  what  con¬ 
stitutes  a  "reasonable  choice"  of  P,  T,  J(»),  {P^},  and  For  example 

[7-2],  in  the  special  case  when  P  is  an  undamped  flexible  structure,  linear 

regulators  designed  from  truncated  normal  mode  models  P  do  not  converge  to  a 

n 

limiting  linear  regulator  for  the  plant  P.  Consequently,  it  is  extremely  im¬ 
portant  to  address  the  following  fundamental  questions: 

(1)  When  does  a  set  S  C  T  of  optimal  designs  for  the  plant  P  exist? 

(2)  When  do  the  corresponding  sets  Sr  exist  for  the  plants  Pn? 


(3)  If  the  sets  exist  for  all  n,  does  the  sequence  {S^}  converge  to 
some  set  S  C  r? 

(4)  If  the  sequence  {S^}  converges,  does  it  converge  to  S  ? 

Sufficient  conditions  are  developed  in  this  section  for  the  answers  to 
each  of  these,  questions  to  be  at  least  approximately  "yes".  The  results 
reported  previously  [7-1]  addressed  only  the  last  of  these  questions,  provid¬ 
ing  sufficient  conditions  for  the  optimality  of  the  limiting  set  of  compensa- 

tors  to  which  the  sequence  (S^J  converged,  assuming  the  sequence  converged. 


Continued  investigation  has  shown  that  some  of  the  key  theorems  employed  in 

establishing  these  results  are  special  cases  of  more  general  results  [7-3]. 

In  particular,  the  results  of  Reference  7-3  are  sufficient  to  guarantee  the 

*  *  *  * 
existence  of  the  sets  S  and  S^,  the  convergence  of  the  sequence  {S^},  and  the 

optimality  of  the  limit.  Further,  these  results  are  applicable  to  infinite- 
dimensional  design  sets  F,  which  suggests  the  possibility  of  developing  con¬ 
vergence  conditions  for  infinite-dimensional  compensators  like  the  linear 
regulator  considered  in  Reference  7-2.  The  price  paid  for  this  increased 
generality,  however,  is  the  necessity  of  dealing  with  weak  convergence  and  re¬ 
lated  Banach  space  concepts.  If  these  results  are  specialized  to  the  original 
problem  for  which  T  is  finite  dimensional,  the  strong  convergence  results  con¬ 
jectured  in  the  conclusions  of  Reference  7-i ,  Section  2  are  obtained. 

Consequently,  the  results  of  this  section  are  presented  in  the  following 
format.  First,  the  notion  of  weak  convergence  and  weak  compactness  are 
introduced  in  Section  7.2,  and  their  importance  is  established  in  the  context 
of  the  problem  considered  herein.  Next,  the  general  results  of  Reference  7-3 
are  presented  in  Section  7.3.  Specifically,  sufficient  conditions  on  J(*), 

iff 

Jn(*),  and  T  are  presented  for  sets  S  and  Sn  of  global  minimizers  of  J(*)  and 
Jn(*)  to  exist  along  with  sufficient  conditions  to  guarantee  that  the  sequence 

•fa  ^ 

(Sn)  converges  weakly  to  a  subset  of  S  as  n  00 .  Once  these  general  results 

are  established,  they  are  specialized  in  Section  7. A  to  the  original  problem 
where  the  set  T  over  which  J(0  is  minimized  is  finite  dimensional.  This 
yields  a  variety  of  existence  and  convergence  results  for  the  optimal  fixed- 
form  compensator  design  problem,  which  culminates  in  Theorem  7-6,  Theorem  7-6 
provides  sufficient  conditions  for  the  answer  to  the  four  questions  posed  pre¬ 
viously  to  be  approximately  "yes".  The  implications  of  the  requirements  and 
conclusions  of  Theorem  7-6  are  discussed  briefly  in  Section  7-5  to  make  these 
answers  more  precise. 


Notion  of  Weak  Convergence  and  Weak  Compactness 


In  the  results  developed  in  Reference  7-1,  the  compactness  of  closed, 
bounded  subjects  of  the  set  T  of  admissible  compensators  played  an  essential 
role.  If  T  is  defined  on  an  infinite-dimensional  Banach  space,  however,  it  is 
no  longer  necessarily  true  that  closed,  bounded  subsets  of  T  are  compact.  For 

example,  the  set  {x  e  b|  ||x||  £  M}  is  compact  if  and  only  if  the  space  B  is 

finite  dimensional  [7-4].  Consequently,  to  extend  the  results  of  Reference 
7-1  to  problems  for  which  T  is  defined  on  an  infinite-dimensional  Banach 
space,  it  is  necessary  to  introduce  the  notions  of  weak  convergence  and  weak 
compactness. 


First,  recall  that  a  Banach  space  is  a  normed  linear  vector  space  that 
is  complete  in  its  norm  (e.g.,  given  a  sequence  {x^}  such  that  x^  e  B  for  all 

n  and  jjx^  -  xj J  0  as  n  -*•  “,  it  follows  that  x  e  B).  Convergence  of  this 
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<  e  for 


type  (e.g.,  given  e  >  0,  there  exists  N  such  that  n _>  N  *>  jjx^  -  x| j 

all  n  >.  N)  is  a  strong  convergence,  denoted  xr  ■?  x,  and  it  corresponds  to  the 

usual  notion  of  convergence  in  the  Euclidean  space  Rn  for  any  finite  n. 

If  B  is  an  infinite  dimensional  Banach  space,  then  a  weaker  type  of  con- 

* 

vergence  is  possible.  Defining  B  as  the  dual  space  of  B  Ci.e.,  the  space  of 

bounded  linear  functionals  mapping  B  into  R  ),  {x}  converges  weakly  to 

x  (x  %  x)  if  |f(xn)  ~  f(x)j  0  as  n  ®  for  all  f  e  B  .  The  implication  that 

convergence  implies  weak  convergence  is  a  standard  result  [7-4,  7-5],  but  the 
converse  is  not  true  unless  B  is  finite  dimensional. 


*  ** 

Since  B  is  a  Banach  space,  the  second  dual  space  B  of  B  is  well  de¬ 
fined  and  it  is  another  standard  result  [7-6]  that  B  may  be  identified  with  a 
'■  ** 

subspace  B  of  B  .  In  general,  this  inclusion  is  proper,  but  a  reflexive 

~  ** 

Banach  space  is  one  for  which  B  *  B  .  This  is  a  very  broad  class  of  spaces, 
which  includes  all  Hilbert  spaces  and,  as  a  consequence,  ail  finite- 

dimensional  Euclidean  spaces  Rn.  For  the  problem  considered  here,  the  most 
significant  feature  of  this  class  of  spaces  is  that  all  bounded  sequences  {xn] 

(i.e.,  all  sequences  such  that  ||xn||  <  M  for  all  n)  have  weakly  convergent 

subsequences  [7-5].  Consequently,  any  countable,  bounded  sequence  on  a  re¬ 
flexive  Banach  space  may  be  decomposed  into  (possibly  a  countably  infinite 
number  of)  convergent  subsequences,  each  of  which  converges  weakly  to  some 
l imi t  in  B . 


This  fact  may  be  used  to  define  the  notion  of  weak  compactness  on  re¬ 
flexive  Banach  spaces.  Specifically,  in  analogy  with  sequential  compactness, 
a  set  S  is  weakly  sequentially  compact  if  every  sequence  (xn)  in  S  has  a  weak¬ 
ly  convergent  subsequence  that  converges  to  an  element  of  S.  Since  compact¬ 
ness  and  sequential  compactness  are  equivalent  on  metric  spaces  [7-4],  the 
term  "weak  sequential  compactness"  can  be  shortened  to  "weak  compactness"  or 
"w-compactness".  Similarly,  weakly  closed  (w-closed)  sets  may  be  defined  as 
sets  containing  all  of  their  weak  limit  points  (i.e.,  S  is  w-closed  if  x^  e  S 

for  ail  n  and  xn  ^  x  together  imply  x  e  S).  It  then  follows  from  the  preced¬ 
ing  result  that  if  S  is  a  weakly  closed,  bounded  subset  of  a  reflexive  Banach 
space  B,  then  S  is  weakly  compact.  Thus,  the  concepts  of  weak  convergence  and 
weak  compactness  play  the  same  roles  in  infinite-dimensional,  reflexive  Banach 
spaces  that  their  strong  counterparts  play  in  finite-dimensional  spaces. 

Finally,  it  is  important  to  note  that  weak  closedness  is  a  "stronger" 
notion  than  ordinary  (i.e.,  strong)  closedness  in  spite  of  its  name.  Speci¬ 
fically,  if  a  set  is  weakly  closed,  it  is  also  strongly  closed  since  if  xr  e  S 

for  all  n  and  x  $  x  then  x^  ^  x,  implying  x  e  S.  The  converse,  however,  is 


not  generally  true.  A  special  case  for  which  the  two  concepts  are  equivalent 
occurs,  though,  when  S  is  convex  also  [7-3]  (i.e.,  given  e  ^  f°H°w8 

that  Xxj  +  (1  -  X)x^  e  S  for  ail  0  £  X  £  1).  Thus,  if  a  subset  S  of  a  reflex¬ 
ive  Banach  space  B  is  closed,  bounded,  and  convex,  it  is  weakly  compact. 


General  Existence  and  Convergence  Results 


The  following  paragraphs  summarize  the  pertinent  results  of  Reference 

•  •  *  .  . 

7-3  concerning  the  existence  of  global  mintmizers  x  of  a  functional  f C  * )  de¬ 
fined  on  a  reflexive  Banach(  space  B  and  the  weak  convergence  of  the  global 
*  * 
miniraizers  x^  of  a  sequence  of  approximating  functionals  { f n ( • ) }  to  x  as 

n  °°.  In  particular,  very  general  sufficient  conditions  are  developed  for 
both  of  these  situations  to  occur.  More  specifically,  the  first  problem  con¬ 
sidered  is  as  follows. 


Given  f(*):  E  R*  for  some  E  C  B,  what  conditions  on  E  and  f(*)  are 

.  .  *  _  *  *  . 
sufficient  to  guarantee  the  existence  of  a  set  S  C  E  such  that  x  e  S  im- 

plies  f(x  )  £  f(x)  for  all  x  e  E?  First,  define  the  following  constructs. 

The  functional  f(*)  is  lower  semi continuous  (Isc)  on  the  set  E  if  [7-3,  7-7] 
the  set  {x  e  Ejf(x)  >  X}  is  open  for  any  real  X  or  the  set  (x  e  Ejf(x)  X}  is 
closed  for  any  real  X.  Alternatively,  define  lsc  functions  sequentially  [7-5, 

7-7,  7-8]  (i.e.,  f(*)  is  lsc  on  E  if  for  any  sequence  {x^}  e  E  with  xr  5  x,  it 

follows  that  f(x)  <  lira  inf  f(x  )).  Strictly  speaking,  this  latter  deacrip- 
—  n+oo  n 

tion  defines  a  sequentially  lower  semicontinuous  function,  but  it  is  noted 
in  Reference  7-3  that  while  considering  strong  convergence,  the  two  concepts 
are  equivalent  [c.f.  7-7].  This  is  not  the  case,  however,  if  weak  convergence 

is  considered.  That  is,  we  can  define  f(*)  as  a  weakly  lower  semicontinuous 

(wise)  function  on  E  if  for  any  real  X,  the  3et  {x  e  E|f(x)  _<  X)  is  weakly 

closed,  and  define  f(*)  as  a  weakly  sequentially  lower  semicontinuous  (wslsc) 

function  on  E  if  for  any  {x  }  e  E  with  x  $  x  it  follows  that  f(x)  <  lim  inf 

n  n  **■“  n+°o 

f(xn).  It  is  noted  in  Reference  7-3  that  weak  sequential  lower  semicontinuity 

is  generally  a  weaker  property  than  weak  lower  semicontinuity,  i.e.,  if  f(») 
is  wise,  it  is  wslsc,  but  the  converse  is  not  generally  true.  . 

Similarly,  as  in  the  case  of  strongly  and  weakly  closed  sets,  it  is 
important  to  note  that  weak  lower  semicontinuity  is  a  "stronger"  concept  than 
(strong)  lower  semicontinuity.  Specifically,  if  f ( • )  is  wise  on  E,  then  the 

set  {x  e  E|f(x)  X)  is  weakly  closed,  hence  closed  for  all  real  X,  which 

implies  that  f(* )  is  lsc  on  E.  Again,  as  in  the  case  of  strongly  and  weakly 
closed  sets,  the  converse  does  hold  if  additional  conditions  are  imposed  on 
f(*).  Specifically,  the  functional  f(«)  is  quasiconvex  [7-3]  on  the  convex 

set  C  if  the  set  {x  e  c|f(x)  _<  \}  is  convex  for  all  real  X.  Thus,  if  E  is 
convex  and  weakly  closed  and  f(*)  is  lsc  and  quasiconvex  on  E,  the  set 
{x  e  E  f(x)  <  X}  is  (strongly)  closed  and  convex,  hence  weakly  .closed, 
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implying  £(•)  is  wise  on  E.  Note  that  the  notion  of  a  quasiconvex  functional 
is  a  generalization  of  the  more  common  notion  of  a  convex  functional,  which  is 
defined  on  a  convex  set  C  as  one  for  which 

f(Xxj  +  (1  -  X)x2)  <  Xf(Xl)  +  (1  -  X)f(x2) 

for  any  x^,  x2  e  C  and  0  <  X  <_  1.  It  is  noted  in  Reference  7-3  that  a  quasi¬ 
convex  functional  is  one  for  which 

f(XXj  +  (1  -  X)x2>  <  tnaxIfCx^  ),f(x2>] 

for  any  Xj,  x^  e  C  and  0  _<  X  £  1.  Since 

Xf(xj)  +  (1  -  X)f(x2>  <  maxlfCxj ),f(x2>] 

for  all  0  j<  X  <_  1,  it  is  clear  that  f(*)  is  qucsiconvex  if  it  is  convex. 

Given  the  definition  of  weak  sequential  lower  semicontinuity  just  pre¬ 
sented  and  the  notions  of  weak  convergence  and  weak  compactness  presented  in 

.  .  .  * 
Section  7.2,  it  is  simple  to  establish  the  following  existence  result  for  S  . 

Theorem  7-1  ([7-3],  Theorem  1.4.1) 

If  f(*):  E  -*■  is  wslsc  on  the  weakly  compact  set  E,  there  exists  a 

'ft  'k  k 

nonempty  set  S  C  E  such  that  x  t  S  implies  f(x  )  *  inf  f(x). 

xeE 

Proof 

Let 

f  *  inf  f(x) 
xeE 

By  definition  of  inf  f(x),  there  exists  a  sequence  {xn)  C  E  such  that 

lim  f(x  )  *  f* 

n+»  n 
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i 


K 


I 


•Jf 

Since  E  is  w-compact,  there  exists  some  x  e  E  and  a  subsequence  {x  }  of  {x^} 

"k 

such  that  x  ^  x  .  Clearly,  since  f(x  )  f  ,  £(x  )  ■*■  f  also.  Thus, 

\  n  \ 

since  f(»)  is  wslsc 


x  +  x  =>  f(x  )  <  lim  inf  f(x  ) 

nv  “  n*00  nk 


*  ,  *. 
f  <  f(x  ) 


Thus,  x  e  S  .  ■ 

Unfortunately,  in  the  applications  considered  here,  the  set  E  (corre¬ 
sponding  to  the  set  T  of  admissible  compensators  in  the  optimal  compensator 
design  problem)  generally  is  not  bounded,  hence  net  w-compact.  If  additional 

constraints  are  imposed  on  f(*),  however,  it  is  possible  to  guarantee  that 

f(*)  exhibits  a  set  of  global  minima  on  some  w-compact  subset  E'  of  E.  Spe¬ 
cifically,  f(*)  satisfies  a  T-property  on  E  if  there  exist  x^  e  E  and  Tq  >  0 

such  that  | | x  -  x^|j  T^  implies  f(x)  >  f(Xg)  for  any  x  e  E.  As  the  follow¬ 

ing  theorem  shows,  if  f(*)  satisfies  a  T-property  on  a  w-closed  set  E,  the  op¬ 
timization  can  be  restricted  to  some  w-compact  subset  of  E. 

Theorem  7-2  (adapted  from  [7-3],  Theorem  1.4.2) 

Suppose  E  is  a  weakly  closed  subset  of  a  reflexive  Banach  space  B.  If 

f(*)t  E  R*  is  wslsc  and  satisfies  a  T-property  on  E,  then  there  exists  a 

+  *  * 

w-compact  subset  C  of  E  and  a  set  S  C  C  such  that  x  e  S  implies 


f(x  )  =  inf  f(x) 

xeE 


Proof 

Let  S(Xq,Tq)  be  a  strongly  closed  sphere  of  radius  Tq  about  x^  s  E. 

Since  the  set  is  strongly  closed  and  convex,  it  ic  weakly  closed  and  bounded, 
hence  weakly  compact. 

Since  E  is  weakly  closed,  the  set 

C  -  { x  e  E |  1 1 x  -  xQ 1 1  <  TQ)  “  S(x0,T0)OE 

is  also  weakly  compact. 
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Since  £(•) 
x  £  C.  Thus 


satisfies  the  T-property, 


f(x)  >  f(xg)  for  any  x  e  E  but 


inf  f(x) 
xeE 


min 


inf  f(x),  inf  f(x) 
xeC  xftC 

xeE 


=  inf  f(x) 
xeC 


•k 

Since  C  is  w-compact,  the  existence  of  a  set  S 
f(*)  over  C  is  guaranteed  by  Theorem  7-1. 


of  global  rainiraizers  of 

£ 


It  is  possible  to  establish  that  f(*)  satisfies  a  T-property  in  a 
variety  of  ways.  One  of  the  simplest  is  to  guarantee  that  f(x)  increases  as 

IMI  increases.  For  example,  f(*)  is  semicoercive  on  E  if  [7-7] 


f  (x) 

INI 


0 


and  coercive  if 


+w 


from  which  it  is  immediately  clear  that  any  coercive  function  is  semicoercive. 
Lemma  7-1 


If  f(»):  E  R*  is  semicoercive  on  E,  it  satisfies  a  T-property. 

Proof 

Since  £(•)  is  semicoercive,  there  exists  some  a  >  0  such  that 


a 


Thus,  there  is  some  constant  M  >  0  such  that 


|  jxj  j  )  M,  x  £  E  *> 


f(x) 

INI 


>  I 
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“>  f(x)  >  f  j  |xj  | 


Now,  choose  any  x^  £  t  and  define 


M' 


max 


2f(x  ) 

H,  — - +  1 


so  that 


HI  >  *e  E  *> 


£(*) 


f(xQ)  + 
>  f(xQ) 


Finally,  note  that,  by  the  triangle  inequality 

||*- -oil  <  INI +  ll-oll  ■  INI +  INI 

->  ||*|l  >  ||*  -  »oll  •  ll*oll 


Thus,  if 


To  ■  ll*oll+M' 


it  follows  that 


1 1  *  -  *0 1 1  i  To  I  Ml  -  M’  =>  £(x)  >  £<Xc,) 
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Once  it  has  been  established  that  some  set  S  of  global  minimizers  of 
f(*)  over  the  weakly  closed  set  E  exist  on  some  weakly  compact  set  C,  we  can 

* 

turn  to  the  problem  of  how  elements  of  the  set  S  are  to  be  constructed.  Spe¬ 
cifically,  consider  the  sequence  { f^( • ) }  of  approximating  functionals  defined 

on  the  approximating  sets  {E^}  for  n  *  1,...  The  basic  question  considered 

here  is  "what  relationship  must  exist  between  the  approximating  sequences 
^  •$( 

{f  (•)}  and  {E  }  and  f(*)  and  E  to  guarantee  that  solutions  x  of  the  approxi- 
n  n  n 

mating  problems  defined  by  f^C*)  and  E^  converge  in  some  useful  sense  to  an 

it  ^ 

element  of  S  ?"  A  very  useful  answer  to  this  question  is  developed  in  Refer¬ 
ence  7-3  by  introducing  the  concept  of  a  "consistent  discretization"  of  the 
original  problem.  The  remainder  of  this  section  is  therefore  devoted  to 
presenting  a  slight  generalization  of  these  results  as  the  basis  for  the 
results  that  follow.  (The  generalization  is  this:  Reference  7-3  assumes  that 
the  set  E  over  which  f(*)  is  defined  and  minimized  is  the  entire  reflexive 
Banach  spa.ee  B.  In  the  problem  considered  here,  E  may  be  any  weakly  closed 
subset  of  B.) 


A  discretization  of  the  minimization  problem  of  f(*)  over  C  may  be  de¬ 
fined  [7-3]  as  a  family  of  subsets  {E^}  of  reflexive  Banach  spaces  {Bn) ,  a 

family  of  functionals  {f  (•)}  defined  on  E  ,  a  family  of  mappings  (p  }  of  E 

n  n  rr  °  rn  n 

into  E,  a  family  of  mappings  (r  }  of  E  into  E  ,  and  a  family  of  subsets  {C  } 

'ft  n  n 

of  E  .  The  basic  idea  is  that  {f  (•)),  {E  },  {B  },  and  { C  }  represent  ap- 
vi  n  n  n  *  n 

pros  itaat  ions,  in  some  sense,  of  the  components  f(*),  E,  B,  and  C  of  the  orig¬ 
inal  minimization  problem,  respectively.  The  mappings  {p^}  and  {r^}  thus 

represent  projections  of  the  approximating  sets  {E  }  onto  the  set  E  and  the 

n 

set  E  the  approximating  sets  {E^},  respectively.  These  definitions  are  gener¬ 
al  enough  to  allow  the  approximating  problems  and  the  original  problem  to  be 
defined  on  Banach  spaces  {B^}  and  B  that  are  qualitatively  very  different. 

In  particular,  the  approximating  spaces  {B^}  may  be  finite  dimensional. 

Using  these  definitions,  Reference  7-3  associates  the  following  approxi¬ 
mate  problem  with  each  discretization.  Given  a  sequence  of  positive 

numbers  converging  to  zero,  consider  the  set  i  of  elements  x  e  C  such  that 

n  n  n 


f  (x  ) 
n  n 


<  inf  f  (x  )  +  £ 
-  x  eC  n  n  n 
n  n 


Note  that  if  there  exists  some  subset  S  of  C  such  that 

n  n 
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f  (x  )  =  inf  f  (x  ) 

n  n  xneCn  n  n 


for  all  x  e  S  ,  then  it  follows  that 
n  n 


f  (x  )  =  inf  f  (x  )  <  inf  f  (x  )  +  5 
nn  xeC  nn  “xeC  nn  "n 

n  n  n  n 


for  any  £n  >  0,  implying  Sn  C  Sn. 

In  order  to  develop  sufficient  conditions  for  the  set  of  to  converge 

*  .... 
in  some  sense  to  S  ,  Reference  7-3  defines  a  consistent  discretization  as  a 

discretization  of  the  original  optimization  problem  that  satisfies  the  follow¬ 
ing  constraints 


(1) 

The 

set  S  exists 
n 

(2) 

lim 

n+“ 

sup  lf(pnxn) 

(3) 

The 

sets  C°  =  p 

n. 

(4) 

If  : 

z  e  C  and 

T> 

ll  • 

1 

(5) 

For 

*  * 
some  x  e  S  , 

(6) 

For 

* 

the  same  x  e 

n  n 


* 

r  x  e  C  . 
n  n 


Note  that  if  p^C^  C  C  for  all  n,  Conditions  (3)  and  (4)  are  satisfied  auto¬ 
matically.  These  conditions  are  sufficient  to  establish  the  following  result. 

Theorem  7-3  ([7-3],  Theorem  2.2.1) 

Suppose  (f  (• ),  E  ,  B  ,  C  ,  p  ,  r  }  define  a  consistent  discretization, 
rr  n  n*  n  n  rn  n 

and  let  K  be  a  subset  of  E  such  that  Cn  C  K  for  all  n  sufficiently  large.  If 
C  is  weakly  compact  and  x^  e  for  all  n,  then 
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lim  f(p  3  )  =  lim  f  (x  )  =  f(x  ) 

21*00  n  n  n*oo  n  n 

** 

where  x  satisfies  Conditions  (5)  and  (6).  Further,  the  sequence  {p^x^}  ®ay 
be  decomposed  into  subsequences,  each  of  which  converges  weakly  to  an  element 

of  S  ,  i.e. ,  p  S  ?  S  CS  . 

*  n 

Before  proceeding  with  the  proof  of  Theorem  7-3,  the  following  lemma, 
required  in  the  proof,  is  established. 

Lemma  7-2 

The  sequence  converges  to  zero  where 

Y  =  f(x  )  -  inf  f(x) 
n  „n 

xeC 

*  * 

and  x  e  S  .  | 

Proof  of  Lemma  7-2 

•flp  j. 

(1)  First,  note  that  since  x  e  CCC  ,  ^  >  0  for  all  n. 

(2)  Next,  for  each  n,  select  some  z n  e  Cn  such  that 

f(z  )  <  inf  f(x)  +  —  a  f(x  )  +  [-*-  Y  ] 

n  —  n  n  ln  nJ 

xeC 

(3)  Now,  choose  N  such  that  n  >.  N  “>  Cn  C  K  C  B  and  note  that  since  K 
is  bounded,  the  sequence  {z^}  is  bounded  for  all  n  ^  N.  Since  B 
is  a  reflexive  Banach  space,  (zn)  rosy  be  decomposed  into  a  collec¬ 
tion  of  subsequences  {z  },  each  of  which  converges  weakly  to  some 

„  nk 

z  e  B. 
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(4)  3y  Condition  (4)  o£  the  consistency  requirements 


z  ^  z  ->  z  t  C 


so  that  f(x  )  <  £(z). 


(5)  Note  also  that  since  £(•)  is  wslsc  on  K 


f(z)  <  lira  inf  f(z  ) 


k+“ 


<  f(x  )  +  lira  inf  [— - y  ] 

« _  in  J 


k+“ 


=>  lim  inf  [- - y  ]  >  0 

k-M»  n,  n,  J 


(6)  Thus 


lim  inf  [- - y  ]  >  lim  inf  —  -  lim  sup  y 

k-voo  nt  —  k*°°  nv  k+»  nt, 


k+» 


">  lim  inf  —  >  lim  sup  y 

nk  “  \ 


~>  0  >  lira  sup  y 

—  k+»  \ 


>  lim  inf  y  >  0 

“  k*00  \  ~ 


=>  lim  y  =0 
k+co 


(7)  Since  this  proof  may  be  repeated  for  all  subsequences  {z^  }  de¬ 
fined  in  Item  3,  it  follows  that  lim  y  **  0,  ^ 
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Proof  of  Theorem  7-3 
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Since  lim  £  *  lim  y  *  0  and,  by  consistency  Conditions  (2) 

n  n-*»  n 

and  (5),  respectively,  it  follows  that 

lim  sup  n  <0 
n+“  n  ” 


and 


lim.  sup  6  <  0 

n  “■ 


(3)  Thus 


lim  (6  +  n  )  “  0 

q->oo  n  n 

Note  that  if  lim  inf  n  <  a  for  some  a  <  0,  there  exists  some  sub- 

n+«>  n 

sequence  {q  }  so  that 

i 

n  <  a/2  <  0 

n .  — 


for  all  i,  implying  there  is  some  subsequence  {6  }  such  that 

i 


6  >  -a/4  >  0 

n.  — 


Howevei ,  this  implies 


lira  sup  6  >  -a/4  >  0 

n+°°  n  — 


which  contradicts  Condition  (5)  of  the  consistency  requirements. 
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Thus 


lim  inf  n  >  0  >  lim  sup  n  *>  lim  n  *  0 

jv*»  n  —  n  n+oo  n 


Similarly,  it  follows  that 


lim  6  3  0 

n+°°  n 


(4)  Thus,  it  fellows  from  the  inequalities  in  Item  (1)  that 


f(x  )  *  lim  f (p  x  )  3  lim  f  (x  ) 

n->oo  n  n  Q-foo  n  n 


(5)  Finally,  since  Cn  C  K  for  all  n  >_  N  with  K  bounded,  the  sequence 
* 

(p^x^)  may  be  decomposed  into  subsequences  that  converge  weakly  to 
limit  points  in  B. 

Since  x  e  C  ,  p  x  e  Cn  for  ail  n  so  that  by  consistency  Condi- 
n  n  n  n  J 

tion  (4),  if  p  x  ^  z,  then  z  e  C. 

Further,  since  f(»)  is  wslsc  on  K 

f(z)  <  lim  inf  f(p  x  )  *  f(x  ) 


* 

so  z  e  S  . 


Thus  p  S 
rn  n 


Application  to  Fixed-Form  Compensator  Design  Problem 


The  results  just  presented  are  very  general;  thus,  they  are  applicable 
to  a  very  broad  class  of  problems  including  the  fixed-form  compensator  design 
problem  of  interest  herein.  Because  of  their  generality,  the  conditions  these 
results  impose  on  arbitrary  optimization  problems  are  somewhat  hard  to  inter¬ 
pret.  However,  if  they  are  specialized  to  the  optimal  fixed-form  compensator 
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design  problem,  these  conditions  simplify  considerably.  Consequently,  the 
remainder  of  this  section  is  concerned  with  the  following  special  case. 

The  arbitrary  functional  f(*)  will  be  replaced  by  the  cost  functional 
J(* )  defined  over  some  set  T  of  admissible  compensator  designs.  It  is  impor¬ 
tant  to  note  that  this  is  a  departure  from  traditional  optimal  control  prob¬ 

lems  in  which  J(*)  is  a  functional  on  some  function  space  U  of  control  inputs 
u(t).  While  it  is  true  that  a  priori  restriction  of  the  control  law  to  some 

prespecified  form  may  result  in  a  performance  degradation  (i.e.,  the  set  of  T 

of  admissible  compensators  may  not  be  large  enough  to  generate  the  control  in- 

put  u  (t)  that  minimizes  J(u)  over  U),  it  is  also  true  that  practical  designs 
will  be  restricted  to  some  set  of  admissible  compensators  anyway.  Thus,  the 
approach  developed  is  to  identify  the  set  T  a  priori  on  the  basis  of  various 
important  practical  considerations  and  to  optimize  performance  within  these 
restrictions.  In  particular,  it  is  assumed  that  once  these  restrictions  have 
been  identified,  a  compensator  form  meeting  these  restrictions  is  selected 
that  reduces  the  design  problem  to  that  of  selecting  values  for  some  finite 
number  i  of  design  parameters  that  minimize  the  cost  J(*).  Thus,  the  set  T  of 

l 

admissible  compensator  designs  becomes  a  subset  of  the  parameter  space  R  . 

The  most  significant  simplification  inherent  in  this  formulation  of  the 
compensator  design  problem  is  that  optimization  is  now  carried  out  over  the 

£ 

finite-dimensional  Euclidean  space  R  rather  than  the  infinite-dimensional 
space  U.  Consequently,  since  weak  and  strong  convergence  are  equivalent  in 
finite-dimensional  Banach  spaces,  the  results  of  the  preceding  section  can  be 
used  to  obtain  sufficient  conditions  for  the  strong  convergence  of  the  ap¬ 
proximate  compensator  designs  to  the  desired  optimal  limit.  That  is,  by  for¬ 
mulating  the  design  procedure  described  in  the  introduction  as  a  discretiza¬ 
tion  of  the  minimization  of  J(«)  over  T,  sufficient  conditions  for  the  ccn- 

★ 

vergence  of  the  optimal  parameter  vectors  k^  for  the  approximating  plants  Pn 

.  *  .  .  .  .  j 

to  the  optimal  parameter  vectors  jk  for  the  infinite-dimensional  plant  P  may 
be  obtained. 


In  particular,  the  discretization  defined  by  the  design  procedures  is 
the  following.  First,  all  minimizations  will  be  restricted  to  the  closed  set 

i  i 

T  C  R  (i.e.,  B  *  B  a  R  and  E  *  E  *  T  for  all  n  in  the  notation  of  the  pre- 

vious  subsection).  Consequently,  the  mappings  of  {p^}  and  {r^}  appearing  in 

the  definition  of  a  discretization  are  all  identity  mappings  in  this  case. 

The  approximate  optimization  problems  associated  with  this  discretization  will 

•  .  .  .  /  \  / .  ^  .  ch 

be  the  minimization  of  J  (•)  over  T  (i.e,,  the  set  S  of  solutions  to  the  n 

n  *  n 

approximating  problem  consists  of  k  e  F  such  that  J  (k  )  =  inf  J  (k).  In 

— n  n  — n  ^ef  n  — 

order  for  the  design  problem  considered  here  to  be  well  posed,  it  is,  of 

course,  necessary  that  the  sets  {S  }  exist  for  all  n.  If  J  (•)  and  J(»)  are 

n  n 

required  to  satisfy  the  conditions  of  Theorem  7-3,  for  example,  the  existence 
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of  the  sets  { S  }  and  S  is  guaranteed  as  are  the  existence  of  the  compact  sets 

if 

{(M  and  C  containing  (Sn)  and  S  ,  respectively.  Note,  however,  that  without 

further  restrictions  on  the  functions  {J  (•)},  the  compact  sets  {C  }  need  not 

n  r  n 

bear  any  special  relation  to  the  set  C.  In  particular,  it  is  entirely  pos¬ 
sible  that  the  sets  {C  }  are  not  uniformly  bounded  and  will  therefore  not  con- 

* 

verge  to  the  set  C  in  any  sense  as  n  +  »,  which  suggests  that  the  sets  {S  } 

*  " 
may  not  converge  into  the  set  S  as  n  ■*•  «.  If,  however,  it  can  be  established 

that  the  discretization  just  defined  satisfies  the  consistency  requirements 
developed  in  the  preceding  subsection,  it  will  follow  from  Theorem  7-3  that 
* 

the  sets  of  optimal  designs  associated  with  the  plants  Pn  will  converge  to 

* 

to  a  subset  of  the  set  S  of  optimal  designs  associated  with  the  plant  P  and 
.  *  *  *  *  *  * 
that  J^k^)  *  J(k  )  for  any  k^  z  Sn  and  k  e  S  .  Specifically,  combining  the 

definition  of  a  consistent  discretization  with  the  conditions  of  Theorem  7-3 
and  specializing  them  to  the  present  problem,  the  following  result  is 
immediate. 

Theorem  7-4 

£  1 

Suppose  T  is  a  closed  set  in  R  and  J(*):  T  R  is  an  lsc  function  ex- 
.  .  * 

hibitmg  some  bounded  set  S  of  global  minima  over  T.  Suppose  further  that 
each  function  Jn<*):  T  -*■  R*  also  exhibits  some  bounded  set  S^  of  global  minima 
over  T  and  that  the  following  conditions  are  met: 


if  if 

(1)  S^  C  K|  S  CK  for  some  compact  set  K  CT. 


(2)  lim  sup  [j(k  )  -  J  (k  )]  <  0  if  k  eS. 
n-ftw  *— n  n  n  •— n  u 


.  ,  ■  K.  ,  *  * 

13)  lim  sup  J  ik  )  <  Jvk  )  for  some  k  e  S  . 
rv+oo  n  —  —  — 


Under  these  conditions,  the  following  results  hold: 


(A)  S*  +  S  C  S*. 
n 


(5)  lim  J  (k  )  =  lim  J(k  )  ■  j(k*) 

tv>«  n  —o  n+00  — -n  — 


if  if  k  ie 

where  k^  e  and  k  e  S  satisfies  Condition  (3). 
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Note  here  that  the  compact  set  K  defined  in  Condition  (1)  takes  the 
place  of  the  sets  {Cn>  and  C  in  the  original  definition  of  a  consistent  dis¬ 
cretization.  This  guarantees  that  Conditions  (3),  (4),  and  (5)  of  the  origin 

*  * 

al  definition  are  satisfied.  Given  the  assumed  boundedness  of  sets  S  and 
individually,  it  is  clear  that  Condition  (1)  of  Theorem  7—4  can  be  satified  if 
and  only  if  the  sets  {S^}  are  uniformly  bounded. 

If  Condition  (1)  of  Theorem  7-4  is  satisfied,  Conditions  (2)  and  (3)  may 
be  replaced  by  a  uniform  convergence  condition,  as  the  following  lemma 
demonstrates. 

Lemma  7-3 

If  Condition  (1)  of  Theorem  7-4  is  satisfied  and  if  J^C*)  +  J(*)  uni¬ 
formly  on  the  set  K,  then  Conditions  (2)  and  (3)  are  also  satisfied. 

Proof 

(1)  Condition  (3)  follows  immediately  since  if  Jn(*)  J(‘)  uniformly 

on  K,  k*  e  S*C  K  ->  Jn(k*)  ♦  J(k*>. 

Thus 

ic  *  "k 

lira  sup  J  (k  )  *  lim  J  (k  )  *  J(k  ) 

n  n+"  n  — 

(2)  Similarly,  since  Jn(*)  ♦  J(*)  uniformly  on  K,  J(k)  -  Jn(k)  *  0 
uniformly  for  all  k  e  K. 

In  particular,  since  k^  e  C  K  for  all  n,  JCk^)  “  +  ®  as 

n  “. 

Thus 

lira  8 up  [J(k  )  -  J  (k  )]  *  lim  [J(k  )  -  J  (k  )]  “  0 

n-Mx>  r  — n  n  — n  n+«  — n  n  n 

and  Condition  (2)  is  satisfied,  ® 
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To  establish  Condition  (1),  either  one  of  two  approaches  may  be  taken, 
resulting  in  either  a  constrained  optimization  problem  or  an  unconstrained 
optimization  problem.  In  the  first  case,  if  J(>)  is  required  to  satisfy  a 
T-property,  then  it  follows  from  Theorem  7-2  that  there  exists  some  compact 

•ic 

set  C  containing  the  set  S  of  all  global  minimizers  of  J(0  over  T.  Conse¬ 
quently,  if  each  function  )  is  lower  semicontinuous  on  T,  the  existence 

** 

of  nonempty  sets  Sn  C  C  that  minimize  J^i*)  over  C  is  guaranteed  by  Theorem 

•kit 

7-1.  Thus,  applying  Theorem  7-4  to  the  sets  SR  leads  immediately  to  Theo¬ 
rem  7-5. 

Theorem  7-5  (Constrained  Optimization) 

i  1 

Suppose  T  is  a  closed  set  in  R  and  J(*):  r  R  is  an  lsc  function 

* 

satisfying  a  T-property  on  T.  This  guarantees  that  S  C  C  for  some  compact 
set  C.  Suppose  {Jn(*)}  is  a  sequence  of  lsc  functions  on  C  that  converges 

'k’fc 

uniformly  to  j(*)  on  C,  and  let  S^  denote  the  set  of  global  minima  of  <1^(0 
over  C,  guaranteed  to  exist  by  Theorem  7-1. 

Under  these  conditions,  the  following  results  hold: 

-t.-l.  A  X 

,  .  XX  _  X 

(1)  S  +  S  C  S  as  n  ♦  * 

n 

(2)  lira  J  (k  )  *  lim  J(k  )  -  J(k  ) 


**  **  *  * 
where  k  t  S  and  k  e  S  . 
— n  n  — 


It  is  extremely  important  to  note  here  that  the  sets  (Sn  }  are  not  the 

•ff 

the  same  as  the  sets  {Sn}  considered  earlier.  For  example,  suppose  J(*) 

satisfies  the  conditions  of  Theorem  7-!>  and  inf  J(k)  >  0.  Then  if  the  ap- 

ker  “ 

proximating  functions  {J  (•)}  are  defined  by 

n 


J  (k) 
n 


)  J(k) 

( ° 


INI  <  " 
INI  >  * 
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for  all  k  £  T,  it  follows  that  for  a  such  that  C  =  (k  e  rj  ||k||  <  n} 

•kle  *  Irk  1c 

Sn  =  S  ,  guaranteeing  that  Sn  S  as  n  “  in  agreement  with  Theorem  7-5. 

if  ~ 

The  sets  S^,  however,  consist  of  all  elements  Vc  of  T  not  in  L^,  ensuring  that 


J  (k)  *  0  =  inf  J  (k) 

n  -  kef  n  ” 


Note  that  J  (k)  is  lsc  on  T  since  J(k)  is  lsc  on  1.  ,  and  0  is  obviously  lsc 
n  —  —  n 

on  T  -  L^.  Thus,  the  only  points  in  question  occur  on  the  boundary  ft  between 

L  and  T  -  L  .  Since  L  is  open,  T  -  L  is  closed  so  ft  C  T  -  L  implying 
n  n  n*^*n  n  r  J 


J(k)  =  0  <  lim  inf  J(k  ) 

**  n+°o 

if  k  e  ft  for  any  sequence  {k  }  in  T.  Note  also  for  any  compact  set  K  C  T, 

there  is  some  N  such  that  J  (k)  “  J(k)  for  all  k  e  K  if  n  >  N.  Thus,  {J  (*)} 

n  —  —  —  —  n 

converges  uniformly  to  J(*)  on  any  compact  subset  of  T. 

To  avoid  these  difficulties,  the  second  approach  to  establishing  that 
Condition  (1)  of  Theorem  7-4  holds  is  to  strengthen  the  conditions  on 

and  J(*)  as  follows.  First,  uniform  convergence  will  be  required  on  all  of  T 
rather  than  just  compact  subsets  of  T.  This  requirement  will  exclude  the 
pathologies  exhibited  by  the  last  example,  although  even  under  this  restric¬ 
tion,  auother  difficulty  can  arise.  In  particular,  the  fact  that  J(»)  satis¬ 
fies  a  T-property  is  not  enough  to  guarantee  that  J(k)  does  not  approach 

inf  J(k)  as  j  |k|  |  ®.  For  example,  consider  the  one  parameter  cost  function 

ke  F 

illustrated  in  Figure  7-1.  Here,  J(k)  satisfies  a  T-property  at  k  *  0,  but 
J(k)  J (0 )  as  k  ■*■  ±«,  To  correct  this  difficulty,  define  a  sharp  T-property 
as  follows. 


Definition  (Sharp  T-Property) 


The  function  J( * ) :  r  ■*  R  satisfies  a  sharp  T-property  if  there  exist 
kqjkj  G  **  atu*  tq  ^  ®  8UC^  t'iat 


Ilk-kJI  >.  Tq  *>  J(k)  >  JCkj)  >  J(k2) 
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Figure  7-1.  J(k)  satisfies  a  T-property  at  k  a  0  but 
J(k)  inf  J(k)  as  k  ±« 

keR1 

It  is  clear  from  this  definition  that  any  J(*)  satisfying  a  sharp  T- 
property  satisfies  a  T-property  and  that  if 

L  -  (l.I  |  ||*  -£,11  <  V 


L  is  a  compact  set  and 

inf  J(k)  >  J(k, )  >  J(k.)  >  inf  J(k) 
k*L  -  -  -1  -2  -  keL  - 

It  is  also  clear  that  if  J(»)  is  a  semicoercive  function,  it  satisfies  a  sharp 
T-property.  With  these  results,  the  following  unconstrained  optimization 
theorem  may  be  established. 

Theorem  7-6  (Unconstrained  Optimization) 

%  1 

Suppose  r  is  a  closed  set  in  R  and  j(*):  T  ♦  R  is  an  lsc  function  sat¬ 
isfying  a  sharp  T-property  on  T.  If  {J  (•)}  is  a  family  of  lsc  functions  on  T 

that  converges  uniformly  to  J(»)  on  T,  then  the  following  conditions  hold: 

(1)  All  global  infiraa  of  -1^(0  and  J(»)  are  achieved  on  some  compact 

'k  k 

set  K,  i.e. ,  Sn,  S  C  K. 

(2)  S  S  C  S  as  n  ■»  ®, 

n 
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MsnsfeJVE  -  - 


(3)  lim  J  (k  )  =  lira  J(k  )  *  J(k  ) 

ft  —ft  n+«  “ii  “* 


*  *  *  * 
where  k  e  S  and  k  e  S  . 
— n  n  — 


Proof 


(1)  Since  J(*)  satisfies  a  sharp  T-property  on  T,  there  exists  some 
compact  set  K  such  that 


inf  J(k)  >  inf  J(k) 
k£K  “  keK  “ 


Further,  since  J(0  is  lsc  on  T,  by  Theorem  7-1  there  is  some  set 
* 

S  C  K  such  that 


*  *  ,  * 
k  e  S  =>  J(k  ) 


inf  J(k)  =  inf  J(k) 
keK  ~  ker  - 


(2)  Let 


d  *  inf  J(k)  -  inf  J(k)  >  0 

k#K  -  keK  ~ 


Since  +  *(•)  uniformly  on  T,  there  is  some  N  such  that  if 

n  >  N 


J  (k)  -  J(k)|  £  d/3  for  all  k  e  Y 


Thus,  if  k  e  K 


Jn(k)  >  J(k)  -  d/3 


>  inf  J(k)  -  d/3 
“  k£K  *“ 


>  inf  J(k)  +  2d/3 
~  keK  “ 


J(k  )  +  2d/3  for  some  k  e  S 
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Sirai Lar ly 


(3) 


J  (k*)  <  J(k*)  +  d/3 

n  —  —  — 


*>  J  (k  )  <  J  (k)  for  any  k  i  K 

n  —  n  —  — 


Thus 


inf  J  (k)  =  inf  J  (k) 

ker  n  keK  n 


for  all  n  >  N. 


Since  J^C*)  is  lsc  and  K  is  compact,  by  Theorem  7-1,  exists 
such  that 


k*  e  S*  *>  J  (k*)  »  inf  J  (k)  *  inf  J  (k) 

n  n  n  -n  keK  n  -  kef 


ic 

Now,  since  sn  C  K  and  S  C  K,  K  compact,  Condition  (1)  of  Theorem 
7-4  is  satisfied. 

Since  J^C*)  *  J(*)  uniformly  on  T,  it  follows  from  Lemma  7-3  that 
Conditions  (2)  and  (3)  of  Theorem  7-4  are  satisfied  also. 

Consequently,  it  follows  from  Theorem  7-4  that 


sn  *  scs’ 


and  that 


lim  J  (k  )  = 


n+» 


n  — n 


lim  J(k  )  *  J(k*) 


-n 


*  *  *  * 
where  k  e  S  and  k  e  S  . 
— n  n  — 
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Since  it  was  the  difficulty  of  dealing  directly  with  the  cost  function 
J(*)  that  motivated  the  design  procedure  considered  here  originally,  it  is  de¬ 
sirable  to  have  an  optimality  theorem  that  establishes  the  results  of  Theo-. 
rem  7-6  on  the  basis  of  the  properties  of  the  approximating  functionals 
(Jn(*)}.  Because  of  the  uniform  convergence  requirement  imposed  on  the  se¬ 
quence  {Jn(*)},  it  is  not  too  difficult  to  establish  such  a  theorem,  on  the 
basis  of  the  following  lemma. 

Lemma  7-4 

If  J  (• )  is  lsc  on  the  closed  set  T  for  all  n  and  the  sequence  {J  (•)} 
n  n 

converges  to  J(* )  uniformly  on  compact  subsets  of  T,  then  J(»)  is  lsc. 

Proof 

(1)  To  prove  J(*)  is  lsc  on  T,  it  must  be  shown  that  given  any  It  e  T 
and  any  C  >  0,  there  exists  q  >  0  such  that 

A  A 

k  e  r,  |  |k  -  k|  |  <  n  ->  J(k)  _>  J(k)  -  C 

(2)  To  accomplish  this,  consider  the  set 

K  *  {k  e  r|  | |k  -  k| |  £  1} 

-  r  0{k|  ||k  -  k||  <  1} 

which  is  compact  since  T  is  closed, 

(3)  Thus,  J  (• )  +  J(»)  uniformly  on  K  so  that  there  is  ’some  integer  N 

such  that  n  N  =>  | Jn Ck_)  -  j(k)|  <  C/3  for  any  k  e  K 

(4)  Since  J  (• )  is  lsc  on  T,  hence  K,  there  is  some  6  >  0  such  that 

n  ’  ’ 

k  e  f,  | jk  -  k| |  <6  =>  JN<k)  >  JN<k)  -  C/3 
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(5)  Thus,  let  n  3  rain  (1,6)  so  that  if  k  e  T,  j  jh  -  kj  j  _<  n  then 
J(k)  >  JN(k)  -5/3  (by  (3) ) 

>  JN(k)  -  25/3  (by  (4)) 

2^  J(k)  -  5  (by  (3)) 

®>  J(0  is  lsc  on  T.  ■ 

Combining  the  results  of  this  lemma  with  those  of  Theorem  7-5  and  the 
definition  of  a  semicoercive  function  yields  the  following  optimality  theorem. 

Theorem  7-7 

If  the  following  conditions  are  met: 

(1)  T  is  a  closed  set. 

(2)  Jn(*)  *  J(‘)  uniformly  on  compact  subsets  of  T. 

(3)  J  (k)  >  Jrt(k)  for  all  k  e  T  for  some  semicoercive  function  J„(*). 

n  —  —  0  —  —  0 

(4)  J  (•)  is  lsc  for  all  n. 

n 

then 

(5)  All  global  infima  of  Jn(0  and  J(* )  occur  on  some  compact  set  K, 

ic  £ 

>  S  j  S  C  K» 
n 

if  A  if 

(7)  S  +  S  C  55  as  n  + 

n 

(8)  lim  J  (k  )  »  lira  J(k  )  3  J(k  )  where  k  e  S*,  k*  e  S. 

n  — n  n.>co  — n  —  — n  n  — 

Proof 

(1)  By  Lemma  7-4,  since  +  J(»)  uniformly  on  compact  subsets  of 

T  and  J  (•)  is  lsc  on  T,  J(‘)  is  lsc  on  T. 
r* 
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(2)  Since  Jg(*)  is  semicoercive 


u,  Vi> 

INI”  Hill  '  ° 

for  some  a  >  0.  In  particular,  there  is  some  B  >  0  such  that  j |lc| j 
>  S  «>  JQ(k)  >  a/2  | |k| | . 

(3)  Now,  pick  any  s  F  and  let 


B' 

1  \ 

-  max  ^B,  Q  ,  lliollj 

so  that  i 

£  Hill 

>  B  * 

then 

Vi> 

if 

Hill  > 

aB'  . 

2  - 

2J(ko> 

Next,  define  K  * 

{_k  e 

rl  Hill 

£  B'}  and 

note  that 

for  any  n 

inf  J  (k) 

m 

inf 

J  (k) 

>  inf 

J„(k)  > 

2J(k-) 

k*K  n  “ 

ll> 

sih 

n  — 

t 

"  llill>6 

0  —  — 

• 

—0 

Further, 

since  J 

t 

*  J(V- 

there  is 

some  integer  N  such  that 

n  >  N  *>  1  2  a>  Jn^0^  —  2 


Thus 


inf  J  (k)  >  J  (k.)  >  inf  J  (k) 

k^K  n  —  n  -0  -  kef  n  - 


so  that  S  C  K  for  all  n. 
n 
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(5)  Similarly,  since  J(*)  is  lsc  and  K  is  compact,  J(*)  achieves  its 

* 

global  infimura  over  K  on  some  set  S  C  K.  Also,  for  any  It  t  K 

J(k)  *  lim  J  (k)  >  2J(k.)  >  J(kJ  >  inf  J(k) 

-  n+c  n  -  -  -0  0  -  kef  - 


so  that  all  global  minima  of  J(*)  occur  on  K. 

(6)  Consequently,  Condition  (5)  is  established,  and  Conditions  (6)  and 
(?)  follow  from  Theorem  7-6,  B 

7,5  Summary  of  Convergence  Conditions 

The  results  presented  in  this  section,  culminating  in  Theorem  7-7,  pro¬ 
vide  sufficient  conditions  on  the  set  T  of  admissible  compensator  designs  and 
the  approximating  costs  {J^C*)}  to  guarantee  that  the  four  fundamental  ques¬ 
tions  raised  in  the  introduction  have  (approximately)  affirmative  answers. 
Specifically,  these  conditions  guarantee  the  following  results: 

(1)  A  set  S  of  optimal  compensator  designs  for  the  distributed  param¬ 
eter  plant  P  exists. 

* 

(2)  For  every  n,  a  set  Sr  of  optimal  compensator  designs  for  the  fi- 
n.  te  dimensional  approximating  plant  exists. 

ic 

(3)  The  sequence  {Sn}  of  sets  of  optimal  designs  for  the  plants  {P^} 

converges  to  a  set  S  of  optimal  designs  for  the  distributed  param¬ 
eter  plant  P. 

Several  important  points  should  be  noted  regarding  these  results.  First 
the  assumption  thnt  the  set  T  of  admissible  compensators  is  closed  is  essen¬ 
tial.  In  particular,  this  requirement  is  imposed  to  avoid  the  possibility 
that  inf  J(k)  exists  but  is  not  achieved  on  T.  A  simple  illustration  of  this 
keT 

is  the  following:  minimize  J(k)  s  -k  over  T  =  (-“,0).  Here 

inf  J(k)  *  0 

ke.r 


even  though  J(k)  >  0  for  all  k  e  T. 


This  example  is  unrealistic  in  certain  respects,  however,  since  it  is  to 
be  expected  that  the  definition  of  r  will  somehow  reflect  a  stabilizability 
requirement  and  that  J(k)  will  approach  •*»  as  jc  approaches  the  boundary  of  T. 

A  more  representative  example  is  therefore  the  following:  minimize 


J(k) 


k2 


over  T 


<  krk2  <  0) 


Here,  for  any  fixed  k^,  J(k)  -*■  +“  as  fc^  *  0“  and  similarly  for  fixed  k2»  J(k) 

+  +°°  as  kj  -*■  0“.  Further,  this  function  is  coercive,  suggesting  it  should 

achieve  its  infimura  c"'er  F  on  some  compact  set  CCr,  Note,  however,  that  for 
any  fixed  k2 


V 


J(k1,k2) 


k 


2 

2 


when  kj  3  k2<  Thus,  the  sequence  {k^}  *  { (-l/n,-l/n)}  yields  the  sequence  of 

2 

costs  { J(k  )}  *  {1/n  }  which  converges  to  0  as  n  However,  for  any  k  e  T, 

*" 

2 

J(k)  >  k„  >  0,  inf  J(k)  «  0  is  never  achieved  since  the  minimizing  sequence 
*  keT 

{(-1/n, -1/n)}  converges  to  a  limit  outride  i.  Consequently,  since  the  set  A 
of  stabilizing  compensators  meeting  any  particular  structural  constraint  is 
usually  an  open  set  (correspond-,.  ^  to  the  fact  that  the  root  locus  is  continu¬ 
ous),  it  is  clear  that  F  mvo*-  he  taken  as  some  closed  subset  of  A.  Physical¬ 
ly,  this  corresponds  to  ' <  ^position  of  some  minimum  stability  margin  or 
some  other  standard  of  /ccep  ability  like  that  considered  in  Reference  7-9. 
Another  possibility  woukU  \c  to  restrict  the  optimization  a  priori  by  placing 
explicit  bounds  on  each  component  k.  of  the  parameter  vector  k.  Such  a 

restriction  would  further  guarantee  that  the  set  T  is  compact,  which  taakes 
the  results  of  Theorem  7-5  directly  applicable. 

The  uniform  convergence  requirement  on  the  sequence  {J  (•)}  is  also  ex- 

n  *, 

tremely  important,  as  illustrated  by  the  example  in  Figure  7-2.  Here,  J^ik)  + 

J(k)  for  all  k  e  R*  but  k  =  2  and  j(k  )  3  1  while  k  3  l/2n  and  {J  (k  )}  *  0 

n  n  n 

for  all  n.  Thus,  the  sequences  {k  )  and  {j  (k  )}  both  converge  to  well 

n  n  n 

defined  limits  that  are  not  k  and  j(k  ),  respectively. 
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The  requirement  that  Jn(*)  b®  isc  for  all  n  is  noC  particularly  restric- 

tive.  Specifically,  in  most  practical  compensator  design  problems,  the  cost 
functionals  are  chosen  to  be  differentiable  so  that  gradient -derived  necessary 
conditions  for  optimality  or  gradient-based  computational  algorithms  may  be 
used  to  compute  the  optimal  compensator  parameters.  Consequently,  J^C*)  will 

normally  be  continuous  and  thus  lsc. 

Assuming  all  of  these  conditions  are  met,  it  is  important  to  note  that 
they  only  guarantee  the  convergence  of  the  sequence  of  sets  (Sn)  into  the  set 
*  *  .  * 

S  .  If  S  consists  of  a  unique  optimal  design  k  ,  then  any  sequence  {k^}  of 

elements  of  (S  }  is  guaranteed  to  converge  to  It  if  the  conditions  of  Theorems 

★ 

7-5,  7-'  or  7-7  are  met.  If,  however,  S  contains  more  than  one  optimal  de- 

it 

sign,  then  the  limit  of  the  sequence  {S^}  need  be  only  a  subset  of  S  .  This 
point  is  illustrated  by  the  following  example: 


J(k) 


k  <_  -1 
-1  <  k  £  0 
0  <  k  £  1 
k  >  1 


J  (k) 
n 


-k  ••  1  -  -  ,  k  <  -1 

n  — 

( 1  +  ~)  k  +  1  ,  -1  <  k  <  0 

i  -  ( !  +  i-)  k  ,  0<k<  1 
1 


k  -  1  - 


2n 


,  k  >  1 


Note 


that  only  (j^Ck)  -  J(k)|  £  1/n  for  all  k  e  r'  so  Jn(0  +  J(*)  uniformly 


on  R  ,  but  while  J(k)  exhibits  minima  at  k  *  ±1,  J  (k)  exhibits  a  unique 

n 

global  minimum  at  k  *  -1  for  all  n.  Thus,  *  {-1},  S  *  {-1,1},  so  that 

•k  *  *  it 

S  +  SCS  buu  S  *  S  . 
n 
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Finally,  note  that  the  convergence  conditions  developed  here  .are  only 
sufficient  and  not  necessary.  In  fact,  it  would  appear  to  be  very  difficult 
to  develop  any  necesary  conditions  for  convergence ,  as  the  following  example 
illustrates. 


Suppose  J(k)  =  |k|  and  take 

( 1  +  e  +  |k|  ,  ifn  odd 

(l-e+'j-j-)  |kj  ,  ifn  even 

for  some  e  >  0.  Clearly  both  J(k)  and  J^Ck)  exhibit  unique  global  minima  at 

k  =  0  for  which  J(k)  «*  J^Ck)  =  0.  Thus,  k^  +  k  and  J^Ck^)  J(k*),  but  the 
sequence  of  functions  Un(*)}  does  not  converge  to  anything  anywhere  except  at 
k  =  0. 


7.6  Conclusions 


The  results  developed  are  general  enough  to  apply  to  a  wide  variety  of 
fixed-fora  compensator  design  problems.  Current  efforts  are  aimed  at  special¬ 
izing  them  to  the  optimal  output  feedback  compensator  design  problem  that 
originally  motivated  these  efforts.  In  particular,  the  following  questions 
are  currently  under  active  investigation. 

(1)  How  can  the  conditions  of  Theorem  7-7  be  interpreted  physically? 

It  is  possible  that  approximation  results  for  strongly  continuous 
semigroups  can  be  used  to  guarantee  that  these  conditions  are 
satisfied  for  output  feedback  control  of  flexible  structures, 
subject  to  certain  conditions  on  the  inherent  damping  in  the 
structure. 

(2)  For  problems  in  which  the  convergence  conditions  are  met,  what 
factors  influence  the  rate  of  convergence?  This  question  is  of 
considerable  practical  importance  because  the  rate  of  convergence 
determines  how  many  modes  must  be  retained  in  a  finite  element 
model  of  a  structure  to  provide  a  reasonable  basis  for  control 
system  design. 
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SECTION  8 


LARGE-ANGLE  SPACECRAFT  SLEWING  MANEUVERS 


8.1  Introduction 

The  problem  of  large-angle  maneuvers  for  flexible  spacecraft  continues 
to  be  of  interest  in  the  aerospace  community.  This  section  develops  a  num-, 
ber  of  extensions  to  the  work  first  reported  in  Reference  8-1.  In  particu¬ 
lar,  this  section  presents  techniques  for  improving  the  optimal  torque  pro¬ 
files  by 

(1)  Allowing  the  solution  process  to  determine  the  optimal  terminal 
boundary  conditions  for  the  maneuvers. 

(2)  Developing  a  control-rate  penalty  technique  for  producing  smooth 
control  profiles. 

Several  example  maneuvers  are  provided  that  demonstrate  the  practical  appli¬ 
cation  and  utility  of  the  techniques  developed  herein. 

This  section  is  presented  in  five  parts.  Section  8.2  treats  the  dis¬ 
tributed  control  problem  and  develops  the  techniques  for  handling  slewing 
maneuvers  when  the  final  angle  for  the  maneuver  is  to  be  determined  as  part 
of  the  solution.  The  results  of  Section  8.2  are  appropriate  for  maneuvers 
where  the  final  angular  rate,  rather  than  the  final  orientation,  is  impor¬ 
tant. 


Section  8.3  treats  the  optimal  distributed  control  problem  and  devel¬ 
ops  a  control-rate  penalty  technique  for  smoothing  the  optimal  control 
torque  profiles.  The  control  designer  gains  the  ability  to  specify  the  ter¬ 
minal  control  state  and  control  rates  as  a  by-product  of  the  control-rate 
penalty  technique.  In  addition,  Section  8.3  combines  the  free  final  angle 
problem  of  Section  8.2  with  the  control-rate  penalty  technique. 

Section  8.4  develops  techniques  for  solving  a  single-axis  retargeting 
maneuver  for  a  rigid  spacecraft.  Computational  algorithms  are  proposed  for 
handling  systems  modeled  by  nonlinear  dynamics  and  time-varying  weighting 
matrices  in  the  performance  index. 

Section  8.5  presents  the  results  of  a  preliminary  study  of  slewing  ma¬ 
neuvers  for  the  ACOSS  Model  2.  All  example  maneuvers  are  presented  in  Sec¬ 
tion  8,6. 


8.2  Free  Final  Angle  Optimal  Slewing  Maneuvers 

The  specific  model  considered  (see  Figure  8-1)  consists  of  a  rigid  hub 
with  four  identical  elastic  appendages  attached  symmetrically  about  the  cen¬ 
tral  hub.  A  single-axis  maneuver  with  the  flexible  members  restricted  to 
displacements  in  the  plane  normal  to  the  axis  of  rotation  is  the  only  case 
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considered.  Furthermore,  the  body  (as  a  whole)  is  assumed  to  experience 
only  antisymmetric  deformation  modes  (see  Figure  8-2).  The  distributed  con¬ 
trol  system  for  the  vehicle  consists  of 

(1)  A  single  controller  in  the  rigid  part  of  the  structure. 

(2)  Four  appendage  controllers,  one  assumed  to  be  located  halfway 
along  the  span  of  each  of  the  four  elastic  appendages. 


Ul 


Figure  8-1.  Undeformed  structure. 


The  extension  to  the  case  of  multiple  controllers  along  each  appendage  is 
straightforward;  however,  in  Section  8.2  only  the  single  appendage- 
controller  case  is  presented. 

The  formulas  for  large-angle  maneuvers  presented  in  this  section  are 
developed  in  physical  space,  because  the  free  final  angle  transversaiity 
conditions  are  understood  most  easily  in  terms  of  physical  coordinates. 


8.2.1  Equations  of  Motion 

For  the  vehicle  shown  in  Figure  8-1,  the  equations  of  motion  can  be 
obtained  from  Hamilton’s  extended  principle  [8-2],  leading  to 


t  T  *  T  •»  •  •t  * 

(i  -  £  M  nje  +  MQnn  -  20n_  M  n 


U1  *  4u2 


(8-1) 


MflJ  +  +  lKnn  +  “  Fu9 

— on  nn—  1  nn  —  2 


(8-2) 
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Figure  8-2.  Antisymmetric  deformation. 


Referring  to  Eg.  (8-15)  in  Reference  8-3,  the  linear  time-invariant 
form  of  Eq.  (8-1)  and  (8-2)  foLlows  as 


MC  +  K£  = 


(8-3) 


M  =  (positive  definite) 

T 

K  =  K  (positive  semidefinite) 

Defining  the  state  variable  subsets  as 

£l  =  &  12  =  L  (8-4) 

leads  to  the  first-order  differential  equations 

"  I2  -2S  "M_1kIi  +  M_lpH.  (8-5) 

T  T  T 

Letting  j  >  the  state  space  equation  becomes 

£  =  Al  +  Bu  (8-b) 

where 
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8.2.3  Fixed  Time  Free  Final  Angle  Optimal  Control  Problem 


8. 2. 3.1  Statement  of  Problem 

The  rotational  dynamics  of  a  flexible  space  vehicle,  restricted  to  a 
single-axis  large-angle  maneuver,  are  considered  where  the  system  dynamics 
are  governed  by  Eq.  (8-6).  The  optimal  control  problem  is  to  find  the 
solution  to  Eq.  (8-6)  that  minimizes  the  performance  index 


-  4  f  fuTW  u  +  ZTW  Z]  dt 

7  *  nil—  —  ss— J 


(8-7) 


and  satisfies  the  specified  terminal  states,  where  W  is  the  control 
r  ’  uu 

weighting  matrix  and  W  is  the  state  weighting  matrix.  The  specified  ter- 

s  s 

minal  states  are  given  by 


[0O  JlT(to)]T  So  *  l®0  iT(to)l 


T-  mT 


(8-8) 


hf  -  n.(t£)  if  *  [e£  n.  (tf)] 


(8-9) 


where  the  constraint 


_n  (tf)  =  _n(tf)  =  0 


is  imposed  in  Eq.  (8-9).  In  addition,  the  final  angle  0f  is  not  specified 
since  it  is  to  be  determined  as  part  of  the  solution. 

8. 2. 3. 2  Derivation  of  the  Necessary  Conditions  from  Pontryagin's  Principle 

In  preparing  to  use  Pontryagin's  necessary  conditions,  the  Hamiltonian 
functional 


4(uTW  u  +  ZTW  Zj  +  AT(AZ  +  Bu) 
2V—  uu—  —  ss — '  —  —  — 


(8-10) 


is  introduced,  where  the  A^'s  are  Lagrange  multipliers  (also  known  as  costate 
or  adjoint  variables).  Pontryagin's  principle  requires,  as  a  necessary 
condition,  that  the  A^'s  satisfy  costate  differential  equations  derivable 
from 
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A 


(8-11) 


3H 

3s 


-WE-  ATA 

S»-  — 


and  that  the  conf  *ol  torque  u  must  be  chosen  at  every  instant  so  that  the 
Hamiltonian  of  Eq.  18-10)  is- minimized;  that  is,  for  continuous  u^  one 
requires  that 


=  0  ■  W  u  +  BTA  (8-12) 

0  u  —  uu—  — 


from  which  the  optimal  torque  is  determined  as  a  function  of  the  costate 
variables  as 


u  -  -w"1  BTA  (8-13) 

—  uu  — 

Substituting  Eq.  (8-13)  into  Eq.  (8-6)  yields  the  state  differential  equa¬ 
tion 


E  -  AE  -  BW_1BTA  (8-14) 

—  —  uu  — 


Because  the  final  angle  is  assumed  to  be  free,  the  transversality 
condition  providing  the  natural  boundary  condition  for  the  problem  follows 
as 


-AT(tfmf  -  0  (8-15) 

where  A_(t)  is  the  costate  vector  in  physical  space  and  SE^  is  the  variation 
of  the  physical  space  state  at  the  final  time.  Since  0£  is  free  and 

0j,  _n£,  and  _n£  are  specified,  6E£  in  Eq.  (8-15)  can  be  written  as 

6£f  =  [60  6n£  60 f  6uJ]T  *  [6ef0T0_0T]r  (8-16) 

Upon  introducing  Eq.  (8-16)  into  Eq.  (8-15),  the  desired  free  final 
angle  transversality  condition  follows  as 

Aj (i f)  60  -  0 
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or 


AxCtf)  «  0  (8--1?) 

since  68  ^  is  arbitrary,  where 

A  -  (Aj  A2  ...  A2n]T  (N  »  n  +  1) 

and  A^  litnotes  the  costate  variable  for  the  rigid-body  rotation  angle  0. 

Thus,  AjCtjJ,  6j,  rip  and  provide  the  2N  final  time  boundary  condi¬ 
tions  necessary  for  specifying  the  optimal  control  solution  of  Eq.  (8-11) 
and  (8-14),  when  fixed-time  free  final  angle  maneuvers  are  of  interest. 

8. 2, 3. 3  Solution  for  the  Initial  Costates 

In  Eq.  (8-11)  and  (8-14),  the  initial  boundary  conditions  for  are 
completely  known  while  for  A.  they  are  unknown,  and  the  final  boundary  condi¬ 
tions  for  and  A.  are  only  partially  known.  Thus,  the  application  of  Pon- 
tryagin's  principle  has  led,  as  usual,  to  a  two-point  boundary-value  problem 
(TPBVP).  To  solve  Eq.  (8-11)  and  (8-14),  the  merged  state  vector  is  defined 
as 

X  «  [ST  AT]T  (8-18) 

from  which  it  follows  that 

X  ■  nx  (8-19) 

where 


a  constant  coefficient  *  ^rix 
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or  equivalently 


Observing  that  the  known  quantities  m  Eq.  (8-22)  are 

A|(t^),  and  E^(0), . . .  and  carrying  out  the  partitioned  matrix  multi¬ 
plication  to  solve  tor  )» ♦  •  •  >57N^Cf t*ie  ecluation  defining  the 

solution  for  A(0)  follows  as 

U]A(0)  -  b  -  [BJ2(0)  (8-23) 

where 


Equation  (8-23)  is  in  the  standard  form  ior  the  linear  algebraic  equa¬ 
tion  Ax  ■  Jj,  which  can  be  solved  via  Gaussian  elimination  for  x  to  yield 
A(0). 

The  optimal  control  time  histories  are  obtained  upon  integrating 
Eq.  (8—1 1 )  and  (8—1 A)  subject  to  £(0)  given  by  Eq.  (8-8)  and  A(0)  given  by 
Eq.  (8-23).  Equation  (8-23.)  provides  the  initial  costates  in  physical 
coordinates.  However,  if  the  maneuver  simulations  are  modeled  in  modal 
coordinates,  the  following  two  transformations  permit  the  use  of  the  results 
of  this  section. 

First,  the  modal  space  state  transition  matrix,  $(t£,0)  is  mapped  to 
physical  space  via  the  transformation 
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*(tf,0) 


(8-24) 


The  derivation  for  the  transformation  matrix,  0,  can  be  found  in  either 
Reference  8-5  or  8-6. 

Second,  introducing  the  required  partitions  of  Eq.  (8-24)  into 
Eq.  (8-23)  and  solving  for  MO),  the  modal  space  initial  conditions  are 
obtained  from  the  transformation  equation 


where  MO)  is  the  modal  space  initial  coatate.  The  optimal  control  time 
histories  for  modal  space  are  obtained  upon  integrating  the  modal  space  form 
of  Eq.  (8-11)  and  (8-14)  subject  to  the  initial  conditions  provided  by 
Eq.  (8-25). 

Case  1  in  Section  8.6  provides  an  example  maneuver. 
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8 . 3  Control-Rate  Penalty  Technique  for  Optimal  Slewing  Maneuvers 

Introducing  control-rate  penalties  into  the  performance  index  accom¬ 
plishes  two  things.  First,  the  jump  discontinuities  in  the  control  time 
history  are  pushed  into  the  higher  control  time  derivatives.  In  other 
words,  the  control-rate  penalties  act  like  control  smoothing  penalties.  A« 
a  result,  the  frequency  content  of  the  resulting  control  profile  is  signif¬ 
icantly  reduced;  thus,  the  vehicle's  higher  frequency  modes  are  excited  only 
mildly  by  the  smoothed  control  profile. 

Second,  the  use  of  control-rate  penalties  permits  the  control  designer 
to  specify  the  terminal  control  state  and  control  time  derivatives.  For 
example,  in  many  reasonable  maneuvers,  it  is  convenient  to  specify  that  the 
slewing  control  system  is  turned  off  both  initially  and  finally.  Of  course, 
depending  on  particular  mission  objectives  (as  shown  in  the  retargeting  ma¬ 
neuvers  of  Section  8.4),  other  terminal  boundary  conditions  for  the  control 
system  are  possible.  To  fully  appreciate  the  significant  gain  in  perform¬ 
ance  chat  the  control-rate  penalty  technique  permits,  compare  the  example 
maneuvers  of  Section  8.6  with  the  results  presented  in  Section  8  of  Refer¬ 
ence  8-1. 

The  vehicle  depicted  in  Figure  8-1  is  assumed  to  be  modeled  for  the 
results  of  this  section.  Only  the  case  of  a  single-axis  maneuver  is  consid¬ 
ered,  where  the  elastic  displacements  occur  in  the  plane  normal  to  the  axis 
of  rotation,  and  only  antisymmetric  modes  are  modeled  (see  Figure  8-2). 

8.3.1  Problem  Formulation 


From  Eq.  (8-3),  the  equations  of  motion  can  be  written  as 


M£  +  K£  -  Pu 


(8-26) 


As  shown  in  Reference  8-1,  Eq.  (8-26)  can  be  cast  into  state  space  form, 
through  transformation  sequences  (t  t  -*■  s_j)  and  (s^jS^  +  s),  leading  to 


s  M  As  +  Bu  (8-27) 


The  slewing  maneuver  problem  is  to  find  the  solul  on  to  Eq.  (3-27)  that 
satisfies  the  terminal  boundary  conditions 
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The  optimal  control  problem  is  to  seek  the  torque  history  u(t)  that  will 
generate  an  optimal  solution  of  Eq.  (8-27),  initiating  at  Eq.  (8-28)  and 
terminating  at  Eq.  (8-29),  that  minimizes  the  performance  index 


(8-30) 


where  Wgg  is  a  weighting  matrix  for  the.  states,  is  a  weighting  matrix 

for  the  control,  and  W  ^  for  i  ■  l,...,k  are  weighting  matrices  for  the 

control  rates.  With  k  equal  to  zero,  the  performance  index  penalizes  only 
the  states  and  controls,  and  thus  the  control  problem  corresponds  to  a  prob¬ 
lem  of  minimum  control  energy,  kinetic  energy,  and  elastic  potential  ener¬ 
gy.  The  solution  for  this  problem  was  addressed  in  Reference  8-1.  With  k 
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greater  than  zero,  the  performance  index  penalizes  the  states,  controls,  and 
time  derivatives  of  the  controls.  The  example  maneuvers  in  Section  8.6  dem¬ 
onstrate  that  including  control-rate  penalities  leads  to  much  smoother  modal 
amplitude  and  control  torque  time  histories.  Additional  example  maneuvers 
can  be  found  in  Reference  8-6.  The  flexible-body  response  is  improved 
because  the  vehicle's  flexural  degrees  of  freedom  are  very  sensitive  to 
discontinuities  in  the  control  time  histories;  hence,  by  using  the  formu¬ 
lation  in  this  section,  the  jump  discontinuities  in  the  conti*ol  at  t  "  0  and 
t  *  tg  are  moved  to  the  higher  control  time  derivatives,  which  improves 

overall  system  performance. 

8.3.2  Necessary  Conditions  for  the  Optimal  Control  Problem 

Adjoining  the  equation  of  motion,  Eq.  (8-27),  as  a  differential  equa¬ 
tion  constraint  to  the  performance  index  in  Eq.  (8-30),  the  augmented  per¬ 
formance  index  is  given  by 


£  I  If  T„  .  .  r  (i)T„  (i) 


/  ItUXs!  +  1  £  1  wiiil  U]  +  i  [As  +  Bu  -  i]}  dt 


(8-31) 


where  the  Xjs  are  Lagrange  multipliers.  Using  standard  calculus  of 

A 

variations  techniques  the  tirst  variation  of  J  follows  as 


Lf  k  .  . 

6J  *  /  {ds^W  s  +  7  6u^^W..u^  +  6X^(As  +  Bu  -  s) 

0  "  SS“  i«0  "  ll~  “  -  “  ~ 


+  X^(A6s_  +  B6u_  -  5j0}  dt 


(8-32) 


♦  9  W  J 

Integrating  the  terms  containing  6^8  and  5ju  by  parts  for  i  *  l,...,k  leads 
to 


fcf  Cf  fcf 

/  T  •  T  t  #T 

/  r6s  dt  -  -X16s  +  /  Xi6s  dt 

0  0  0 


(8-33) 
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/  6£l)TW.  .u(l)  dt  *  l  (-l)j“16u(l"j)1IV..u(l+j"1) 

0  ll~  j-1  11-  0 


(-1)1  J  6uTW.  .u(2l)  dt 
0  11 


(8-34) 


since  £  and  £  (for  i  *  0,l,...,k-l)  are  specified  both  initially  and 
finally  in  Eq.  (8-28),  the  terminal  variational  terms  in  Eq.  (8-33)  and 
(8-34)  are 


6jB(tQ)  ■  6£(t^)  *  0 


^(tQ)  *  ^(tg)  *  £  (j  *  l,...,i) 


(8-35) 


where  t^  *  0.  Substituting  Eq.  (8-33),  (8-34),  and  (8-35)  into  Eq.  (8-32) 
yields 


f  k 


6J  *  /  {[  l  (-l)lW..u(2l)  +  BTXjT  6u  +  [W  £  +  ATX_  +  x]T  fig 


0  i«0 


+  [/£  +  B£  -  oX}  dt 


(8-36) 


Since  Pontryagin’s  principle  requires  the  first  variation  of  J  to  be 
zero,  the  necessary  conditions  follow  as 


*  As  +  Bu 


(8-37) 


X  ■  -W  s  -  A  X 

—  6S —  — 


(8-38) 


-  -5TA 

11—  — 


(8-39) 
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In  order  to  cast  Eq.  (8-39)  into  the  standard  space  form,  the  control 
state  is  defined  as 


n  -  f  T  T 

U  “  Ijij  u2 


T  il 
•  — 2kr 


(8-40) 


where 


u.  =  u(i  _  (i  8  1,2 . 2k) 

— l  — 


and  u.  is  an  N  *  1  vector  for  i  =  1,2,...  ,2k,  U  is  a  2kN  x  1  vector,  ai 
—l  c  »  »  >  »  _  c  > 

denotes  the  number  of  independent  controls.  Referring  to  Eq.  (8-40), 
Eq.  (8-39)  can  be  written  as 


U  =  CU  +  DX 


(8-41) 


where 


C21  C22 


°21  ’  (-1)k“\iW00 


(-i)kt2w-‘[o  :  «„  0  -«22  o  ...  ;.  o  ;.  (•i)\liW:»] 


(-I)k+1W"JBT 

kk 


where  0*  is  a  (2k  -  1)N  x  2N  null  matrix 

c 
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Defining  the  E  matrix  as  foLlows 


I" 

I 

I 


i 


E  *  IB  (oV) 


(8-42) 


and  substituting  Eq.  (8-40)  and  (8-42)  into  Eq.  (8-37),  the  state,  costate, 
and  control  differential  equations  follow  as 


State  equations 


s_  »  A  +  E  U_ 


(8-43) 


4 

% 

I 


Costate  equations 


X 


-w 


T 

s  -  AX 
si —  — 


(8-44) 


Control  equations 


U  «  CU  +  DX 


(8-45) 


8.3.3  Solution  for  the  Initial  Costatei  and  Control  Time  Derivatives 

In  Eq.  (8-43),  (8-44),  and  (8-45),  the  boundary  conditions  for  £  are 

known  both  initially  and  finally.  Those  for  £  are  split  so  that  half  are 

known  at  the  initial  time  and  half  are  known  at  the  final  time.  Those  for  £ 

are  totally  unknown.  Thus,  application  of  Pontryagin's  principle  has  led  to 
a  TPBVP.  To  obtain  the  solution  for  Eq.  (8-43),  (8-44),  and  (8-45),  the 
merged  state  vector  is  written  as 


X.  -  [sT  XT  HT]T  (8-46) 

from  which  it  follows  that 


X  -  AX  (8-47) 

Mi 


81 


where 


ft  «  -W  -A 

S3 


0  *  constant  coefficient  matrix 


Since  ft  is  constant,  it  is  well  known  that  Eq.  (8-43)  possesses  the  solution 


ft(t-tQ) 

X(t)  =  e  u  X(tQ) 


(8-48) 


where  is  the  (4N  +  2kN  )  x  (4N  +  2kN  )  exponential  matrix  which  is  ef- 

c  c 

ficiently  obtained  from  a  Psde  series  expansion.  Upon  setting 
ft(t-t0) 

e  *  4>(t,tg)  and  writing  out  Eq.  (8-48)  in  partitioned  form,  one  finds 


8(tf) 


*ss  *sX  *sua  *subl  /^(0) 


X(tf) 


*Xs  *xx  ^xua  *xuj  I2L(°) 


Sa(tf> 


4>  4>  .  41  $  Ju  (0) 

Yu  8  Yu  X  u  u  Yu  I -a 


4*  <t>  . 

V  V 


UKU» 

b  bi 


A<0) 


where 


(8-49) 


r  T  T 

hi 


l^k+1  --k+2 


m  .  fri 

**  -2k^ 
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where  is  known  both  initially  and  finally  and  u^  is  unknown  both  ini¬ 
tially  and  finally. 

Observing  in  Eq.  (8-49)  that  the  known  quantities  are  _s(0) ,  jKt^), 

u^(0),  and  u^(tg)  and  carrying  out  the  partitioned  matrix  multiplication, 

solving  for  s(t,)  and  u  (t,),  one  obtains 
—  f  — a  f 


(3-50) 

Setting  iMO)  *  u^t^)  *  £.  (i.e.,  turning  the  control  system  off  initially 

and  finally)  in  Eq.  (8-50)  (see  Eq.  (8-28)  and  (8-29;)  and  solving  Eq.  (8-50) 
for  \X0)  and  u^(0),  one  finds 


Equation  (8-51)  is  in  the  standard  form  for  the  linear  algebraic  equation 
Ajx  •  _b,  which  can  be  solved  via  Gaussian  elimination  for  x,  to  yield  X ( 0 ) 
an<^  H^®^*  optimal  control  time  histories  are  obtained  by  integrating 

Eq.  (8-43),  (8-44),  and  (8-45)  subject  to  the  given  state  and  control 
boundary  conditions  _s(C)  and  u  (0)  and  the  costate  and  control  rate  bound- 

ary  conditions  provided  by  Eq.  (8-51).  Cases  2  and  3  in  Section  8.6  present 
example  maneuvers  using  the  formulations  of  this  section. 


8*3*4  Free  Final  Angle  Maneuvers  with  Control-Rate  Penalties  Included 

Since  the  analytical  developments  for  solving  the  problem  of  free 

final  angle  maneuvers  are  similar  to  those  shown  in  Section  8.2,  only  kev 
equations  are  given.  3  y 
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As  shown  in  Section  8.2,  the  boundary  condition  for  the  free  final 
angle  corresponds  to  setting  the  final  value  of  the  costate  for  the 
rigid-body  angle  to  zero  (see  Eq.  (8-15)  and  (8-16)).  As  a  result,  the  set 
of  prescribed  boundary  conditions  follow  as 


io  -  [e„nX»T  io  ’ 


u 

-a. 


(8-52) 


Jlf  =  nttf) 


if  *  [efiT<tf)] 


u  =  u  (t.) 
— a,.  —a  f 


(8-53) 


and 


Aj  (t  f)  -=  0 


(8-54) 


where  we  impose  the  requirements  that  n_(tf)  *  n(tf)  *  0_  in  Eq.  (8-53)  and 
u^tg)  **  ^(tj )  a  £  in  Eq.  (8-52)  and  (8-53). 

Using  the  performance  index  of  Eq.  (8-31)  and  formulating  the  solution 
in  physical  spece,  the  state,  costate,  and  control  differential  equations 
follow  as: 

State  equation 


E  -  AE  +  EA 


(8-55) 


Costate  equation 


-WE-  A  A 

89—  — 


(0-56) 


Control .equation 


U  *  CU  +  DA 


(8-57) 
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where 


l  -  [c1  f]T 


Writing  the  merged  state  vector  as 


X  *  [e1  at  UT]T 


(8-58) 


the  solution  for  Eq.  (8-55),  (8-56),  and  (8-57)  can  be  written  as 


n(t-t  ) 

X(t)  *  e  °X(t0) 


(8-59) 


«<t-t  > 

where  9  is  defined  in  Eq.  (8-47).  Setting  e  “  $(t,tg)  and  performing 

th  partitioned  matrix  multiplication,  solving  for  £2(tg) . E2N^fcf^* 

A^(t^)  in  Eq.  (8-59),  the  linear  equation  to  be  solved  for  A_(0)  and-u^(O) 
follows  as 


U]a  =>  b  -  [B]c 


(8-60) 


where 


E2NA1 


.  .  .  ,  <{> 

E2A2N  E2Ub, 


.  .  .  .  0- 

2N  2N  2NUb, 


<f> 

«  •  •  v 


E2NUb 


A1A2N  Alub. 

l 


.  .  .  *, 


>  .  .  .  .  $  .  $ 
u  A.  u  A0  vi  u. 

a^  1  a^  2N  a^  b^ 


.  .  .  * 


u  u. 
al  bv 


u  A, 
*v  1 


...  $ 


«  U„  Uk 

av  2N  av  bj 


<f> 

*  •  •  * 


u  u, 
a  b 

v  v 
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[B] 


<f>  A 

E2E!  '  *  ‘  E2E2N 


E2NE1 

Vl 


u  E . 
si  1 


u  E , 
a  1 

v 


•  •  • 


£  E 
2N  2N 

K  r 

A1E2N 


\Z7* 


$ 


Ua  E2N 
\> 


£  *  [Aj CO)  ...  M„.(0)  M  (0)  Ub  ^  ***  % 

D1  v 


rj 

®  [Z2(t  f)  ...  E^Ct^)  A^(tf)  Ua  ^  f  ^  •••  ua  £>  ] 


c  =  [Z^O)  ...  ^2N-1(0)  E2N(0)  ua  (0)  •••'ua  (0)^ 


v  =  kNc 


Eq.  (8-60)  is  in  the  standard  form  for  the  linear  algebraic  equation 
Ax  *  b,  which  can  be  solved  via  Gaussian  elimination  for  x.,  yielding  MO) 
and  n^O) . 

The  optimal  control  time  histories  are  obtained  upon  integrating 
Eq.  (8-55),  (8-56),  and  (8-57)  subject  to  £(0),  u^(0)  given  by  Eq.  (8-52) 

and  MO),  u^(0)  given  by  Eq.  (8-60). 

Eq.  (8-60)  provides  the  initial  costate  and  control  rates  in  terms  of 
physical  space  coordinates.  However,  if  the  maneuver  simulations  are 
modeled  in  terms  of  modal  space  coordinates,  the  following  two  transforma¬ 
tions  permit  the  use  of  the  results  of  this  section.  First,  the  modal  space 
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state  transition  matrix,  i}>(t^,0),  is  mapped  to  physical  space  via  the  trans 
formation 

$(tf,0)  «  04>(tf,O)0-1  (8-61) 

where 


Second,  introducing  the  required  partitions  of  Eq.  (8-61)  into 
Eq.  (8-60)  and  solving  for  A(0)  and  u^(0),  the  modal  space  initial  condi¬ 
tions  are  obtained  from  the  transformation 


The  modal  space  optimal  control  time  histories  are  obtained  upon  inte¬ 
grating  the  modal  space  form  of  Eq.  (8-55),  (8-56),  and  (8-57'  subject  to 
the  initial  conditions  provided  by  Eq.  (8-62),  Case  4  in  Section  8.6  pro¬ 
vides  an  example  maneuver  using  the  formulations  of  this  section. 
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8.3.5  Selection  Technique  for  Weighting  Matrices  Wgs  and  W^.  in  the 
Control-Rate  Penalty  Formulation 

It  has  been  found  that  the  particular  choice  of  weighting  matrices 
significantly  affects  the  computational  effort  and  resulting  maneuvers.  A 
brief  study  of  the  effects  has  been  performed. 

In  the  example  maneuvers  of  Section  8.6,  the  following  block  diagonal 
forms  of  the  weighting  mr.tr ices  have  been  used 


W 

ss 


(2N  x  2N) 


W.  . 

11 


(N  x  N  ) 
c  c 


where  w  ,  w. .  are  scalar  quantities  and  p  is  defined  in  the  text  that 
ss  11 

follows. 

Since  the  problem  being  solved  is  a  slewing  maneuver,  the  rigid  body 

angular  displacement  is  not  penalized.  Hence,  the  first  element  of  W  is 

ss 

set  two  orders  of  magnitude  less  than  the  other  diagonal  elements.  The 

first  .element  of  W  is  not  set  to  zero  because  this  causes  repeated  .zero 
'  88  • 

eigenvalues  in  51  of  Eq.  (8-47),  which  can  potentially  lead  to  numerical 

diffuclties  in  calculating  the  state  transition  matrix. 

The  value  of  p  in  W„  determines  the  participation  of  each  set  of 

independent  controllers  in  the  slewing  maneuver.  Since  it  is  desirable  to 
have  the  rigid-body  torque  execute  most  of  the  slewing  maneuver  while  the 
appendage  torques  function  mainly  as  vibration  controllers,  p  is  set  to  2  or 
3  to  provide  a  smaller  pe’halty  o.i  the  rigid-body  torque. 

In  a  preliminary  parametric  study,  it  has  been  found  that  by  adjusting 

tne  norms  of  the  weighting  matrices  W  ,  W„-,  and  W. .  (for  i  *  l,...,k-l)  to 

ss  00  li 

several  orders  of  magnitude  smaller  than  the  norm  of.W  ,  the  eigenvalue 

bandwidth  of  51  in  Eq.  (8-47)  can  be  reduced  significantly;  thus,  the  sensi¬ 
tivity  of  the  optimal  control  problem  decreases.  However,  the  decrease  in 


88 


eigenvalue  bandwidth  of  fl  is  obtained  at  the  expense  of  increased  peak  flex¬ 
ural  deformations  during  the  maneuver,  though  the  time'  histories  of  the 
modal  amplitudes  are  considerably  smoother.  An  assejcwnt  of  the  impact 
this  weighting  matrix  scheme  has  on  the  unmodeled  defies  of  freedom  is  a 
topic  of  further  interest. 

8.4  Single  Axis  Retargeting  Maneuvers  for  a  Rigid  'Spacecraft 

The  problem  addressed  in  this  section  is  a  s' .t-le-axis  maneuver  of  a 
rigid  spacecraft  which  is  slewing  to  engage  a  mov-.u;,  target.  The  retarget¬ 
ing  maneuver  is  complicated  by  the  fact  that  the  ?. llowing,  parameters  and 
boundary  conditions  are  unknown: 

(1)  The  maneuver  time,  (t^  -  t^). 

(2)  The  final  engagement  angle,  0^. 

(3)  The  final  engagement  angular  rate, 

(4)  The  final  control  torque,  u^.. 

(5)  The  final  control  torque  tore,  u^. 


The  performance  index  selected  for  the  optimal  control  problem  repre¬ 
sents  a  tradeoff  between  elapsed  maneuver  time,  state,  and  control-rate 
penalties.  A  rigid  spacecraft  model  is  used  instead  of  a  flexible  vehicle 
to  simplify  the  calculations  and  provide  useful  insight  for  solving  the 
flexible  body  case.  In  the  developments-  Chat  follow,  the  target  trajectory 
is  assumed  to  be  known.  Since  work  on  tftis  topic  is  preliminary,  only  the 
solution  strategy  is  presented. 

The  solution  procedure  consists  of  the  following  three  steps.  First, 
the  linearized  free  final  time  problem  is  solved.  Second,  the  nonlinear 
solution  is  obtained  using  the  final  time  and  initial  costates  from  the 
first  step  as  initial  iteratives.  Third,  a  grid  search  is  performed  in  time 
to  find  the  final  time  that  minimizes  the  performance  index  for  the  non¬ 
linear  problem. 


The  target  trajectory  is  assumed  to  be 


8(t,tg)  =  tan 


-1  %  +  V(t  “  Cn)} 


(8-63) 


which  represents  a  linear  flyby  target  motion,  where  x^,  y^,  and  V  are 

constants  that  correspond  to  the  target's  initial  position  and  velocity  (see- 
Figure  8-3). 


t  CONSTANT  TARGET 
VELOCITY  AND 
DIRECTION 


s' 

s 

s 


S’ 

s 

INITIAL  TARGET  ANGLE 


INITIAL  TARGET 
POSITION  (x0,  yQ) 


-  X 


Figure  8-3.  Model  of  a  linear  flyby  target  motion. 


8.4.1  Solution  for  the  Linearized  Free  Final  Time  Problem 


The  solution  of  the  linearized  problem  provides  estimates  of  the  op¬ 
timal  final  time,  initial  costates,  and  control  time  derivatives  for  the 
nonlinear  solution. 

The  state  space  form  of  the  equation  of  motion  for  the  linear  problem 
is  given  by 


(8-64) 
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where  Sj  is  the  angular  displacement,  8.^  is  the  angular  velocity,  u  is  the 

control  torque,  and  I  is  the  moment  of  inertia  of  the  vehicle.  The  optimal 
control  problem  is  to  seek  the  solution  to  Eq.  (8-64)  that 

(1)  Minimizes  the  performance  index 

J  ■  /  (<*  +7woou2  +TWll“2  +TW22“21  dt  (8-65) 

where  a  is  a  weight  penalizing  elapsed  time,  and  Wqq,  Wjj, 

are  scalar  weights  on  the  control,  first,  and  second  control 
time  derivatives,  respectively. 

(2)  Satisifies  the  terminal  constraints 


"l^o1 

1 

0 

*2<t0> 

o 

CM 

QQ 

H 

(8-66) 

“V 

*  0 

“V 

-  0 

(8-67) 

Si  (tf ) 

-  5<tf,t0> 

(8-68) 

u(tf  ) 

•  “W 

u(t^) 

-  ie(tf,t0) 

(8-69) 

free 

(8-70) 

and  ett.tg)  is  the  angular  target  motion  given  by  Eq.  (8-63). 

The  final  boundary  conditions  given  by  Eq.  (8-69)  are  imposed  in  order 
to  establish  a  terminal  torque  and  torque  rate  that  will  allow  accurate 
pointing  and  tracking  of  the  target  at  the  end  of  the  maneuver.  The  target 

motion  in  terms  of  8(t,t0),  0(t,tQ),  §(t,tg),  and  e’U.tg)  is  assumed  to  be 

available  from  some  estimation  process,  for  practical  applications. 

Application  of  Pontryagin's  principle  leads  to  the  following  necessary 
conditions 
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S1  -  s2 


(8-71) 


s2  11  u/I 


(8-72) 


X1  "  0 


(8-73) 


X2  a  “X1 


(8-74) 


W00u  -  WnU  +  W22u  =  — 


(8-75) 


On  imposing  the  terminal  constraints  listed  in  Eq.  (8-68)  and  (8-69),  the 
transversality  conditions  governing  the  optimal  maneuver  follow  as 


G1  5  -X1(tf)6(tf ,t0)  -  X2(tf)8(tf ,tQ)  +  X1(tf)s2(tf)  +  X2(tf)u(tf)/1  +  a 


+  2  W00u  +  2  Wll^  +  2  W22U  *  0 


(8-76) 


G2  =  "  6(tf,t0)  =  0 


(8-77) 


s2(tf)  -  9(tf,t0) 


(8-78) 


=  u(t£ )  -  ie(tf,t0) 


(8-79) 


G5  =  u(tf)  -  I9’(tf,t0)  «  0 


(8-80) 


To  cast  Eq.  (8-75)  into  state  space  form,  the  control  state  is  defined 


lul  u2  u3  u4)  "  «  «  «] 


(8-81) 


Subject  to  Eq.  (8-81),  Eq.  (8-71)  to  (8-75)  can  be  written  in  the  merged 
state  vector  form 


where 


x  -  [« 


x  =  Ax 


•#  J 


(8-82) 


and 


0  10  0 


0  0  O'  0 


A 


0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

Q 


0 

0 

-1 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


1 


i 

I 

0 

0 

0 

0 

0 

woo 


0 

0 

0 

1 

0 

0 

0 


0  0 


0 

0 

0 

1 

0 


0 

0 

0 

0 

1, 

0 


Since  A  is  constant,  it  is  well  known  that  Eq.  (8-82)  possesses  the 
solution 


Alt- ty ) 

x(t)  =  e  x(tQ) 


(8-83) 


where 


A(tf"t0)  ‘  ' 

e  is  the  exponential  matrix 


Setting 


A(tf-t0) 

e 


93 


introducing  the  appropriate  partitions  of  $(tf,t0)  into  Eq.  (8-50)  and  re¬ 
arranging,  the  linear  system  defining  the  solution  for  X^(tQ>, 
u^Uq),  and  follows  as 

!(tf)  -  4»ss8(0) 

— a(tf }  “  *u  s^(0) 
a 

where 

S  -  [S1  82]T 

K  -  [\  *2]T 

"a  "  h  u2^t 

%  "  fu3  UJT 


The  optimal  control  time  histories  follow  on  integrating  Eq.  (8-82)  subject 
to  8(0)  and  u  (0)  given  by  Eq.  (8-66)  and  (8-67)  and  X(0)  and  u,  (0)  given 

by  Eq.  (8-84). 


Due  to  the  nonlinear  character  of  the  transversality  conditions,  the 
necessary  equations  are  solved  by  iteration.  Starting  iteratives  for  the 
transversality  conditions  are  obtained  by  assuming  that  the  final  time  is 
known.  As  a  result,  a  number  of  final  time  grid  points  are  selected.  At 
each  final  time  grid  point,  the  optimal  control  solution  is  obtained 
subject  to  satisfying  the  transversality  conditions  of  Eq.  (8-77)  through 
(8-80),  and  the  performance  index  is  computed.  Once  a  minimum  is  found  for 
the  performance  index  within  three  consecutive  final  time  grid  points,  the 

k  k  k  2 

performance  index  is  modeled  by  J^t^)  *  Cq  +  Cjt^  +  ^(tf)  w^ere  *7^  de- 

notes  the  k*"*1  performance  indtx  and  t ^  denotes  the  kC^  final  time  grid 
point.  Upon  solving  for  Cq,  Cj ,  and  C^,  the  starting  iterative  for  t^  is 
obtained  by  setting  the  derivative  of  J^tt)  equal  to  zero,  yielding 


t 


f 


-V2v 


Integrating  Eq<  (8-82)  and  substituting  the  final  states,  costates, 
and  control  states  into  Eq.  (8-76)  to  (8-80)  provides  the  errors  in 
satisfying  the  necessary  conditions.  The  corrections  to  the  initial 
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costates,  control  time  derivatives,  and  final  time  are  obtained  from  the 
following  equation 


and  the  desired  corrections  are  denoted  by 

ACx  for  i  =  1,...,5 


Using  the  updated  final  time,  the  terminal  errors  for  Eq.  (8-76)  through 
(8-80)  are  computed  on  integrating  Eq.  (8-82),  and  if  necessary,  the  process 
is  repeated  until  the  final  time  converges. 

8.4.2  Retargeting  Maneuvers  with  Nonlinear  System  Dynamics 

A  retargeting  maneuver  is  considered  where  the  system  dynamics  are 
nonlinear  (this  example  is  similar  to  the  nonlinear  flexible  problem  shown 
in  Reference  8-1)  and  time-varying  weights  are  included  in  the  performance 
index.  Specifically,  the  equation  <~f  motion  is  assumed  to  be  given  by 


(8-86) 


and  the  new  performance  index  follows  as 


'tan 
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The  terminal  boundary  conditions  listed  in  Eq.  (8-66)  through  (8-69) 
are  assumed  for  this  problem,  and  the  initial  estimates  for  A(0)  and  u^(0) 

are  given  by  Eq.  (8-84)..  Written  in  merged  state  vector  form,  Eq.  (8-88) 
through  (8-92)  become 

a  -  Ax  +  a.  b  (8-93) 

—  —  k— 

where 

0 

-<82/1 

-y(t)(s1  -  0) 

-B(t)(s2  -  6)  -  2A2(y)s2 

0 
0 
0 
0 


and  o^  is  the  k  relaxation  parameter,  where  the  sequence  of  parameters 

{oq  =  0  <  Oj  <  ...  <  cr_j  <  o^  =  l}  has  been  introduced  and  r  is  preset. 

The  solutr-on  for  Eq.  (8-93)  is  refined  iteratively  using  the  following  equa¬ 
tion 


-fD  "  if! 


(8-94) 


where 

x^D  is  the  desired  final  state  vector 
is  the  integrated  final  state  vector 

T  T 

[9Xf/3Xo]  is  the  state  transition  matrix 
Ax»  is  the  correction  vector 
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T  T 

Assuming  that  Ux^Xq]  in  Eq, 


(8-94)  can  be  expressed  in  the  form 


*(tf,t0)Q(tf,t0) 


(8-95) 


where  ^Ct^t^)  is  the  state  transition  matrix  for  the  linear  time-invariant 

part  of  the  problem,  and  Q(tf,tQ)  is  the  state  transition  matrix  for  the 

nonlinear  part  of  the  problem,  it  can  be  shown  using  Lagrange’s  variation 
of  parameters  method  that  Q(tpt^)  satisfies  the  equation 


Q(tf,t0)  =  1  +  0,,/ 


4>-1  (t  »tQ) 


W—  m  — 1  T 

3bT(t) 


3x 


4> ( T  ,tQ)Q(T  >tQ)  dT 


(8-96) 


where  the  relaxation  parameter,  o^,  has  been  introduced  artificially.  Ob¬ 
serving  that  Eq.  (8-96)  is  an  integral  equation  for  Q(t,tQ),  and  using  the 
method  of  successive  substitutions,  one  finds  the  following  equation 


*• 

Q(tf »tQ)  *  i +  \  !  ♦  (T>to) 

co 


3bT(x 

3.x 


-hT 


4) (t , tQ )  dt 


(8-97) 


where  the  higher  order  terras  in  have  been  assumed  to  be  small. 

The  solution  procedure  is  summarized  as  follows.  First,  the  relaxa¬ 
tion  parameters  are  chosen  such  that  {o^  =  0  <  . . .  <  =  1 } . 

Second,  the  initial  costates  and  control  time  derivatives  for  the  linear 
solution  are  used  as  initial  conditions  for  the  integration  of  Eq.  (8-93) 
with  k  *  1  and  F.q.  (8-97).  The  final  states  and  controls  are  compared  with 
the  desired  values  at  the  final  time,  and  if  the  terminal  errors  exceed  a 
prescribed  tolerance,  then  corrections  to  the  initial  costates  and  control 
time  derivatives  are  calculated  from  Eq.(8-94),  (8-95),  and  (3-97).  The 
integrations  and  corrections  are  repeated  until  the  norm  of  the  terminal 
error  vector  is  sufficiently  small.  The  index  on  the  relaxation  parameter 
is  increased,  and  the  iteration  process  continues  until  the  solution  for 
J  is  found.  The  solution  for  o  **  1  yields  the  desired  nonlinear 
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solution.  It  is  anticipated  that  Eq.  (8-97)  will  lead  to  significant  compu¬ 
tational  saving3  when  [3bJ(f)/8jt]  is  sparse. 

The  optimal  final  time  for  the  nonlinear  solution  is  obtained  by  a 
grid  search  for  the  final  time  that  minimizes  the  performance  index  of 
Eq.  (8-87).  That  is,  a  series  of  nonlinear  solutions  and  performance 
indices  are  calculated  for  different  final  times.  A  quadratic  fit  is  then 
performed  to  estimate  the  final  time  that  minimizes  the  performance  index. 

8 . 5  Preliminary  Study  of  Slewing  Maneuvers  for  the  ACOSS  Model  2 

Several  single-axis  slewing  maneuvers  have  been  computed  for  a  modi¬ 
fied  version  of  the  ACOSS  Model  2.  The  modifications,  which  simplify  the 
calculations,  include: 

(1)  Locked  isolator  springs. 

(2)  A  symmetric  mass  distribution  for  the  vehicle  about  its  center 
of  mass. 

The  "rigid-body"  torque  is  applied  in  the  z  direction  at  the  center  of  the 
equipment  section  (node  44).  The  formulation  of  Section  8.3  is  used  with 
k  *  2  in  Eq.  (8-30)  for  the  performance  index.  The  mass,  stiffness,  and 
actuation  matrices  are  computed  by  the  Dynamic  Interaction  Simulation  of 
Controls  and  Structures  program  (DISCOS)  which  uses  NASTRAN  finite-element 
data. 


The  solution  procedure  is  as  follows.  A  number  of  modes  to  be  con¬ 
trolled  is  chosen.  The  optimal  control  time  history  is  computed,  using  the 
formulation  of  Section  8.3  and  the  mass,  stiffness,  and  actuation  matrices 
from  DISCOS.  The  control  time  history  is  then  input  to  DISCOS  which  simu¬ 
lates  the  maneuver  and  treats  all  nonlinear  kinematic  effects. 

Results  from  the  DISCOS  simulations  show  that  the  uncontrolled  sym¬ 
metric  modes  have  small  vibrations.  Though  the  controlled  antisymmetric 
modes  have  larger  peak  amplitudes,  the  resuits  show  that  the  modal  ampli¬ 
tudes  are  zero  at  and  beyond  the  final  time  of  the  maneuvers,  which  is  the 
expected  response.  An  example  maneuver  is  shown  in  Section  8.6,  Case  5. 

8.6  Example  Maneuvers 

This  section  presents  example  maneuvers  for  the  formulations  discussed 
in  Sections  8.2  through  8.5.  For  Cases  1  to  4,  the  geometry  of  Figure  8-1 
is  assumed  with  the  following  configuration  parameters:  the  moment  of 
•  •  2 

inertia  of  the  undeformed  structure,  I,  is  7000  kg-m  ;  the  mass  per  unit 
length  of  the  four  identical  elastic  appendages,  p,  is  0.004  kg/m;  the 
length  of  each  cantilevered  appendage,  L,  is  150  m;  the  flexural  rigidity  of 

3  2 

the  appendages,  El,  is  1500  kg-m  /s  ;  and  the  radius  of  the  rigid  hub,  r,  is 
1  m.  In  the  integrations  over  the  mass  stiffness  distributions,  the  radius 
of  the  hub  is  not  neglected  in  comparison  to  the  appendage  length.  The  fol¬ 
lowing  comparison  functions  are  used  as  assumed  modes 

99 


* 


4>  (x  —  r) 
P 


1  -  cos 


pir(x  -  r) 
L 


(.P?(.xL.-.r.>)2 


P 


1,2,., 


which  satisfy  the  geometric  and  physical  boundary  conditions 
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x=r 


♦••I  - 

P!x3r+L  p  x=r+l 
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of  a  clamped-free  appendage.  For  Case  5,  the  ACOSS  Model  2  is  used  where 
the  isolation  springs  are  locked  and  the  mass  distribution  has  been  modi¬ 
fied. 


In  Cases  1  to  4,  all  the  W  matrices  are  set  to  diagonal  matrices  in 

S  8 

physical  space,  and  then  mapped  into  modal  space  using 


[w 

L  q  e 


s  s ^ modal 


— 

— 

et 

0 

[w  ]  . 

1  88 J physical 

E  0 

0 

EI 

0  E 

For  Case  5,  the  W  matrix  is  set  to  a  diagonal  matrix  in  modal  space.  Re- 
s  s 

ferring  to  Table  8-1  and  Figures  8-4  to  8-14,  the  example  maneuvers  are 
described  as  follows. 

8.6.1  Case  1  (Figure  8-4) 

Case  1  is  a  60-second  free  final  angle  rotation  reversal  maneuver 
using  the  formulation  of  Section  8.2.  The  structure  returns  to  its  original 
angular  position  at  the  end  of  the  maneuver.  The  modal  amplitudes  are  nega¬ 
tive  and  the  control  torques  are  positive  throughout  the  maneuver.  In  addi¬ 
tion,  the  control  torques  have  jump  discontinuities  at  the  initial  and  final 
time. 


8.6.2  Cases  2  and  3  (Figures  8-5  through  8-7) 

Cases  2  and  3  are  60-second  rest-to-rest  maneuvers  using  the  formula¬ 
tion  of  Section  8.3.  The  states,  controls,  and  control  rates  are  included 
in  the  performance  inden  for  Case  2,  while  the  second  and  third  time  deriva¬ 
tives  of  the  controls  are  included  additionally  for  Case  3.  For  both  cases, 
the  control  torques  at  the  initial  and  final  times  are  all  zero;  therefore, 
they  are  continuous  at  those  points.  The  slope  of  the  maneuver  angle  versus 
time  plots  are  very  small  at  the  initial  and  final  times  which  is  consistent 
with  a  smooth  modal  amplitude  history.  Case  3  differs  from  Case  2  in  that 
there  is  a  30-to-l  reduction  in  the  first  mode  peak  amplitude,  a  2-to-l  re¬ 
duction  in  the  second  mode  peak  amplitude,  and  a  20-to-l  reduction  in  the 
peak  tip  deflection.  On  the  other  hand,  Case  3  shows  a  l-to-1,3  increase  in 
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Table  8-1.  Parameters  used  in  test  case  maneuvers. 


- - 

Case 

No. 

No.  of 
Modes 
Controlled 

No.  of 
Controls 

1 

1  ° 

!  CD  U 

9o 

(rad) 

ef 

(rad) 

8f 

(rad) 

W 

ss 

W00 

W11 

W22 

W33 

1 

2 

5 

0 

-0.5 

free 

0.5 

10 

i 

— 

— 

— 

2 

2 

5 

0 

0 

IT 

0 

-3' 
10  J1 

10"9I 

i 

— 

— 

3 

2 

5 

0 

0 

TI 

0 

10“3I 

10-9I 

-9 

10  ’i 

-9 

10  *1 

.  i 

4 

2 

5 

0 

-0.5 

free 

0.5 

i<r3i 

10~9I 

io"9i 

i 

— 

5 

2 

1 

0 

0 

tt/36 

0 

IQ'7! 

10*Ai 

10"A1 

I 

— 

A  -2 
I  is  an  identity  matrix  with  the  first  element  set  to  10  ,  setting  a  lower 

weight  on  the  maneuver  angle. 

-3 

I  is  an  identity  matrix  with  the  first  element  set  to  10  ,  setting  a  lower 

weight  on  the  rigid-body  control  or  control  time  derivative. 


the  rigid-body  peak  torque  requirement  and  a  l-to-8.4  increase  in  the 
elastic  appendage  peak  torque  requirement  when  compared  with  Case  2.  The 
slopes  of  the  control  time  history  plots  for  Case  3  are  zero  at  the  initial 
and  final  times  because  of  the  higher  order  time  derivatives  in  the  perform¬ 
ance  index. 

8.6.3  Case  4  (Figures  8-8  and  8-9) 

Case  4  is  a  60-second,  free  final  angle,  rotation  reversal  maneuver 
using  the  formulation  of  Section  8.3.4  with  the  performance  index  penalizing 
the  states,  controls,  and  the  first  and  second  time  derivatives  of  the  con¬ 
trols.  Compared  with  Case  1,  Case  4  shows  a  19-to-l  reduction  in  the  first 
mode  peak  amplitude  and  a  2-to-l  reduction  in  the  second  mode  peak  ampli¬ 
tude.  The  rigid-body  peak  torque,  on  the  other  hand,  increases  by  a  1-to- 

1.3  ratio  and  the  appendage  peak  torque  by  a  l-to-13  ratio.  The  initial  and 
final  control  torques  are  zero,  and  therefore  continuous,  while  in  Case  1 
they  are  discontinuous  at  the  initial  and  final  times.  The  modal  amplitude 
and  control  torque  histories  are  much  smoother  in  Case  4  than  in  Case  1. 

From  these  results,  it  can  be  concluded  that  the  inclusion  of  ccntrol-rate 
penalties  in  the  performance  index  produces  torque  and  modal  amplitude  time 
histories  which  are  much  smoother  than  those  produced  by  the  formulation 
without  the  control-rate  penalty. 
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8.6.4  Case  5  (Figures  8-10  through  8-14) 


Case  5  is  a  5-second,  5-degree,  rest-to-rest  maneuver  of  the  ACOSS 
Model  2  about  the  z  axis.  The  maneuver  is  simulated  for  7  seconds  in  order 
to  observe  the  residual  vibrations.  The  formulation  of  Section  8.3  is  used 
with  the  performance  index  penalizing  the  state,  control,  and  first  and 
second  time  derivatives  of  the  control.  The  first  two  antisymmetric  modes, 
Modes  10  and  11,  are  controlled  with  the  others  uncontrolled.  The  results 
show  that  the  symmetric  modes  are  only  slightly  excited.  Mode  7,  which  is  a 
vibration  in  the  x-z  plane,  and  Mode  8,  which  is  a  symmetric  vibration  in 
the  x-y  plane,  are  excited  only  slightly.  Modes  9,  12,  and  18,  which  are 
mainly  vibrations  in  the  x-z  plane  with  some  antisymmetric  components  on  the 
x-y  plane,  are  excited  somewhat  more  than  Modes  7  and  8  because  of  the 
antisymmetri  ,  component. 

Modes  12  and  18  have  small  amplitude  residual  ringing  at  the  end  of 
the  maneuver,  but  Mode  9  has  a  relatively  large  amplitude  at  the  terminal 
time.  Modes  10  and  11,  which  are  the  controlled  antisymmetric  modes,  are 
excited  but  are  forced  to  have  essentially  zero  amplitude  at  the  end  of  the 
maneuver,  as  expected.  Mode  15,  which  is  a  twisting  of  the  solar  panels,  is 
hardly  excited.  Mode  17,  which  is  purely  antisymmetric  in  the  x-y  plane,  is 
excited,  but  the  amplitude  is  greatly  reduced  at  the  terminal  time  with 
small  amplitude  residual  ringing.  The  large  residual  amplitude  of  Mode  9 
suggests  that  it  should  be  controlled  in  future  simulations. 
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Figure  8-11.  Case  5,  time  response  of  Modes  7,  8,  and  9. 


AMPUTUOE  (ml 


TIMEUi 


- - - - - - . - - - 1 - 1 - 1 - 1 - 1 - - - - - - - - - A- 

2. HO  3.20  H.OD  U. 80  5.60  6. HO  7.20  6. 00 

TIMEM 


TIME(t) 


Case  5,  time  response  of  Modes  15,  17,  and  18. 
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Figure  8-14.  Case  5,  body  £ramo>  maneuver  angle 


SECTION  9 


ORDER  REDUCTION  BY  IDENTIFICATION — 
SOME  ANALYTICAL  RESULTS 


9.1  Introduction 


Considerable  attention  in  the  community  studying  large  space  structures 
(LSS)  has  been  de/oted  to  the  design  of  (stable)  controllers.  The  major  dif¬ 
ficulties  encountered  are  the  dimensionality  problem  and  the  spillover  prob¬ 
lem.  Various  devices  have  been  suggested  to  overcome  these  difficulties. 
Roughly  speaking,  the  common  philosophy  of  all  the  designs  is:  select  n  domi¬ 
nant  modes  (n  relatively  small),  design  a  controller  for  the  n  modes  while 
preventing  (or  minimizing)  the  spillover  effect.  Although  this  is  a  plausible 
philosophy,  it  suffers  from  a  major  flaw:  the  success  of  the  corresponding 
design  depends  on  knowledge  of  a  precise  high-order  model.  Thus,  the  burden 
of  dealing  with  a  high-order  model  is  transformed  from  a  control  design  prob¬ 
lem  to  a  modeling  problem.  Unfortunately,  the  identification  of  high-order 
models  is  even  more  imperilled  with  numerical  difficulties  than  controlling 
such  systems.  Thus,  it  seems  that  a  compromise  between  these  difficulties  is 
imminent.  Actually,  such  a  compromise  is  practiced  regularly  by  engineers  as 
follows: 

(1)  A  model  structure  is  selected  (e.g,,  linear  of  order  n) . 

(2)  Using  experimental  data,  the  best  model  of  the  selected  structure 
is  estimated. 


(3)  A  controller  is  designed  for  the  "best  model". 

(4)  The  performance  of  the  controller  when  applied  to  the  plant  is 
tested. 


Although  a  common  practice,  no  proof  exists  that  such  a  scheme  will  work,  ex¬ 
cept  for  Step  (4).  The  purpose  of  this  section  is  to  evaluate  the  feasibility 
of  this  approach  theoretically.  Specifically,  an  attempt  is  made  to  charac¬ 
terize  control  designs  that  will  guarantee  stability  despite  the  order  reduc¬ 
tion  induced  by  the  identification.  Since  the  least  squares  (LS)  method  is 
one  of  the  more  robust  identification  schemes,  and  since  analytical  expres¬ 
sions  for  order  reduction  exist  for  the  LS,  this  method  is  analyzed  herein. 

Section  9.2  discusses  general  properties  of  LS  identification,  espe¬ 
cially  those  pertaining  to  order  reduction.  A  general  stability  theorem  for 
the  discrete-time  system  and  its  relation  to  LS  is  presented  subsequently. 

Sections  9.3  and  9.4  characterize  controllers  which  are  stable  when  de¬ 
signed  based  on  reduced-order  models.  It  is  shown  that  stability  enhancing 
controllers  (i.e.,  controllers  which  increase  damping)  are  robust  when  applied 
to  systems  with  reduced-order  models.  In  the  last  section  conclusions  are 
drawn  and  extensions  are  outlined. 
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Nomenclature 


All  vectors  are  denoted  by  underlined  lower  case  letters,  e.g.,  ,  _a^, 

n 

a_  . 

The  j  entry  of  x  is  x^ . 

Matrices  and  operators  are  denoted  by  capital  letters,  e.g.,  P  ,  X,  E,  J 

The  dimensions  of  the  matrices  (operators)  are  not  specified  if  they  can 
be  deduced  from  the  context. 

9.2  Least-Squares  Identification 

To  facilitate  the  analysis  of  the  subsequent  section,  we  present  first 
some  basic  results  in  LS  theory. 

Let  {xq.Xj.--.}  be  a  set  of  vectors  in  a  Hilbert  space  H. 

Let 


L 

n 


A 


Span  { x^ , . . . ,x^} 


and 


*  n  A  _  A  .  .  ,  . 

=  P  x.  ''  projection  of  x_  on  L 
-0  n-0  -On 


(9-1) 


r  n 

=  )  a.x. 

.  .  L— 1 

1  =  1 


Then  the  following  results: 
Result  9-1 


n+1 


l 

i=l 


(9-2) 
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(9-3) 


An+1  _  An  An+1  „  .  _  , 

3  ■  cL  •  "  3  i  Cl  •  *  i  X  1  .  •  •  •  .U 

i  i  n+1  n+l,i 


where 


PiAi+l 


11 

7  a  . ,  .  x. 
. L.  n+1,  l—i 
i=l 


and  where  an  .  x  .is  Che  projection  of  x.  on  x  x  .  *  x  ..  -  P  x  ... 

n+1  -n+1  r  J  -0  —n+1  -n+1  —n+1  n-n+1 

Result  9-1  is  best  illustrated  as  shown  in  Figure  9-1. 


/.n+1  - 
an+i  5n+1 


Pn*n+1 


Figure  9-1. 


The  algebraic  relations  which  correspond  to  Result  9-1  can  be  expressed 
as  follows 


P  x_  -  X  (XTX  )“1XTxn 
n-0  n  n  n  n-0 


(9-4) 


116 


where 


X  is  the  matrix 
n 


X 

n 


•in 


] 


and  the  superscript  T  denotes  the  transpose  operation.  Defining  the  vector 


and  the  truncation  operator  E 


Ex 


we  have 


\-l  -T 

x  . .  J  x  .  ,  x„ 


-o+l 


-n+1-0 


(9-5) 


E<;n+1 


"n+1 

Ai+1 


(xTx  ) 

v  n  n  ’ 


-1 


T 

X  x  . 
n-n  +  1 


(9-6) 


The  application  of  Eq.  (9-5)  and  (9-6)  resuLts  in  the  system  identi¬ 
fication  theory  as  follows.  Let 


n  m 

y(k)  -  \  a.y(k  -  i)  +  \  bu(k  -  1)  +  w(k)  (9-7) 

i=l  i=0  1 


be  an  (n,m)  AR1MA  model  to  be  identified,  where  (y()},  (u()},  and  (w()}  are 
output,  input,  and  noise  (residuals)  sequences,  respectively.  Define 
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5i 


2i 


■y(n  -  i) 
y(N*-  i) 


i  *  0) » •  *  )H 


['  u(2n  +  1  -  i)  1 
u(N  +  n  +  1  -  i) J 


i  »n  +  l,.,.,n  +  n 


where  N  is  the  length  of  the  data. 

The  LS  estimate  of  the  (n,m)  ARIMA  model  can  be  transformed  to  the 
(n  +  l,m)  or  the  (n,m  +  1)  order  ARIMA  model  using  the  relations  above. 

Assuming  that  the  number  of  unstable  poles  is  known  and  choosing  n 
larger  than  that  number,  it  can  be  shown  [9-1]  that  the  unstable  modes  can  be 
identified  exactly  using  the  LS  algorithm.  Thus,  in  the  sequel,  only  the  ef¬ 
fect  of  order  reduction  on  stationary  (stable)  systems  is  considered. 

In  the  stationary  case,  assume  that  N  is  large  enough  such  that 


i  l  y(k)y(k  -  Z)  *  E[y(k)y(k  -  £)]  -  v(£)  (9-8) 

k-1 


Thus,  the  LS  identification  of  the  (n,0)  model  is  given  by  the  last  column  of 

V  where 
n 

Vn  =■  Toep  [v(0)  ,v(l) , . . .  ,v(n)  ] 


A 


v(0)  v(l)  ...  v(n) 

\ 

v(l)  v(0)  \ 

X  X 
'  \  X 

Lv(n)  \  >. 


(9-9) 


It  follows  from  the  Toeplitz 
can  be  written  as 


structure  of  V  that  the  relation  in  Result 
n 


9-1 


118 


An+1 


An  An+1  ,'n 

a  -  a  ,  Ja 
—  n+i  — 


An+1 

an+2 


(9-10) 


where 


J  =A 


0  1 

r  **o 


i.e. 


_xr 

"x 

n 

= 

• 

• 

X 

n 

Equation  (9-10)  is  the  basis  fo..  the  celebrated  lattice  structure  of  stable 
systems  (see,  e.g.,  [9-2]).  Note  that 


1  -  ( 


n+1  \2 
an+l  J 


(Ea[ 


n+1  An+1  An+li 
+  a  JEa  ] 
n+1  —  1 


(9-11) 


The  intimate  relation  between  Vn,  a_  ,  and  stability  theory  has  been 

discussed  in  a  number  of  articles  [9-3,  9-4,  9-3]  and  is  summarized  in  Theo¬ 
rem  9-1 . 

Theorem  9-1 


Statements  (1)  through  (5)  regarding  the  polynomial 


P  (z)  =  zn  -  £  a.z 

n  i-1  1 


n-1 


n 

b,i 
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and  the  corresponding  companion  matrix 


A  = 


“o 

“ 

• 

• 

• 

1 

0__ 

a 

•  •  •  v  « 

n 

1 

are  equivalent. 


(1)  P  (z)  =  det[zl  -  A]  has  its  roots  inside  the  unit  circle. 

n 

(2)  The  Toeplitz  matrix  solution  o£  the  Ricatti  equation 


AV  AT  + 
n 


[0  ...  0  1] 


n 


is  positive  definite, 

(3)  Mj,,  defined  in  Eq.  (9-13)  is  positive  definite 


M 


1\  0 
-a. 

-a  '  •  •  “a,  ^  1 

n  1 


’V"  "an 

0 

\! 


0\  0 
-a  \ 

VX 

L"al  _an 


0  \  v  " 
\  n 


A1AI  -  A2A? 


(9-12) 


(9-13) 
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(4) 


The  Schur-Cohn  matrix  SC^,  defined  beLow,  is  positive  definite 


A 1  A  2  A  n 

(5)  Let  _a  ,  a  ,...,a  be  the  LS  parameter  vectors  of  the  first  to 

the  n  order  AR  models  obtained  from  the  stationary  data  (y(i)}. 

a  ^ 

then  ja^  j  <  1 ,  i  *  1 , . . . ,n. 


(6) 


Let  a0,  a°  \.,.,a*  be  a  set  of  vectors  obtained  using  Eq.  (9-10) 
and  (9-11)  repeatedly.  Then,  |a^|  <  1,  i  =  l,...,n. 


Theorem  9-1  shall  serve  in  assessing  the  effect  of  order  reduction  on 
the  stability  of  a  feedback  design. 


Note  in  particular  that  if  P  (z)  is  stable,  V  of  Statement  2  and  M  *  of 

n  n  n 

Statement  3  are  identical  to  the  covariance  matrix  of  the  (n,0)  AR1MA  model 

excited  by  unit  variance  white  noise.  The  main  benefits  in  the  stability 

statements  of  the  theorem  above  are  that  they  tie  stability  to  model  reduction 

and  LS  estimation  and  that  explicit  expressions  (the  matrices  M  or  V  )  to  be 

r  n  n 

used  in  stability  analysis  are  given  in  terms  of  the  mcdel  parameters  (as  op¬ 
posed  to  the  implicit  Lyapunov  conditions).  These  features  will  be  exploited 
in  the  following  analysis. 

9 .3  Second-Order  System 

In  this  section,  we  shall  consider  the  effect  of  designing  a  stable 
feedback  for  a  second-order  plant  when  the  designer  assumes  a  model  of  order 
one.  This  simple  case  will  give  some  insight  into  the  more  general  case  to  be 
addressed  in  Section  9.4. 

Consider  the  stable  (2,1)  system 


y(k)  B  a^y(k  -  1)  +  a^y(k  -  2)  +  u(k) 


(9-14) 
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inode  Led  by  the  first-order  model  (obtained  through  LS  identification) 


y(k)  =  ajy(k  -  1)  +  u(k) 


(9-15) 


By  the  previous  discussion 


1  -  a 


and  |aj|  <  1,  | a ^ j  <  1  by  assumption,  Consider  the  feedback 


u(k)  =  -pa|y(k  -  1) 


A  stable  design  restricts  p  to  satisfy 


1  -  p|  < 


(9-16) 


However,  to  guarantee  closed-loop  stability  of  the  plant,  we  must  satisfy 


2  1 
a2  " 

1  -  a  l 


<  1 


or  equivalently 


,  2 
1  -  a2  -  p 


1  -  a; 


(9-17) 


The  constraints,  Eq.  (9-16)  and  (9-17),  are  plotted  in  Figure  9-2. 
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Figure  9-2. 


(C.1) 

(C.2) 


a 


2 

2 


O 

It  folLows  that  for  at  <  0,  any  control  law  that  satisfies  F;q.  (9-16)  is 

1  2 

stable  when  applied  to  the  full-order  plant.  If  a^  >  0,  there  exist  no  p  *0 
for  which  stability  can  be  guaranteed.  Note  that  if  Aj,  A2  are  the  plant 

eigenvalues,  =  -X^X^.  Thus,  a2  >  0  only  if  X^  X2  <  0,  which  is  possible 

only  if  both  eigenvalues  are  real  and  sign  (  X^ )  =  -sign(X2).  The  failure  to 

guarantee  stability  using  a  reduced  order  model  of  such  a  system  is  not  sur¬ 
prising.  We  note  that  the  result  above  is  also  applicable  in  the  case  of  a 

second-order  system  with  one  unstable  mode,  say  x| ,  and  where  a^  =  X^.  This 

discussion  is  summarized  in  Lemma  9-1. 


Lemma  9-1 


Let  Eq.  (9-14)  be  the  system  equation  and  Eq.  (9-15)  be  a  model  obtained 
from  Eq.  (9-14)  by  LS  definition.  Then,  the  feedback  u(k)  =  Ey(k  -  1),  satis- 

1  2 
fying  |a|  +  b|  <  1,  is  stable  if  <  0. 

2 

Conversely,  there  always  exists  a2  >  0  such  that  this  feedback  is 
unstable. 
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& 


t 


2 

Note  in  particular  that  for  -1  >  a^  <  0  and  jaj|  <1,  0  <  p  <  2  guaran¬ 
tees  stability,  and  that  | ( 1  -  p)j  <  |a!|  in  this  case.  We  shall  refer  to 

such  a  control  law  as  stability  enhancing  control  (see  Definition  9-2  in  Sec- 
t ion  9.4) . 

It  is  reasonable  to  assune  that  the  analysis  of  the  residual  energy  will 
provide  tighter  bounds  on  ja^j  thus  allowing  for  a  more  judicial  choice  of  p 
2 

to  accommodate  a^  >  0.  In  particular,  note  that  for  0  <  p  <  1,  stability  is 

2  1  +  O  ~  p)|aj| 

guaranteed  if  -1  <  a.  <  — — 

1  + 

9.4  n  Order  Systems 

In  this  section,  the  results  of  the  previous  section  are  generalized. 

To  enable  this  generalization,  a  few  definitions  and  relations  must  first  be 
presented. 

Definition  9-1 : 


S 


4 


fo 


S  is  called  the  shift  operator  since 


*0 

X1 

x. 

• 

1 

• 

s 

. 

• 

• 

X 

x  . 

n 

n-l_ 

(9-18) 


Notation 

The  vector  of  coefficients  in  the  polynomial  Pn(z)  of  Theorem  9-1  will 

be  denoted  by  o  ,  i.e.  , 

— n 
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T  h  n 

Ai  U’~ 


fl.-aj.---.-aJ 


(9-19) 


It  follows  that  the  matrices  and  A2  defined  in  Eq.  (9-13)  satisfy 


A.  =  'fa  Sa  Sna  ] 
1  L— n  — n  — n  J 


(9-20) 


A.  =  S  [jot  SJa  ...  SnJa  1 
2  L  — n  — n  — nJ 


(9-21) 


or  equivalent  ly,  from  Eq.  (9-13) 


Mn  -  l  ((s1”,)^)1  - 

i=0 


(9-22) 


By  Theorem  9-1,  P(z)  is  stable  if  and  only  if  M  >0,  i.e. 


Tu  r  r  t  T_l  i2  r  T„i+1  .  \2  -i 

x  M  x  =  )  lx  SaJ  -  x  S  Jo  1 

—  n—  .  ‘v.  L — n '  K  — n 1  J 

1=0 


>  0  -V-  xeR 


(9-23) 


Let  j2n+l  related  to  _an  via  Eq.  (9-10)  and  (9-11),  and  define 


Then  W"  obtain 


eRn+2  ,  xeRn+1 


iT»n+ii  -  [» -  (•£{  )2  i\i] 


/  -T  f-nl  n+l-T  f-nl  \2 


(9-24) 
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Thus,  M  .  >  0  if  M  >0  and  |an+*|  <  1  as  in  Theorem  9-1.  Finally,  the  fol- 
*  n+1  n  I  n+1 1 

lowing  concept  is  needed  for  the  generalization  of  Section  9.3. 

Let  Vn  =  Toep [vq, . . . ,vn+^ ]  be  the  solution  of  Eq.  (9-12)  corresponding 

to  a  .  Consider  the  stable  feedback  law 
— n 


and  define 


u(k) 


y(k  -  1) 
y(k  -  n) 


(9-25) 


a 

— n 


[1  .;  -(a"+  *x) 


and  correspondingly,  define  V  ,  P^Cz),  etc. 

Definition  9-2:  8  is  said  to  be  a  stability  enhancing  control  law  if  V  <  V 

— -  "  — n  — .  . «-  ■  ^  '•  ■  ■■  ■-  ■  n  —  n 

or  equivalently  Mn  M^. 

Using  the  concept  in  Definition  9-2,  the  main  result  of  this  study  is 
obtained. 

Theorem  9-2 

Let  a  and  a  be  related  via  Eq.  (9— 1 U)  and  (9-11),  Let  8  be  a  sta- 
— n+1  — n  — n 

bility  enhancing  control,  designed  for  the  system^.  Then  the  control  law, 
Eq.  (9-25),  is  stable  when  applied  to  the  system  if  it  satisfies  -V-  jkeR 


1(1  -  (vl)2X  *  Cl  R„]  5  >  0  (9-26) 


n+1 
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Remarks 


(1)  a  in  Eq.  (9-26)  can  be  replaced  by 
n 


(2)  The  matrix  Rn  is,  in  general,  nondefinite.  However,  the  sign 

definiteness  of  fl  -  1  i  M  +  an*)R  as  a  function  of  the 

v  n+1  '  J  n  n+1  n 

value  of  an+!  can  be  established  for  a  given  choice  of  stability 
n+1 

enhancing  control. 


2 

(3)  For  n  =  1,  it  is  easily  verified  that  Eq,  (9-26)  implies  -1  <  a^ 

<  0,  which  is  consistent  with  Lemma  9-1. 

(4)  Note  that  1  -  (an+|  f  M  +  an+jR  +  M  as  an+J  ^  0,  which  is  con- 

L  v  n+1  ’  J  n  n+1  n  n  n+1  * 

sistent  with  reasoning. 

(5)  The  proof  of  the  theorem  is  obtained  by  using  Eq.  (9-10),  (9-11), 
and  (9-12);  the  notion  of  stability  enhancing  control  and  using 
the  identities 


xVj 


a 

— n 


T  i 

x  S  Jet 
—  — n 


9.5  Conclusions 


The  effect  of  order  reduction  on  the  stability  of  feedback  design  has 
been  addressed.  The  assumption  made  was  that  the  reduction  of  order  is  based 
on  the  results  of  system  identification  rather  than  on  analytical  reduction  of 
known  full-order  systems. 

The  stability  analysis  was  based  on  the  relation  between  LS  estimation 
and  stability  theory.  Sufficient  conditions  for  a  "stability  enhancing 
controller"  to  be  stable  when  applied  to  a  full-order  system  have  been 
established. 

This  chapter  presents  only  a  first  stab  at  understanding  the  effect  of 
model  reduction  induced  by  identification  on  stability.  The  generalization  of 
Theorem  9-2  to  an  arbitrary  difference  in  orders  between  the  model  and  the 
plant  yields  similar  conditions  to  that  of  Eq.  (9-26).  However,  explicit 
constraints  on  the  ignored  parameters  are  cluttered  by  the  algebra.  Further 
studies  are  necessary  to  alleviate  this  algebraic  problem. 
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In  the  discussion,  it  was  assumed  that  the  ignored  parameter  a^j  is  ab¬ 
solutely  bounded  by  one.  Actually,  one  should  expect  that,  by  choosing  n  suf¬ 
ficiently  large,  tighter  bounds  on  a^*j  can  be  found  (e.g.  ,  by  relating  a^*j 

to  the  energy  in  the  residual  sequence).  Such  an  analysis  in  the  time  domain 
is  analogous  to  the  often  used  assunptlon  of  poor  knowledge  of  the  frequency 
response  of  a  plant  in  the  high-frequency  range  (see  e.g.,  [9-6]). 


The  results  presented  in  this  section  are  the  first  to  give  theoretical 
justification  to  a  common  practice,  success  of  simplified  modeling  has  been 
reported  in  various  fields  (e.g.,  power  generation,  which  is  very  similar  to 
LSS  in  its  characteristics  [9-7]).  It  is  our  intention  to  demonstrate  the 
practicality  of  this  approach  on  Draper  Model  # 2  in  the  near  future. 
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